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Summary 

This  manual  describes  the  use  of  and  relevant  equations  programmed  in  a  computer  code  designed 
to  carry  out  a  comprehensive  linear  analysis  of  multilayered  fiber  composites.  The  analysis  contains 
the  essential  features  required  to  effectively  design  structural  components  made  from  fiber 
composites.  The  program  is  an  outgrowth  of  two  in-house  computer  codes,  MFCA  (Multilayered 
Filamentary  Composite  Analysis)  and  INHYD  (Intraply  Hybrid  Composite  Design).  The  inputs  to 
the  code  are  constituent  material  properties,  factors  reflecting  the  fabrication  process,  and  composite 
geometry.  The  code  performs  micromechanics,  macromechanics,  and  laminate  analysis,  including 
the  hygrothermal  response  of  fiber  composites.  The  code  outputs  are  the  various  ply  and  composite 
properties,  composite  structural  response,  and  composite  stress  analysis  results  with  details  on 
failure.  The  code  is  in  Fortran  IV  and  can  be  used  efficiently  as  a  package  in  complex  structural 
analysis  programs.  The  input-output  format  is  described  extensively  through  the  use  of  a  sample 
problem.  The  code  manual  consists  of  two  parts.  The  mechanics  for  using  the  code  are  described  in 
the  first  part,  the  pertinent  equations  programmed  in  the  code  are  described  in  the  second  part. 


Introduction 

The  importance  of  and  need  for  a  multilevel  analysis  used  for  designing  structural  components 
made  of  multilayered  fiber  composites  are  documented  in  reference  1.  A  multilevel  analysis,  which 
was  efficient  in  predicting  the  structural  response  of  multilayered  fiber  composites  (with  the 
constituent  material  properties,  the  fabrication  process,  and  the  composite  geometry  known),  is  also 
documented  in  reference  1 . 

The  multilayered  analysis  presented  in  reference  1  consists  of  (1)  micromechanical  theories  for  the 
thermoelastic  properties  and  the  stress-level  limit  of  the  single  ply  as  functions  of  constituent  material 
properties  and  the  particular  fabrication  process,  (2)  the  combined  stress-strength  criterion  for  the 
single  ply,  and  (3)  multilayered  composite  structural  response  and  analysis  (macromechanical  or 
laminate  analyses),  where  the  interply  layer  effects  are  taken  into  account.  A  computer  code  designed 
to  carry  out  this  multilevel  analysis,  supplemented  as  noted  by  references  2  to  10,  has  been  developed 
at  the  Lewis  Research  Center.  This  code  is  identified  as  MFCA  for  Multilayered  Filamentary 
Composite  Analysis  (ref.  11). 

Intraply  hybrid  composites  are  a  logical  sequel  to  conventional  and  interply  hybrid  composites. 
Recently,  theoretical  and  experimental  investigations  have  been  conducted  on  the  mechanical 
behavior  of  intraply  hybrids  at  the  Lewis  Research  Center  (refs.  12  to  14).  The  theoretical  methods 
and  equations  described  in  these  references,  together  with  those  for  hygrothermal  effects  (ref.  15), 
have  been  integrated  into  a  computer  code  for  predicting  hygral,  thermal,  and  mechanical  properties 
of  intraply  hybrid  composites.  This  information  can  then  be  used  in  designing  these  composites.  This 
code  is  identified  as  INHYD  for  Intraply  Hybrid  Composite  Design  (refs.  16  and  17). 

The  present  computer  code  is  a  synergistic  combination  of  the  aforementioned  computer  programs 
MFCA  and  INHYD  together  with  several  significant  enhancements.  The  code  is  referred  to  as  ICAN 
for  Integrated  Composite  Analyzer.  It  utilizes  the  micromechanical  design  of  INHYD  and  the 
laminate  analysis  of  MFCA  to  build  a  comprehensive  analysis  capability  for  structural  composites. 
Additional  features  unique  to  ICAN  are  the  following: 


(1)  Ply  stress-strain  influence  coefficients 

(2)  Microstresses  and  microstress  influence  coefficients 

(3)  Stress  concentration  factors  around  a  circular  hole 

(4)  Calculation  of  probable  delamination  locations  around  a  circular  hole 

(5)  Poisson’s  ratio  mismatch  details  near  a  straight  free  edge 

(6)  Free-edge  stresses 

(7)  Material  cards  for  finite-element  analysis  using  NASTRAN  or  MARC 

(8)  Failure  loads,  summary  based  on  the  maximum  stress  criterion  and  laminate  failure  stresses, 
and  summary  based  on  first -ply  failure  and  fiber  breakage  criteria 

(9)  Transverse  shear  stresses  and  normal  stresses 

In  addition  to  the  above,  ICAN  has  its  own  data  base  of  material  properties  for  commonly  used 
fibers  and  matrices.  The  user  needs  to  specify  only  the  coded  names  for  the  constituents.  The 
program  searches  and  selects  the  appropriate  properties  from  its  library.  Furthermore,  input  data 
preparation  has  been  simplified  substantially  by  the  introduction  of  a  partial  free-field  format.  The 
output  formats  have  also  been  improved  significantly  to  ease  user  interpretation  of  the  results.  These 
enhancements  make  ICAN  significantly  more  user  friendly  than  its  predecessors.  The  computer  code 
has  been  programmed  in  Fortran  IV  and  has  been  tested  in  UNIVAC  1108,  IBM  370,  and  CRAY  1 
computers. 

Since  this  report  is  to  serve  as  a  users  manual,  the  code  is  divided  into  two  parts,  the  users  manual 
and  the  programmers  manual.  The  Users  Manual  describes  the  mechanics  of  using  the  code  with 
respect  to  program  format,  input  and  output,  and  sample  input  data  sets.  The  descriptions  are 
extensive  enough  sc  that  even  designers  and  analysts  with  little  or  no  programming  experience  can 
easily  use  the  code. 

The  programmers  manual  gives  the  various  subroutine  descriptions  and  the  equations  programmed 
therein,  with  details  on  the  input  and  output  and  the  global  storage  locations.  This,  along  with  the 
listing  of  the  source  program,  allows  the  user  to  make  his  own  modifications  to  the  code  as  they 
become  appropriate  for  further  enhancements. 

The  Fortran  variables  are  defined  in  appendix  A.  Included  is  information  such  as  which  part  of  the 
program  of  each  global  variable  is  generated.  Table  I  provides  a  summary  of  details  for  preparing 
data  cards,  and  the  input  data  given  in  table  II  provide  for  immediate  testing  of  the  code.  Properties 
for  a  few  commonly  used  fibers  and  matrix  materials  are  listed  in  appendix  B.  Appendix  C  shows 
sample  input  and  output  data  for  a  specific  case. 


Symbols 


^cx 

composite  axial  stiffness 

Ar 

™cx 

reduced  axial  stiffness 

BIDE 

boolean,  true  if  interply  effects  are  included 

C cx 

composite  coupling  stiffness 

Ce  t 

string  with  force  variables 

ce2 

string  with  displacement  variables 

COMSAT 

boolean,  true  if  laminate  analysis  is  wanted 

CSANB 

boolean,  true  if  membrane  and  axial  symmetry  exists 

DC,D( 

moisture  diffusivity 

composite  flexural  rigidity 

dr 

^cx 

reduced  bending  rigidity 

Dv 

displacement  vector 

df 

filament  equivalent  diameter 

filament  elastic  constants 

Ef\b  etc* 

fiber  normal  modulus 

ply  elastic  constants 

etc. 


Em,  Ecm 
Em  11,  etc. 

&mtetc 

F 

0/12  >  etc. 
Gm>  etc- 
Gmi2,  etc. 

Hj 

Hkc 

K 

U 

^cll,c22,c33 
^cxx,  cyy,  cxy 
KfU 

^/«1  1 
kv 

kf 

kfl 


kmV 

kvl 

Lsc 

Mz 

Mcx 

McM(x 

McT(x 

m 


Nr. 


N, 


cM(x 


N, 


cT(x 


Nf 

Nf 

Nfc 

NONUDF 


N, 


pc 


N, 


'pl 


cp 


nP 

Qf,i,p,r,s 

R 

RINDV 


ply  normal  modulus 

matrix  elastic  constants 

matrix  normal  modulus 

matrix  failure  strain  allowables 

combined  stress-failure  function 

fiber  shear  modulus 

ply  shear  modulus 

matrix  shear  modulus 

interply  distortion  energy  coefficient 

array  of  constituent  heat  conductivities 

composite  heat  capacity 

index,  generally  ply  or  interply 

composite  three-dimensional  heat  conductivities 

composite  two-dimensional  heat  conductivities 

fiber  heat  conductivity 

ply  heat  conductivity 

matrix  heat  conductivity 

apparent  void  volume  ratio 

actual  fiber  volume  ratio 

ply  apparent  fiber  volume  ratio 

actual  matrix  volume  ratio 

ply  apparent  matrix  volume  ratio 

ply  apparent  void  volume  ratio 

array  of  limiting  conditions 

ply  moisture 

applied  moment 

hygral  moment 

thermal  moment 

load  condition  index 

applied  membrane  loads 

hygral  force 

thermal  force 

number  of  fibers  per  end 

number  of  plies 

number  of  load  conditions 

boolean;  true  if  detailed  Poisson’s  ratio  differences  chart  is  to  be  suppressed 

string  PROPC  length 

string  PROP  length 

composite  properties  array 

string  PROPC 

ply  properties  array 

string  PROP  main  program 

indices  to  print  out  string  PROP 

transformation  matrix 

boolean;  true  if  displacements  are  known 


Sc 

Sfii  T>  etc- 

T( 

h 

wcb 

x,y,z 


ac 

af 

a( 

am 

0C 

0e>0e 

0h 

0(,0f,0  m 

0S 

0V 

5/ 

6f 

ecs 

ecsx 

ee 

v/\ 2,  etc. 
vm>  etc- 

pm\2 »  etc* 

Pf,m 

Pmw 

1,  2,  3 


composite  failure  stress 
ply  limit  failure  stresses 
ply  temperature 
ply  thickness 

composite  local  curvature  changes 

structural  reference  axes 

composite  coefficient  of  thermal  expansion 

fiber  thermal  coefficient  of  expansion 

ply  thermal  coefficient  of  expansion 

matrix  thermal  coefficient  of  expansion 

moisture  expansion  coefficients 

correlation  factors  for  ply  thermoelastic  properties 

correlation  factor  for  ply  heat  conductivity 

moisture  expansion  coefficients  for  ply,  fiber,  and  matrix 

correlation  factor  for  ply  strength 

matrix  strain  magnification  due  to  ply  strain  in  the  presence  of  voids 
interfiber  spacing 
interply  layer  thickness 
interfiber  spacing 

angle  between  composite  material  and  structural  axes 
reference  plane  membrane  strain 
ply  strain 

angle  between  ply  material  and  composite  axes 

fiber  Poisson’s  ratio 

ply  Poisson’s  ratio 

matrix  Poisson’s  ratio 

fiber  and  matrix  weight  density 

density  of  matrix  with  moisture 

ply  stresses,  fiber  stresses,  and  matrix  stresses 

material  reference  axes 


Users  Manual 

The  mechanics  required  to  use  this  code  for  the  analysis  of  multilayered  fiber  composites  are 
described  in  this  part  of  the  report.  The  theory  on  which  the  code  is  based  is  described  in  the  second 
part  of  the  report  (Programmers  Manual). 

The  physical  representations  of  the  constituents  used  in  the  code  are  illustrated  in  figure  1 .  This 
figure  shows  a  complete  integration  schematic  starting  with  the  constituent  materials,  fiber  and 
matrix.  The  required  input  properties  and  computed  properties  at  various  levels  are  summarized  in 
symbolic  form  as  follows: 

(1)  Properties  required  by  code  as  input  for  a  fiber:  £/n(22,335  ^/I2,23,l3i  G/12,22,135  «/ll,22,33» 
^/l  1,22, 33 >  Hcf\  pp  Np  df,  and  Sft. 

(2)  Properties  required  by  code  as  input  for  a  matrix:  Em 1U2,335  ^12,23,135  G/12, 23,13  amll, 22,335 

11,22,335  Prri>  ^mc>  ^ mpV>  ^ mpc 5  ® mps 5  an<^  &mpTOR- 

(3)  Properties  required  by  code  as  input  for  a  single  ply:  fiber  and  matrix  properties  and  ply 
characteristics  /3e ,  0n ,  0S9  and  7). 


(c)  Single  ply. 


(d)  Composite. 

Figure  1. — Schematic  of  typical  multilayered  fiber  composite  and  some  basic  components. 


(4)  Properties  computed  by  code  for  single  ply:  Em> 22,33;  "02,23,13;  Gfl2, 23,131  “01,22,331  ,22,331 

Hct\  pf,  te;  8t;  Sflir,iic,22r,22C,l25,23Sl  ^fl2l  an^  stress  analysis  factors  1,22,121  ^o  1,22,12!  an<3 
1.0 -F{a,S,K(n). 

(5)  Properties  required  by  code  as  input  for  a  composite:  ply  properties  and  composite 
characteristics  0#,  Hp  2a(3>  Ncx>  Mcx  or  Ucx ,  and  Wcx. 

(6)  Output  computed  by  code  for  a  composite:  (ecJ  =  [iic](<jc)  +  [Ec]  ~ KCXXyyjXy;  Hc ; 


i 

A  C 
^cx^cx 

I 

1 

Vex 

1  ! 

1  +  < 

N cxTf 

\MCX\ 

u 

C  D 

'-'CXr'CX 

1 

Wcx 

1  1 

McxTt 

and  the  inverse  A<pdelj. 

Figure  2(a)  shows  the  subroutines  and  sequence  in  the  code.  The  subroutines  between  the  Main 
program  and  the  input  data  may  be  arranged  in  any  desired  order.  The  user  should  refer  to  figure  2(b) 
for  the  logic  flow  of  the  analysis. 

The  following  four  steps  are  required  to  use  the  code  in  the  user’s  computer  facility: 

(1)  Obtain  the  code 

(2)  Make  it  operational  in  the  user’s  computer  facility 

(3)  Supply  the  input  data 

(4)  Interpret  the  code  output  results 

Obtain  the  Code 

The  code  may  be  obtained  in  cards.  If  this  is  not  convenient  or  possible,  the  cards  can  be  punched 
from  the  compiled  listing  (contact  COSMIC,  The  University  of  Georgia,  Athens,  GA  30602, 
concerning  the  availability  of  this  program). 


5 


JCL  cards 

Main  program,  ICAN 

2 

Micromechanics,  INHYD 

5 

Composite  properties,  FIBMT 

4 

Hygrothermomechanical 
response,  HTM 

i 

|5 

— U 

Hybrid  composite  properties, 
COMPP 


Maximum  of  an  array,  AMAXF 


Three-dimensional  composite 
properties,  GACD3 


Two-dimensional  composite 
properties,  GPCFD2 


Flexural  strengths,  FLFXX 


Composite  stress  analysis, 
COMSA 


10 


Failure  mode  and  failure  load, 
FLBLD 


20 


Microstress  influence 
coefficients,  MINCOF 


Name  strings,  BLOCK  DATA 


22 


Print  the  emblem  ICAN,  LOGO 


23 


11 

User  input,  Input  data 

24 

Read  constituent  properties, 

BANKRD 

unit  5 

Calculate  free  edge  stresses, 

12 

Database,  FBMTDATA.  BANK 

EDGSTR 

unit  8 

Property  generator  for 
INHYD,  IDGER 


(a)  Schematic  showing  subroutines  and  sequence  of  ICAN  code. 
Figure  2. — Subroutines,  sequence,  and  logic  flaw  of  ICAN  code. 


Make  it  Operational 

A  prerequisite  to  the  program  is  the  availability  of  a  Fortran  compiler  in  the  user’s  computer 
facility.  To  run  the  program,  certain  computer-system-dependent  control  cards  (Job  Control 
Language  (JCL)  cards)  may  also  be  necessary.  The  computer  system  personnel  should  be  consulted 
about  these  items. 

Once  the  deck  of  cards  has  been  assembled  (except  input  data)  with  the  proper  control  cards  as 
shown  in  figure  2,  the  user  is  ready  to  compile  the  code  in  his  facility.  The  compilation  will  indicate 
whether  any  additional  modifications  are  needed.  Most  modifications  will  be  minor  and  will  usually 
deal  with  certain  Fortran  statements  peculiar  to  each  compiler. 

Supply  the  Input  Data 

The  physical  arrangement  of  the  input  data  cards  is  illustrated  in  figure  3.  Details  for  preparing  the 
input  data  cards  are  summarized  in  table  I.  A  detailed  description  of  these  cards  is  given 
subsequently.  A  sample  for  preparing  input  data  for  a  four-ply  symmetric  laminate  is  presented  in 
table  II.  This  laminate  has  two  different  material  systems.  The  0°  plies  are  of  AS  graphite 
fiber/intermediate-modulus,  low-strength  epoxy  matrix  composite.  The  90°  plies  are  made  of  a 
hybrid  composite.  The  primary  composite  is  S  glass/high-modulus,  high-strength  epoxy.  The 
secondary  composite  is  AS  graphite/intermediate-modulus,  high-strength  epoxy. 


ICAN  (Main)  read  and  write 

Emblem  printer,  LOGO 

input  data. 

Call  various  subroutines. 

* 

Select  fiber  and  matrix 
properties  from  the  data 

Generate  input  data  deck 
for  INHYD 

bank,  BANRD 

Microstresses,  MCRSTR 


X 


^  Read  input  from  unit  7. 
Call  various  subroutines  to 
compute  ply  level  properties 
and  return  to 
ICAN  (Main). 


Three-dimensional  properties, 
GACD3 


Free  edge  stresses,  EDGSTR 


Two-dimensional  properties, 
GPCFD2 

. . i  ~ 


Composite  stresses,  COMSA 


Microstress  influence 

_ ^ 

Stress  concentration  factors, 

Poisson's  ratio  differences, 

coefficients,  MINCOF 

STRCNF 

NUDIFS 

Failure  load  analysis,  MSCBFL 


Compute  influence  coefficients. 
Call  COMSA 


(b)  Schematic  showing  logic  flow  of  ICAN  code. 
Figure  2.— Concluded. 


Figure  3. — Physical  arrangement  of  input  data  cards. 


TABLE  I.— SUMMARY  OF  DETAILS  FOR  PREPARING  INPUT  DATA  CARDS 


Card 

group 

Identification 

Code  symbol 

Number  of 
entries 

List  of 
entries, 
sequential 
order 

Card 

field 

columns 

Format 

Comments  and 
engineering  units 

1 

Title  card 

TITLE 

80 

Alphabetic 

characters 

1  to  80 

(10a8) 

2 

STDATA 

NL,NLC,NMS 

3 

Nf,NfC,Nms 

9  to  36 

(a8,3I8) 

Composite  geometries 

3 

Boolean  for  input 
displacement 
Boolean  for  interply 
layer  energy 
contribution 

Boolean  for  mem¬ 
brane  and  bending 
symmetry 

Boolean  for  laminate 
analysis 

Boolean  for  Poisson 
ratio  chart 

RINDV 

BIDE 

CSANB 

COMSAT 

NONUDF 

1 

1 

1 

1 

1 

I  to  6 

1  to  6 

1  to  6 

1  to  6 

1  to  6 

T  (true)  if  displays  are  inputs; 

otherwise  F  (false) 

T  (true)  if  contributions  are 
desired;  otherwise  F  (false) 

T  (true)  if  symmetry  exists; 
otherwise  F  (false) 

T  (true)  if  laminate  analysis  is 
desired;  otherwise  F  (false) 

T  (true)  if  Poisson  ratio  differ¬ 
ences  chart  is  not  desired; 
otherwise  F  (false) 

(L6) 

4 

PLY  (ply  desired) 

INP1,IP1,TU,TCU 

DELM,TETA,THCKNS 

7 

TCU,AM, 

dt,t( 

1  to  64 

(a8,2I8,5F8.3) 

Ply  layup  and  temperature 
and  moisture  conditions 

5 

MATCRD  (material 
system  details) 

CODES(l,U), 

CODES(l,2,J), 

VFP,VVP, 

CODES(2,l,J), 

CODES(2,2,J), 

VFS,VVS,VSC 

9 

Primary 

composite 

code 

names  for 
fiber  and 
matrix,  kp 
kv; 

Secondary 

composite 

code 

names 

for  fiber 

and 

matrix, 
ksr,  kp,  kv 

1  to  64 

(a8,2a4,2E8.2) 

Description  of  material 
systems  to  be  used 

6 

PLOAD  (loading 
details) 

NX,NY,NXY,THCS 

MX,MY,MXY 

DMX,DMY 

PRSSU,PRSSL 

4 

3 

4 

6cs,Mx,My, 

Mty 

dM,  Id., 
dMy  Idy, 

Puf,  ' 

1  to  32 

1  to  32 

1  to  40 

(a8,7E8.4) 

(a8,7E8.4) 

(a8,7E8.4) 

Loading  conditions  (inplane) 
Angle  of  inclination  to  the 
structural  jr-axis 

Loading  conditions  (bending) 
Loading  conditions 
(tranverse) 

Input  data  for  additional  composite  systems  may  be  easily  prepared.  This  is  done  by  selecting  the 
desired  fiber  and  matrix  from  the  available  materials  listed  in  appendix  C  using  the  variable 
FBMTDATA.BANK  and  modifying  the  appropriate  entries  in  the  input  data  sample  illustration. 

After  the  input  data  have  been  properly  assembled  (as  shown  in  fig.  3),  they  are  placed  in  their 
physical  position  (fig.  2)  and  the  code  is  ready  to  be  run. 

Detailed  Description  of  Input  Data 

The  card  group  numbers  referred  to  here  are  given  in  figure  3  and  table  I.  The  sequential  order  of 
the  entries  in  each  card  group  is  given  in  table  I.  Note  that  most  data  cards  are  identified  by  a 
mnemonic  to  indicate  the  card  group  in  which  it  belongs  in  the  input  data  deck.  Also,  most  data  cards 
are  divided  into  fields  of  eight,  with  one  entry  per  field  being  allowed.  The  mnemonic  is  entered  in 


TABLE  II.— ICAN  SAMPLE  INPUT  DATASET 


FOUR  PLY  SYMMETRIC  LAMINATE. 

ICAN  SAMPLE  INPUT 

DATA. 

STDATA 

4 

1 

2 

T 

COMSAT 

F 

CSANB 

F 

BIDE 

F 

RINDV 

T 

NONUDF 

PLY 

1 

1 

70.00 

70.0 

0.0 

0.0 

0.010 

PLY 

2 

2 

70.00 

70.0 

.0 

90.0 

.005 

PLY 

3 

2 

70.00 

70.0 

.0 

90.0 

.005 

PLY 

4 

1 

70.00 

70.0 

.0 

.0 

.010 

MATCRDAS— IMLS 

0.55 

0.02 

AS— IMLS 

0.0 

0.57 

0.03 

MATCRDSGLAHMHS 

.55 

.01 

AS— IMHS 

.4 

.57 

.01 

PLOAD  1000. 

0.0 

0.0 

0.0 

NX,NY,NXY,THCS 

PL0AD  0.0 

.0 

.0 

MX,MY,MXY 

PLOAD  .0 

.0 

DMX/QX, DMY/QY , PRSS 

format  A8,  and  the  integers  are  entered  in  format  18.  The  real  numbers  may  be  entered  anywhere  in 
the  appropriate  field.  The  following  is  a  brief  description  of  each  card  group  together  with  examples 
taken  from  table  II: 

Title  card. 


Any  title  of  length  up  to  80  characters 

FOUR  PLY  SYMMETRIC  LAMINATE.  ICAN  SAMPLE  INPUT  DATA. 


As  shown,  any  title  of  length  up  to  80  characters  including  blanks  may  be  supplied  on  this  card. 
Starting  data  card . 


1  8,9  16,17  24,25  a32 


Mnemonic 

N( 

Nic 

^ ms 

STDATA 

4 

1 

2 

This  card  has  a  mnemonic  STDATA.  It  contains  the  overall  laminate  and  loading  details.  These 
details  are  the  number  of  plies  Afy,  the  number  of  loading  conditions  N&,  and  the  number  of  different 
material  systems  Nms. 

Booleans . 


1  6,7 


Boolean 

T  or  F 

This  space  may  be  used  for  comments 

T 

COMSAT 

F 

RINDV 

F 

BIDE 

F 

CSANB 

T 

NONUDF 

A  set  of  booleans,  COMSAT,  RINDV,  BIDE,  CSANB,  and  NONUDF  is  defined  through  these 
cards.  These  are  five  cards,  one  per  each  logical  variable.  The  format  is  L6.  The  variables  have  the 
following  functions: 

(a)  COMSAT.— The  letter  T  in  the  card  will  direct  the  program  to  perform  a  complete  laminate 
analysis.  A  letter  F  would  terminate  the  program  at  this  stage. 


(b)  RINDV.— The  letter  T  is  entered  in  the  card  if  the  displacements  are  inputs;  otherwise,  the 
letter  F  is  entered. 

(c)  BIDE.— The  letter  T  is  entered  in  the  card  if  the  interply  layer  contributions  on  the  composite 
are  desired;  otherwise,  the  letter  F  is  entered. 

(d)  CSANB.— The  letter  T  is  entered  in  the  card  if  the  composite  has  both  membrane  and  bending 
symmetry;  otherwise,  the  letter  F  is  entered. 

(e)  NONUDF.— The  letter  T  is  entered  if  the  detailed  Poisson’s  ratio  difference  chart  is  to  be 
suppressed;  otherwise,  the  letter  F  is  entered. 

Ply  details  card  group. 


1  8,9  16.17  24.25  32,33  40,41  48,49  56,57  64 


Mnemonic 

Ply 

Material 

MID. 

T 

1u 

T 

1  cu 

M 

e( 

h 

PLY 

1 

1 

70.00 

70.0 

.0 

0.0 

.015 

PLY 

2 

2 

70.00 

70.0 

.0 

90.0 

.005 

PLY 

3 

2 

70.00 

70.0 

.0 

90.0 

.005 

PLY 

4 

1 

70.00 

70.0 

.0 

0.0 

.010 

All  the  cards  in  this  group  have  the  mnemonic  PLY.  The  number  of  cards  is  Afy,  with  eight  entries  on 
each  card.  The  first  entry  is  PLY.  The  second  and  third  entries  are  identification  numbers  for  the  ply 
and  the  material  system,  respectively.  The  fourth  and  fifth  entries  are  the  use  temperature  Tu  and  the 
cure  temperature  Tcu,  respectively.  The  sixth  entry  is  the  percentage  of  moisture  M.  The  seventh  and 
the  eighth  entries  are  the  orientation  angle  0  of  the  ply  and  the  thickness  of  the  ply,  respectively.  A 
default  value  of  0.005  is  taken  for  the  thickness  if  this  entry  is  missing.  The  material  system 
identification  number  should  be  different  not  only  for  different  composite  systems  but  also  for 
varying  use  temperature  or  moisture  content  from  ply  to  ply. 

Material  system  details . 


1  8,9  16,17  24,25  32,33  40,41  48,49  56,57  64 


Mnemonic 

Fiber, 

matrix 

kf 

K 

Fiber, 

matrix 

V 

r  sc 

kf 

K 

MATCRD 

AS— IMLS 

.55 

.02 

AS— IMLS 

0.0 

.57 

.03 

MATCRD 

SGLAHMHS 

.55 

.01 

AS— IMHS 

0.4 

.57 

.01 

All  the  cards  in  this  group  have  the  mnemonic  MATCRD.  The  number  of  cards  is  Nms  with  10  entries 
in  each  card.  The  first  entry  is  MATCRD.  The  second  and  the  third  entries  are  coded  words  for  fiber 
and  matrix  material  of  the  primary  composite.  The  code  words  are  entered  in  format  2A4.  For 
example,  the  code  for  AS-type  fiber  is  AS—  and  epoxy  matrix  is  EPOX.  A  dictionary  of  codes  for 
several  fibers  and  matrices  is  provided  in  appendix  C.  The  user  may  choose  any  combination  of  fiber 
and  matrix  for  a  composite  system.  The  fourth  and  the  fifth  entries  pertain  to  the  details  of  the 
primary  composite  system.  They  are  the  primary  fiber  volume  ratio  and  the  primary  void  volume 
ratio,  respectively.  The  next  two  entries  refer  to  the  secondary  composite  system  which  is  applicable 
for  the  case  of  the  hybrid  composite  ply.  They  should  be  the  same  as  the  second  and  third  entries  for 
standard  composite  systems.  The  next  entry  is  the  secondary  composite  system  volume  ratio.  This  is 
zero  for  the  standard  composite  systems.  The  last  two  entries  are  the  fiber  volume  ratio  and  the  void 
volume  ratio  for  the  secondary  composite  system.  These  values  are  entered  when  applicable. 

Load  cards . 


1  8,9  16,17  24,25  32,33  40 


Mnemonic 

Nx 

Nit 

Tkc 

PLOAD 

1000. 

0.0 

0.0 

0.0 

PLOAD 

0.0 

0.0 

0.0 

PLOAD 

0.0 

0.0 

0.0 

0.0 

All  the  cards  in  this  group  start  with  the  mnemonic  PLOAD.  There  are  three  cards  for  each  loading 
condition.  Thus,  the  total  number  of  cards  is  3A^.  The  first  card  under  each  loading  condition 
contains  entries  Nx,  Ny ,  and  Nxy  for  the  membrane  loads  and  Thcs  for  the  orientation  of  the  loads 
with  respect  to  the  structural  axes.  Similarly  the  second  card  contains  the  bending  resultants  Mx,  My, 
and  Mxy.  The  last  card  contains  the  transverse  shear  resultants  DMX  and  DMy  and  the  transverse 
pressures  Pu  and  P 

The  user  input  data  are  read  from  I/O  unit  5.  Apart  from  this,  ICAN  uses  two  more  units,  7  and  8, 
for  its  I/O  operations.  Unit  8  is  used  to  store  the  material  properties  data  base.  Unit  7  is  used  as  a 
scratch  file  by  ICAN.  These  I/O  units  must  be  appropriately  defined  by  using  the  operating  system 
JCL. 

Output 

The  following  items  are  printed  out  by  the  program: 

(1)  ICAN  logo 

(2)  ICAN  coordinate  systems 

(3)  ICAN  input  data  echo 

(4)  Input  data  summary 

(5)  Fiber,  matrix,  and  ply  level  properties  of  primary  and  secondary  composites 

(6)  Composite  three-dimensional  strain-stress  and  stress-strain  relations  about  the  structural  axes; 
MAT9  card  for  MSC/NASTRAN  solid  elements 

(7)  Composite  properties  generated  in  array  PC 

(8)  Composite  constitutive  equations  about  the  structural  axes 

(9)  Reduced  bending  and  axial  stiffnesses 

(10)  Data  for  finite-element  analysis 

(11)  Displacement-force  relations  for  the  current  load  condition 

(12)  Ply  hygrothermomechanical  properties/response 

(13)  Details  of  Poisson’s  ratio  mismatch  among  the  plies 

(14)  Free  edge  stresses 

(15)  Microstresses  and  microstress  influence  coefficients  for  each  different  composite  material 
system 

(16)  Stress  concentration  factors  around  a  circular  hole 

(17)  Locations  of  probable  delamination  around  circular  holes 

(18)  Ply  stress  and  strain  influence  coefficients 

(19)  Laminate  failure  load  analysis  based  on  the  first-ply  failure/maximum  stress  criteria 

(20)  Summary  of  the  laminate  failure  stress  analysis  based  on  two  alternatives,  first-ply  failure  and 
fiber  breakage 

The  printout  of  the  input  data  summary  (app.  B  item  4)  shows  details  regarding  composite 
geometry,  constituent  specifications,  temperature  and  moisture  conditions,  and  the  loading 
conditions. 

The  next  few  pages  of  the  output  are  generated  by  the  INHYD  program  package.  They  show  the 
fiber -matrix  properties  for  the  different  composite  systems  and  the  ply  level  properties  of  the 
composites  (app.  B,  item  5). 

The  output  of  the  composite  three-dimensional  strain-stress  temperature  and  moisture  relations 
and  composite  stress-strain  relations  about  the  structural  axes  are  printed  under  the  following 
headings: 

(a)  3-D  COMPOSITE  STRAIN  STRESS  TEMPERATURE  MOISTURE  RELATIONS- 
STRUCTURAL  AXES 

The  matrices  and  {(3C}S  in  the  equation 

fecJ s  ~  [Ec) s  *  facls  ”  A  T$ (ac]s  —  Mglflcls 


where  ATg=  T?-  Tcu 


are  printed  by  the  subroutine  GACD3. 


(b)  3-D  COMPOSITE  STRESS  STRAIN  RELATIONS-STRUCTURAL  AXES 
The  matrix  [isjj  in  the  equation 


is  printed  out  by  the  subroutine  GACD3. 

The  subscripts  in  the  preceding  equations  indicate  that  the  relations  are  written  about  the  structural 
axes.  It  is  noted  that  these  properties  are  only  local  to  subroutine  GACD3.  They  can  be  made  global 
if  needed.  The  properties  needed  to  prepare  the  MAT9  card  of  MSC/NASTRAN  are  printed  out  next 
under  the  heading  MAT9  CARD  FOR  MSC/NASTRAN  SOLID  ELEMENTS  (app.  B,  item  6). 

The  output  of  the  composite  properties,  generated  in  array  PC,  are  printed  under  the  following 
heading  (app.  B,  item  7): 

COMPOSITE  PROPERTIES— VALID  ONLY  FOR  CONSTANT  TEMPERATURE  THROUGH 
THICKNESS  LINES  1  to  31:  3-D  COMPOSITE  PROPERTIES  ABOUT  MATERIAL  AXES 
LINES  33  to  62:  2-D  COMPOSITE  PROPERTIES  ABOUT  STRUCTURAL  AXES 
Sixty-two  entries  are  printed  under  this  heading  as  shown  in  the  following  list: 


Code  name 

Notation 

Explanation 

PC(D 

Pc 

weight  density 

PC(2) 

A 

thickness 

PC(3)  to  PC(ll) 

[Ec\ 

three-dimensional  stress-strain  relations  about  material  axes 

PC(12)  to  PC(14) 

K) 

three-dimensional  coefficients  of  expansion  about  material 
axes 

PC(15)  to  PC(18) 

[KC},HC 

three-dimensional  heat  conductivity  and  heat  capacity  along 
material  axes 

PC(19)  to  PC(30) 

EcU>Gc\2’vcn 

three-dimensional  constants  about  material  axes 

PC(31) 

Zc 

distance  to  reference  plane  from  bottom  of  composite 

PC(32) 

blank 

PC(33)  to  PC(38) 

[Ec\~l 

two-dimensional  stress-strain  relations  about  structural  axes 

PC(39)  to  PC(47) 

Eq  1 GcYbvc\2 

two-dimensional  elastic  constants  along  structural  axes 

PC(48)  to  PC(54) 

exp  Kg,  Hq 

two-dimensional  coefficients  of  thermal  expansion,  heat 
conductivity,  and  heat  capacity  along  structural  axes 

PC(55)  to  PC(58) 

Dc 

moisture  diffusivities 

PC(59)  to  PC(62) 

He 

moisture  expansion  coefficients 

Array  PC  and  its  corresponding  string  and  headings  are  controlled  by  the  formats  in  subroutine 
GPCFD2.  For  nonuniform  temperature/moisture,  the  bending  equivalent  and  the  membrane 
equivalent  elastic  constants  may  be  obtained  by  utilizing  the  reduced  bending  rigidity  matrix  and  the 
reduced  stiffness  matrix  which  are  also  regular  output  of  ICAN. 

The  output  for  the  composite  constitutive  equations  are  printed  in  the  following  manner  (app.  B, 
item  8): 


FORCES  FORCE  DISPLACEMENT  DISPL  T-FORCES  H-FORCES 

RELATIONS 


( (Ay  j 

[A  cx\  i^cxl 

\ECsx\ 

(  {NCTtx\ 

[HcM(X) 

(  [Mcx]  ) 

[Qd  1 'Pcx\ 

1 

{wcb}\ 

|  (  ] 

!  lMCMtx\ 

The  elements  of  matrices  Acx ,  Cex,  Dcx ,  Ncj(X3  NcM(X,  McT?x,  and  McM(X  are  printed  out  by  the 
subroutine  GPCFD2. 

The  output  for  the  reduced  bending  rigidities  is  printed  under  the  heading  (app.  B,  item  9): 
REDUCED  BENDING  RIGIDITIES.  The  elements  of  fD^J  are  printed  out  as  a  matrix. 

Similarly,  the  output  for  the  reduced  axial  stiffness  [A^x]  is  printed  out  under  the  heading 
REDUCED  STIFFNESS  MATRIX.  The  corresponding  formats  for  the  above  two  outputs  are  in 
subroutine  GPCFD2  (app.  B,  item  10). 


The  next  printout  comes  from  the  main  program  under  the  heading:  SOME  USEFUL  DATA  FOR 
F.E.  ANALYSIS.  This  information  is  useful  for  preparing  material  data  cards  for  finite  element 
codes  NASTRAN  and  MARC. 

The  inverse  of  the  constitutive  equations  is  printed  out  in  the  following  manner  (app.  B,  item  11): 


DISP  DISPLACEMENT  FORCE  FORCES 
RELATIONS 


l^caxl 

-1 

[NCx\  \.NcT,x)  lNcMtxl 

+  1  +  1 

l^cb)  ] 

1 

[^cd  [®cd 

1 

{Mcxl  [McT(X !  \mcM(X\ 

The  elements  of  this  inverse  are  printed  out  in  the  subroutine  COMSA. 

The  current  values  for  the  loads  and  corresponding  set  of  ply  properties  generated  in  array  PL  are 
printed  out  next  (app.  B,  item  12).  The  explanations  of  the  75  entries  in  the  PL  property  array  are 
given  in  the  following  list: 


Code 

name 

Notation 

Explanation 

PL(1,I) 

kv 

ply  void  volume  ratio 

PL(2,I) 

kfi 

ply  apparent  fiber  volume  ratio 

PL(3,I) 

kf 

ply  actual  fiber  volume  ratio 

PL(4,I) 

kmi 

ply  apparent  matrix  volume  ratio 

PL(5,I) 

km 

ply  actual  matrix  volume  ratio 

PL(6,I) 

Pe 

ply  weight  density 

PL(7,I) 

h 

ply  layer  thickness 

PL(8,I) 

ply  and  interply  layer  thickness 

PL(9,I) 

Hj 

interply  layer  distortion  energy  coefficient 

PL(10,I) 

Z( 

distance  from  bottom  of  composite  to  ply  centroid 

PL(11,I) 

ZCg 

distance  from  reference  plane  to  ply  centroid 

PL(12,I) 

<>cs 

angle  from  structural  axes  to  composite  material  axes 
(same  for  all  plies)  (fig.  2) 

PL(13,I) 

angle  from  ply  material  axes  to  composite  material  axes 
(fig.  2) 

PL(14,I) 

angle  from  ply  material  axes  to  composite  structural  axes 
(fig.  2) 

PL(15,I)  to 

PL(23,I) 

ply  stress-strain  relations 

PL(24,I)  to 

PL(26,I) 

ply  thermal  coefficients  of  expansion 

PL(27,I)  to 

PL(29,I) 

m 

ply  heat  conductivities 

PL(30,I) 

Hd 

ply  heat  capacity 

PL(34,I)  to 

PL(42,I) 

Em<vt\2’Gn2 

ply  elastic  constants 

PL(43,I)  to 

PL(48,I) 

De  and  /3f 

moisture  diffusivities  and  expansion  coefficients 

PL(49,I) 

Pfi  del 

interply  delamination  factor 

PL(50,I) 

Tt 

ply  temperature 

PL(51,I)  to 

PL(60,I) 

Smr>  etc. 

ply  limiting  stresses 

PL(61,I) 

K(\2  ap 

coefficient  in  combined  stress-strength  criterion 

PL(62,I) 

combined  stress-strength  criterion 

PL(63,I) 

interply  delamination  criterion 

PL(64,I)  to 

PL(69,I) 

ply  applied  strains  and  stresses 

PL(70,I) 

Apj 

adjacent  ply  relative  rotation 

PL(71,I) 

Hoffman's  failure  criterion 

PL(72,I) 

M( 

ply  moisture 

PL(73,1) 

fff13 

transverse  shear  stress 

PL(74,I) 

<7I’23 

transverse  shear  stress 

PL(75,1) 

ffl’33 

thickness  stretch  stress 

The  next  printout  shows  Poisson’s  ratio  differences  between  the  plies  and  the  composite  (app.  B, 
item  13).  They  are  printed  out  by  the  subroutine  FESTRE  under  the  heading  DETAILS  OF 
POISSON’S  RATIO  MISMATCH. 

The  stress  peaks  near  the  free  edge  region  are  printed  out  next  by  the  subroutine  EDGSTR  under 
the  heading  (app.  C,  item  14)  FREE  EDGE  STRESSES. 

Item  14  of  appendix  B  shows  ply  stresses  in  the  structural  coordinate  system  and  the  through-the- 
thickness  stresses  uzv  axz ,  and  ayz.  The  boundary  layer  decay  length  is  also  shown  in  the  table  under 
the  heading  YDCAY  LENGTH.  Care  must  be  exercised  in  interpreting  the  results.  They  are  based  on 
approximate  engineering  theories  and  give  good  qualitative  information  regarding  the  relative 
magnitudes  of  the  peaks  in  the  individual  plies.  This  printout  is  suppressed  in  the  case  of  combined 
loading. 

The  microstresses  in  each  ply  are  printed  out  next  by  the  subroutine  MCRSTR  (app.  B,  item  15(a)). 
Two  regions  of  interest  are  considered  for  the  computations,  the  region  between  the  fibers  composed 
entirely  of  matrix  (A)  and  the  region  consisting  of  fibers  as  well  as  matrix  (B).  The  stresses  are  given  a 
descriptive  notation.  Thus,  SM2AL  means  stress  in  matrix  along  the  transverse  (2)  direction  in  region 
A  due  to  a  ply  stress  along  the  longitudinal  direction  of  the  material.  Figure  4  shows  the  definitions 
for  regions  A  and  B.  The  printout  also  shows  microstresses  resulting  from  moisture  and  temperature 
differences  if  nontrivial  M%  and  7}  are  present. 

The  microstress  influence  coefficients,  stresses  due  to  unit  applied  stresses  in  direction  11,  22,  12, 
13,  and  23  (app.  B,  item  15(b));  unit  temperature  difference  Tf,  and  unit  moisture  content  Afyare 
output  from  the  subroutine  MINCOF.  These  variables  are  printed  out  following  the  microstresses. 

Under  the  heading  STRESS  CONCENTRATION  FACTORS  (app.  B,  item  16)  are  printed  out  the 
factors  Klxx,  K\yy ,  and  KXxy  which  are  due  to  inplane  loading  around  a  circular  hole  at  5 0  intervals  by 
the  subroutine  STRCNF.  Cumulative  stress  concentration  due  to  combined  loading  may  be  estimated 
by  simple  addition  of  the  respective  stress  concentration  factors. 

The  next  output  (app.  B,  item  17)  is  under  the  heading  POISSON  RATIO  DIFFERENCES  and 
results  from  the  subroutine  NUDIFS.  For  each  ply,  the  Poisson’s  ratio  differences,  and 

(v\-  vc),  and  the  products  Klxx  {v\-  vc ),  Klyy  (yj-  vc),  and  Klxy  (y\~  vc)  are  printed  out  at  6  intervals 
of  5°  around  a  circular  hole.  This  is  suppressed  if  the  boolean  NONUDF  is  set  to  TRUE.  This  item 
shows  the  locations  of  probable  delamination  for  each  ply.  These  are  the  locations  where  products 
such  as  KXxx  ( v\—v'i ),  for  example,  are  maximum. 


Figure  4. — Definitions  of  regions  for  ply  microstress  calculations. 


The  next  item  in  appendix  B  shows  the  ply  stress  and  strain  influence  coefficient  arrays  and  ply 
stress  influence  coefficient  arrays  (app.  B,  item  18).  These  are  computed  in  the  subroutine  COMSA 
and  are  printed  by  the  main  program  ICAN.  The  first  table  gives  the  influence  coefficients  based  on 
unit  loads  or  moments/inch.  The  second  table  gives  the  influence  coefficients  in  terms  of  unit  applied 
stresses.  Explanations  of  usage  of  these  tables  are  provided  at  the  end  of  each  table. 

The  output  from  the  subroutine  MSCBFL  is  printed  out  next  under  the  heading  LAMINATE 
FAILURE  STRESS  ANALYSIS  (app.  B,  item  19).  The  analysis  is  based  on  first -ply  failure  criteria. 
Results  are  printed  in  a  tabular  form  for  each  ply,  and  a  summary  of  the  analysis  is  shown  in  the  end 
(app.  B,  item  20).  The  summary  shows  the  critical  ply,  the  failure  mode,  and  the  load  for  each  of  the 
applied  load  types,  ocxxT,  ocxxC ,  ocyyT,  ocyyC>  and  ocxyS„  respectively.  The  first  table  shows  results 
based  on  first-ply  failure,  and  the  second  table  shows  results  based  on  fiber  failure  by  breakage. 


A  Typical  IBM  Terminal  Session 

To  run  ICAN,  the  user  must  first  install  and  compile  the  program  on  his/her  computer  according 
to  the  system  to  be  used.  The  procedure  used  on  the  Lewis  Research  Center  IBM  370  is  described  in 
detail  here  starting  from  log  on.  The  computer  prompt  signals  are  identified  with  uppercase  letters. 
User  entries  are  in  lowercase  letters.  The  following  are  prerequisites  for  the  user  to  be  able  to  run 
ICAN: 

(1)  A  knowledge  of  how  to  compile  and  store  the  object  module  in  the  public  storage  space. 

(2)  A  knowledge  of  redit  or  tedit  processors  so  as  to  be  able  to  create  vs  datasets  of  the  input  data 
deck.  The  details  of  the  input  data  format  have  already  been  described  in  earlier  paragraphs. 

(3)  A  knowledge  of  commands  like  rmds,  mds,  ddef ,  libdef,  print,  and  erase.  These  are  a  few  of  the 
commands  commonly  used  in  running  a  program  on  the  IBM  370. 

The  user  is  advised  to  migrate  the  object  deck,  the  input  dataset,  and  the  material  property  data 
base  so  as  to  conserve  his/her  permanent  storage.  The  object  deck,  which  is  a  binary  version  of  the 
compiled  source  program,  is  referred  to  here  as  OBJ. ICAN.  The  data  base  of  material  properties  is 
referred  to  as  FBMTDATA.BANK. 

The  session  is  started  by  logging  on  at  the  terminal.  This  is  achieved  by  typing  logon,  userid,  and 
password.  The  system  replies 

TSS/370  RELEASE  3.0  PRPZ3  FTF18 

SOME  MESSAGE  TASKID  =  OBD7  POOLID  =  LRCFM  - 

LOGON  AT  11:30  ON  01/15/84 

The  user  is  now  ready  for  the  session.  The  first  phase  of  the  session  consists  of  restoring  the 
necessary  data  sets  to  temporary  storage.  This  is  achieved  by  the  following  commands: 

rmds  obj.  ican,  aaa 

SUCCESSFUL  (TEMP)  RESTORE  OBJ. ICAN  AS  (AAA) 
rmds  fbmtdata.bank,  ccc 

SUCCESSFUL  (TEMP)  RESTORE  FBMTDATA.BANK  AS  (CCC) 
rmds  ican. sample. input,  bbb 

SUCCESSFUL  (TEMP)  RESTORE  ICAN.SAMPLE.INPUT  AS  (BBB) 

At  this  point,  the  user  has  all  the  necessary  data  sets  to  run  ICAN  in  his/her  temporary  storage. 

The  input/output  fortran  units  that  are  utilized  by  ICAN  for  its  various  input/output  operations 
need  to  be  defined  next.  This  forms  the  second  phase  of  the  session  and  is  achieved  by 

ddef  ft05f001,  vs,  bbb 
ddef  ft06f001,  vs,  icanout.bbb,  ret=t 
ddef  ft08f001,  vs,  ccc,  ret  =  t 
ddef  ft07f001,  vs,  T7,  ret  =  t 

During  these  operations,  the  system  usually  responds  by  the  minimum  prompt,  the  underscore  (_). 

The  third  phase  consists  of  loading  and  executing  the  object  deck  and  printing  out  the  results.  This 
is  done  by  the  following  commands: 


libdef  Ids,  aaa 
load  gpcom$$$ 
ican 

TERMINATED:  STOP 
print  icanout.bbb,  prtsp  =  edit 
PRINT  BSN  =  8835,  1200  LINES 

The  last  phase  involves  cleaning  up  the  user’s  storage  place  and  is  achieved  by  issuing  the 
commands 

release  Ids 
release  ft 
erase  aaa 
erase  bbb 
erase  ccc 
erase  t7 

The  user  may  either  logoff  or  proceed  to  execute  another  run  for  a  different  set  of  input  data  after  the 
preceding  set  of  commands. 


Programmers  Manual 

A  brief  description  of  the  main  program  (or  control  program)  and  theoretical  equations 
programmed  in  the  code  are  presented  in  this  portion  of  the  report.  The  subroutine  descriptions 
follow  the  order  of  execution  as  shown  in  the  flowchart  (fig.  2(b))  rather  than  the  physical  sequential 
order  (fig.  2(a)).  It  is  assumed  in  the  following  discussion  that  the  user  has  a  working  knowledge  of 
computer  programming  and  that  he/she  is  familiar  with  the  terminology  appropriate  to  multilayered 
composite  mechanics. 

The  assumptions  and  details  leading  to  the  derivation  of  the  equations  programmed  in  the  code  are 
not  included  here.  However,  they  are  described  in  the  references  cited.  It  is  suggested  that  the 
interested  user  have  these  references  available  to  him/her. 

The  information  provided  in  this  portion  of  the  code  together  with  the  source  program  listing 
enables  the  user  to  modify,  implement,  and  extend  the  code  according  to  need. 


Main  Program 

The  main  program  contains  the  global  variables,  the  various  subroutines,  the  input  data  and 
format,  the  various  program  control  statements,  and  the  output.  These  are  discussed  subsequently. 
The  flowchart  of  the  program  is  shown  in  figure  5. 

The  global  variables  are  given  in  the  following  list: 
boolean  CSANB,  BIDE,  RINDV,  COMSAT,  NONUDF 

integers  Ne,  Np(,  Npc,  Np  Nt,  M,  Qj,  Qs, 

Qp,  Qr*  Qf 

real  ^cs*  Pf*  Pm*  ^f*  ^  ^  6/  ft  Pit  tt 

real  arrays  Kv(,  Kp,  6{c,  /,•(  1 , 1 000) , 

(75, 1000),  Pc(l,62) 

maximum  dimensions  Eci,  Ecp  Ecm,  Acx,  Ccx,  Dcx,  D^x, 

v4^.(3,3),  ap  am,  ae,  NcTfX,  MdT(X, 

NcM(X>  McM,X *  ecsz>  edw0>3)>  Lsc(l,6), 

MCX,NCX  (3,Nt),  Dv(  10,6),  AINF  (6,1000,8), 

(fy)p,S’  (fx)p,S*  (4 )p,S* 

( iy)p,S *  (41000) 
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Figure  5. — ICAN  program  flowchart. 


string  arrays  title  (80  characters)  read  in. 

Freight  spaces  per  field,  Npi  fields) 

Ce\  (six  spaces  per  field,  six  fields) 

Ce 2  (six  spaces  per  field,  six  fields) 

Pcp  (six  spaces  per  field,  Npc  fields) 
codes  C 

current  dimensions  Ng,  Np$  Npc,  Np  Nms 

real  arrays  Kv^  K ^ 

current  dimensions  Pc(\,Npc), 

AINF  (6,Afy,8),  \yfXtPjs(l>Ne), 

tx,y,P,S,(l>Ne) 

The  subroutines  are  as  follows: 

INVA  inverse  of  an  array 

GACD3  generates  composite  three-dimensional  elastic  and  thermal  properties  and 

the  two-dimensional  thermal  properties 
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BLOCK  DATA  DISP  (String)  and  RESF  (String) 

GPCFD2  generates  composite  two-dimensional  elastic  constants  and  constitutive 

equations 

COMSA  generates  the  ply  strain  and  stress  states  due  to  applied  loads,  checks  for 

ply  failure  and  interply  delamination,  and  generates  the  ply  stress  and 
strain  influence  coefficients 

INHYB  generates  ply  level  properties  with  the  aid  of  subroutines  FIBMT,  HTM, 

COMPP,  and  FLEXX 

BANKRD/IDGER  generates  constituent  properties  by  using  the  data  base 

FBMTDATA.BANK  and  arranges  them  in  a  proper  format  so  as  to  input 
to  INHYD 

FESTRE  computes  Poisson’s  ratio  mismatch  between  the  plies  and  the  composite 

EDGSTR  computes  interlaminar  free  edge  stresses 

MCRSTR/MINCOF  generates  the  microstresses  and  the  corresponding  influence  coefficients 
STRCNF  generates  the  stress  concentration  factors  around  a  circular  hole 

NUDIFS  generates  the  Poisson’s  ratio  differences  within  the  plies  and  the  probable 

locations  of  delamination  around  the  free  edge  of  a  circular  hole 

MSCBFL  performs  failure  load  analysis  based  on  first  ply  failure/maximum-stress 

criteria  and  prints  the  summary 

AMINF  minimum  value  of  an  array 

AMAXF  maximum  value  of  an  array 

FLRLD  determines  the  failure  load,  failure  mode,  and  the  ply  location 

These  subroutines  are  described  in  detail  in  the  next  section. 

INPUT  title,  Nb  Nk,  Nms ,  CSANB,  BIDE,  RINDV, 

COMSAT,  NONUDF,  tb  0b  7),  Mb  fiber 
name,  matrix  name,  kvb  kjb  ksc,  Ncx,  Mcx , 

DMCx->  Pu>  P( 

(For  substitution  and  definition,  see  appendix  A.) 

Subroutine  Description 

Subroutine  INVA  (N,A,C), — This  subroutine  computes  the  inverse  of  a  square  matrix  A  by  Gauss 
elimination  and  stores  it  in  array  C.  The  check 

|A|*0 

is  made  and,  if  satisfied,  the  program  continues;  otherwise,  the  message  SINGULAR  MATRIX  is 
displayed.  The  subroutine  inputs  are  N,  the  matrix  order,  and  the  matrix  A.  The  output  is 

A-i-C 

Subroutine  GACD3(C).—  This  subroutine  generates  the  three-dimensional  hygrothermoelastic 
properties  of  the  composite  about  its  structural  (x,y,z)  and  material  (1,2,3)  axes.  The  angle  6  is 
measured  from  x  of  the  structural  axes  system.  (See  fig.  6.)  In  figure  6,  replace  xy  etc.  by  1 1  etc.  and 
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Figure  6.— Ply  orientation  geometry.  Composite  structural  axes,  x,y,z ;  composite  material  axes,  1,2,3;  ply  material  axes 
(coincides  with  fiber  direction),  l/,2/,3/. 

measure  6  from  the  material  axes  to  obtain  properties  about  the  material  axes.  These  composite 
properties  are  generated  from  the  following  equations: 


i  r Ni 

7  L  (Zg. 

lc  L/=i 


^-zs)mTm[Rt]+  £  his;] 


1 

-  [Ec]  Yh  {Zg+\-Zs)[RS\T[E^{o 


lc  ;  =  1 


1  f 

(/3c)  =  7[EJ  D  (zg+i-z^[RgiT[Eg\{pg) 

lC  /  =  1 


The  arrays  (ac),  [/3C),  (a#},  and  {ft$}  in  the  preceding  equations  are  given  by 

fa!c]  =  [a.cxxa.cyy  OLCZZ  Oicyz  OLCZX  OiCXy]  ^ 

l  &  cl  -  [ft  cxx  ft  cyy  ftczz  ft  cyz  ft  czx  ft  cxy\  ^ 


M  =  [ani«m«f33  0  0  0]r 
lft&]  =  [ftnxftmftm  0  0  0]r 

For  all  practical  purposes,  the  two-dimensional  thermal  coefficients  of  expansion  about  the 
composite  structural  axes  are  the  same  as  acxx ,  acyy ,  and  acxy  in  the  array  {ac}  for  the  three- 
dimensional  case. 
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The  matrix  [Ec]  1  is  given  by 


Note  that  for  the  case  of  an  anisotropic  material,  the  elements  (1,6),  (2,6),  (3,6),  and  (4,5)  and  their 
symmetric  parts  will  not  be  zero. 

The  matrices  [Eq]  ~ 1  and  [7?#]  “ 1  are  given  by 
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cos2  6 

sin2  6 

0 

0 

0 

1  sin  28 

2 

sin2  6 

cos2  8 

0 

0 

0 

1  •  ™ 
-  -  sin  20 

0 

0 

1 

0 

0 

0 

0 

0 

0 

cos  6 

sin  8 

0 

0 

0 

0 

-sin  8 

cos  8 

0 

-sin  2d 

sin  28 

0 

0 

0 

cos  20 

where  6  =  6$  for  properties  about  the  composite  material  and  0  =  6$+6cs  for  properties  about  the 
composite  structural  axes.  (See  fig.  6.) 

The  matrix  [5^]  is  given  by 

r  A2  -A2  0  0  0  -AB~I 


—  A2  A2  0  0 


0  AB 


0 

0 

0 


-AB 


0 

0 

0 

AB 


0  0 
0  0 
0  0 
0  0 


0  0 

0  0 

0  0 

0  B2 


Here  A  =  sin  20 j -sin  20/ and  B  =  cos  20/ -cos  20/ _  j  where  />1  and  denotes  the  ply  index. 

The  angles  0/  and  0/_j  (fig.  6)  are  given  by 
0/  =  0#  +  0C5 
0/_l=0^’_l+0C5 

The  composite  heat  capacity  is  the  same  for  both  the  two-  and  the  three-dimensional  cases.  It  is  given  by 
1  Nf 

hc=  7  ^  ^ & 
lc  /=  1 

and  tc  is  given  by 


Ng 

^  C  tft 

/-I 
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The  composite  three-dimensional  heat  conductivities  along  the  composite  material  axes,  assuming  an 
orthotropic  composite,  are  given  by 


Kcn  =  -  S  k(Kn\  cos2  fy+^22  s'm2  fy)/ 
lc  i- 1 

1  ^ 

K&1-  7  S  %C^£11  sin2  ^+^22  cos2 

/-I 

Ng  /  . 

f-',1  2  Gr-) 

ac33  * c  i~  i  \Am/  i 

The  angle  0^is  measured  from  the  material  axes  (fig.  6) 

The  composite  two-dimensional  heat  conductivities  along  the  composite  structural  axes  are  given 
by  (see  ref.  9  for  the  transformation  equations) 

i  Nl 

Kcxx=  7  L  kWrn  cos2  9+K( 22  sin2  0)/ 

lc  1 

N, 

Kcyy  =  -  £  /*(*:,  11  sin2  0+*m  cos2  0)f 

*  C  /  =  1 

1 

^cyx  ”  Acxy  ~  7  ^  C^£22  1 )/  ^in  20/ 

/  =  1 

Aczz  =  A"c33 


The  angle  0  in  the  last  set  of  equations  is  measured  from  the  composite  structural  axes  and  is  equal  to 
0C5  +  0^.  The  inputs  to  the  subroutine  are  Afy,  1,  z#,  0C5,  0#,  [£/],  ///,  (a#),  %,  and  [i^],  which  are 
all  global.  The  variable  N \  is  input  data.  The  remaining  quantities  are  either  generated  or  are 
transferred  from  information  stored  in  PL(11,I),  PL(13,I),  PL(15,I-23,I),  PL(8,I),  PL(24,I)  to 
PL(26,I),  PL(30,I),  PL(27,I),  and  PL(29,I).  The  outputs  are  tc  and  the  arrays  are  [Ec]~~l,  [ac],  [Ec], 
hc,  and  [Kc).  The  composite  thickness  /cis  stored  in  PC(2).  The  arrays  [Ec]  ~ J,  [ac],  and  [Ec\  for  both 
composite  material  and  structural  axes  are  printed  out  under  the  headings  3-D  COMPOSITE 
STRAIN  STRESS  TEMPERATURE  MOISTURE  RELATIONS-STRUCTURAL  AXES  and  3-D 
COMPOSITE  STRESS  STRAIN  RELATIONS-STRUCTURAL  AXES. 

The  composite  material  axes  properties  [E^[  and  (ac)  are  stored  in  PC(3)  to  PC(14)  as  global 
variables.  The  corresponding  moduli  are  stored  in  PC(19)  to  PC(30).  The  three-dimensional  heat 
conductivities  and  heat  capacity  along  the  material  axes  are  stored  in  PC(15)  to  PC(18).  The  two- 
dimensional  thermal  coefficients  of  expansion  along  the  structural  axes  are  stored  in  PC(48)  to 
PC(50).  The  two-dimensional  heat  conductivities  and  heat  capacity  along  the  structural  axes  are 
stored  in  PC(51)  to  PC(54).  Note  that  the  heat  capacity  is  a  scalar  quantity  and  is  independent  of  the 
reference  axes.  Therefore,  PC(54)  equals  PC(18).  The  moisture  diffusivities  and  expansion 
coefficients  are  stored  in  entries  PC(55)  to  PC(62). 

Subroutine  BLOCK  DATA.— In  this  block,  the  strings  Cel  and  Ce2,  which  are  printed  out  with  the 
composite  constitutive  equations,  are  defined.  The  string  Ce\  contains  the  resultant  force  notation 
Ncx,  Ncy,  Ncxy,  Mcx,  Mcy ,  and  Mcxy .  The  string  Ce2  contains  the  notation  for  the  corresponding 
displacements. 

Subroutine  GPCFD2  (RESF,  DISP,  PROPC).— This  subroutine  generates  the  required  section 
properties  and  the  force-deformation  temperature-moisture  relations  for  a  two-dimensional 


multilayered  composite.  It  also  generates  the  plane-stress  elastic  constants  for  the  composite.  The 
force-deformation  temperature-moisture  relations  generated  in  this  procedure  are  defined  in  the 
following  equation: 


(Ay 1 

WcJ  tQxl 

ecsx 

,  Krfx) 

1 

1NcM,x ) 

{Mcx}\ 

[C^cJ  tf>cJ 

\  wcbx  ] 

[  WcTtxi  I 

{McM,xl 

The  generic  equations  for  the  elements  in  the  arrays  [Acx],  [CcJ,  [T>c J ,  {Mcf(X), 

and  \McM(X\  are 

Nt  N,- 1 

lAcxi=  £  (zs+1-ZeiMii]TlEBi-W+  £  Hj[Sj] 

/= 1  7=1 

N,  Nt- 1 

[Ccx]  =  i  £  (4+ 1  -4)[^Jr[^]  - 1 

z  (=1  7=1 

Nt  N[- 1 

[t>j = l  £  (4 + 1  - 4>uy  ■ 1  i*«] + 1  £ 

3  /-i  z  7=1 

Nt 

CT tx\  =  £  A  Tg  (zg  + 1  —  Zg)  \R(j\  *  [aff] 

(=1 

N, 

{NCMtx\ =  £  tee + 1 "  *#)  ~ 1 

(=i 

lMcT;Xi  =  9  £  ATe(zl  +  1  -  ~~  1  (“fff 

z  1=1 

1 

z  /= 1 

where  A  Tg=Tg-  Tcui 


The  arrays  {ag\,  [%},  [£#],  and  [Sy]  are 

[««)  =  [«  11  a22  0]f 

[/y  =  [£il  /322  0]f 


[%]  = 


COS2  0 

sin2  6 

~  sin  26 

sin2  0 

cos2  6 

-  ^  sin  2# 

-sin  26 

sin  26 

cos  26 
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Sj22  —  ^ y 1 1  —  ^  (sin  26  j  sin  2 6,_  j 
5y21=Syl2-  _Sy'll 

S/32  =  S7-23  =  \  (sin  26  j  -  sin  26,  _  i)(cos  26  j  -  cos  26,  _j) 

S/31  =5yi3=  -  S/23 

5y33=  i(cos  20 1 -cos  20/ _j)2 

Here  0/  equals  the  0C5  +  0^  (fig.  6).  The  reduced  bending  rigidities  (ref.  6)  are  generated  in  this 
procedure  according  to  the  equation 

&CX  “  \-^CX  ~~  ^cx^cx  Q:d 

The  reduced  axial  stiffnesses  are  generated  in  the  procedure  according  to  the  equation 

A  fx  =  CX  “  ^cx^cx  ^cd 

The  two-dimensional  composite  elastic  constants  are  generated  from  the  following  equation 
(assuming  7#  =  7}  for  /  —  1  to  Afy  and  for  /=  1  to  A/^): 

/  Nt  Nt  \ 

[£CJ-1=  i  (  E  (zft+i-z«rt£ff]-1[^]+  E  «,[«/]) 

tc  \  /  —  1  7-t  ' 

where 

N* 

tc=  E  t(i 
(=1 

The  inputs  to  this  subroutine  are  tg,  T$,  M&)  6j  (relative  to  composite  structural  axes),  Hj,  and  the 
ply  elastic  constants.  These  quantities  are  global  and  are  located,  respectively,  in  PL(7,I),  PL(50,I), 
PL(72,I),  PL(1 4,1),  PL(9,I),  and  PL(31,I)  to  PL(42,I).  The  arrays  [Re]T,  [£ff]_1,  [^&]>  and  [SJ\  and 
the  dimensions  Zy  are  generated  within  this  subroutine. 

The  outputs  are  the  force-deformation  temperature-moisture  relations,  which  are  stored  in  the 
global  arrays  A CX=Acx,  RAC  =A§X,  CPC  =  Ccx,  FL X=Dcx,  RDC  =D§X>  NSDT -Ncj(X, 
MSDT  =McT(X,  NSDH=7VcMfJC,  and  MSDH  =McM(X.  These  are  printed  out  under  the  heading 


24 


FORCES  FORCE  DISPLACEMENT  RELATIONS  DISPL  T-FORCES  H-FORCES.  The 
reduced  bending  rigidities  are  printed  out  under  the  heading  REDUCED  BENDING  RIGIDITIES. 
The  reduced  axial  stiffnesses  are  printed  out  under  the  heading  REDUCED  STIFFNESS  MATRIX. 
The  inverse  of  the  constitutive  equations 


[Acxi 

[Ccxi 

[Ccx] 

[Dcx\ 

are  printed  out  under  the  heading  DISP  DISPLACEMENT  FORCE  RELATIONS  FORCES. 
The  distances  zc,  ze,  and  za  are  stored  in  PC(31,I),  PL(10,I),  and  PL(ll.I),  respectively.  The  two- 
dimensional  composite  stress-strain  relations  are  stored  in  PC(33)  to  PC(38),  and  the  two- 
dimensional  composite  moduli  and  Poisson’s  ratios  are  stored  in  PC(39)  to  PC(47).  The  two- 
dimensional  thermal  properties  are  stored  in  PC(48)  to  PC(54),  as  is  described  in  the  section 
subroutine  GACD3. 

Subroutine  COMSA  (M).  — In  this  subroutine  the  stress  and  strain  states  of  each  ply  are  computed 
given  the  edge  membrane  forces,  the  ply  temperature,  and  the  changes  in  curvature.  In  addition,  two- 
ply,  combined  stress-strength  criteria  and  the  interply  delamination  criterion  are  generated.  Also 
generated  are  the  ply  stress-strain  influence  coefficients.  The  equations  programmed  for  the  z'th  strain 
and  stress  states  are 


M  =  m[Acx]  -  '((Ncc)  +  WcTgxl  +  {NcM,xl  +[CcJ[wcto))  ~z[R(i\{^cbx) 

(ffffj  =  [-£#] _  1  [Acx]  ~ 1  {{Ncx}  +  [NcT(X]  +  [NcM(X]  +  [Ccx]  [wcfof}) 

-  [Ee\  ~ 1  (wc&c)) 


The  reference  plane  strains  ecsx  and  the  curvature  changes  are  computed  from 


_ 

[ACx\  [CJ  ~ 

1 

{NCX\  1 

+ 

+ 

V^cMtx\ 

K  bx) 

[Ccx]  Wcx\ 

i 

Wcx] )  1 

WcT,x)  ] 

|  WcM,x J 

when  either  the  membrane  force  or  the  moments  or  both  are  given. 

The  strains  are  generated  locally  in  EPSL  and  SIGL,  respectively,  and  are  stored  in  PL(64,I)  to 
PL(69,I).  The  matrices  [Rft]  and  [£>,]  are  generated  locally  from  information  transferred  from 
PL(14,I)  and  PL(31,I)  to  PL(42,I).  The  distance  zB,  the  ply  temperature  Te,  and  the  ply  moisture  Mn 
are  transferred  from  PL(1 1,1),  PL(50,I),  and  PL(72,I),  respectively.  The  remaining  matrices  are 

Acx  -  ACX 
Ccx  -  CPC 
NcTlX  -  NSDT 
»cM,x  -  NSDH 
Ncx  -  NSBm 
McT(x  -  msdt 
Mcm(X  -  MSDH 
Mcx  MSBm 

and  wcbx- WXXm  (local  curvature  from  bending  analysis),  where  m  denotes  the  load  condition. 

It  is  important  to  note  that  the  stress  analysis  in  the  coded  form  also  handles  the  case  where  both 
the  reference  plane  membrane  strains  and  the  local  curvatures  are  given.  In  this  case  the  ply  strains 
are  given  by 


lecxi\  —  lecsx\  ziwcbx\ 


where  (eCA7-j  is  the  / th  ply  strain  along  the  structural  axis,  [ecsx\  is  the  reference  plane  membrane  strain, 
z  is  the  distance  from  the  reference  plane  to  the  centroid  of  the  /th  ply,  and  [wcbx]  is  the  local 
curvature.  These  variables  are  read  in  the  array  Dvnr  where  m  denotes  the  load  condition. 

The  corresponding  /th  ply  stresses  are  given  by 

(a,-|  =  [Eg]  - » {[Rg]  [tcxil  -  A  Tg\ag]  -  M g{(3j) 

Tfj-  Tcuj 

where  [oa\  is  the  /th  ply  stress  along  the  material  axes,  [Eg]  is  the  /th  ply  elastic  constant  about  the 
material  axes.  [Rs]  is  the  transformation  matrix  of  the  /th  ply,  (ecx,j  is  the  /th  ply  strain  along  the 
structural  axes  as  given  by  a  previous  equation,  Tg  is  the  temperature  of  the  /th  ply,  Tcui  is  the  cure 
temperature  of  the  ith  ply,  [ae]  is  the  thermal  coefficient  of  expansion  of  the  /th  ply  along  the 
material  axes,  Mg  is  the  moisture  content  of  the  /th  ply,  and  is  the  moisture  expansion  coefficient 
of  the  /th  ply  along  the  material  axes. 

The  displacement  force  relations  are  printed  out  in  the  following  format: 


DISPLACEMENT  DISPLACEMENT  FORCE  RELATIONS  FORCES 


lUcxl 


[^cxl  [Ccx] 
[Qx I  \Pcx^ 


Two  similar  sets  are  printed  out.  In  the  first  set,  the  displacement  and  force  vectors  are  in  symbolic 
form.  In  the  second  set,  the  displacement  and  force  vectors  have  their  numerical  values.  (See  outputs 
of  trial  cases,  app.  B.) 

The  failure  criterion  may  be  determined  by  either  of  the  following  methods: 

(1)  Modified  distortion  energy 


F=  l- 


v  \a  °m  , 

—  + 


3nia 


S(22 


PL(62,I) 


The  parameters  a  and  13  are  specified  as  follows: 


a  = 


0  = 


°n\ -0 
am  <° 


°(22—  0 
Of!22<0 


5ma- 


Sl22ee  ~ 


[SMT 

minCSfnc.  S(hcd) 

{S(22T  (3=T 

y.S(22C  ft  =  C 


a=r 
a  =  C 
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Kn2a$-K'n2a& 


(1  +4v(i2-V(i3)E(22  +  0-  "^23)^11 

\E(  1 1  £(>22(2  +  12  +  *'fl3)(2  +  ^21  +  ^23)] 1/2 


IBET(1,  7) 
BET(2,  7) 
BET(1,  8) 
BET (2,  8) 


a,  13=  T 
a  =  C,  &=T 
a=T,  0  =  C 
a,  13  =  C 


The  multiplier  of  K'n2ap  was  generated  in  the  main  program  and  is  stored  in  PL(61,I).  The  constant 
^2^13  constitute  theory-experiment  correlation  factors.  These  are  set  as  unity  in  COMSA.  However, 
the  user  can  modify  the  correlation  factors  if  he/she  wishes,  by  redefining  the  matrix  BET  in  the 
subroutine  COMSA. 

(2)  Hoffman’s  criterion  (ref.  9) 


5fiic-m'n(^mc>  %i cd) 


F=  1  - 


gfl  1  ~  (rUlaP22  +  gl22 
5fllc5nir  S122CS(22T 


smcsmr 


am  + 


S122C~S122T„ 

—q - o - °122  + 

^(22C;ie22T 


-^-1  -PL(71,I) 
Sms-*1 


F>0  no  failure 
F  =  0  incipient  failure 
F<0  failure 


The  interply  delamination  criterion  for  the  jth  interply  layer  at  the  mth  load  condition  is  governed  by 


[i-(^)l-pw63’,)  ”hen'>i 

A<PJ  ~  2  cyy  ~~  ~  ^0/ _  j)  +  ^  ^0/  ~~  20/  —  l)  tecjc]  =  l^cxl 

{i^cxl  “b  "b  {NcM(x}  “b  [QJ  [^cfor)} 


or  by  the  displacement  force  equation  described  previously. 

The  inputs  to  the  subroutine  are  the  ply  angle  measured  from  the  structural  axes  (0/,  from 
PL(14,I));  the  distance  from  the  reference  plane  to  the  centroid  of  the  ply  (z&,  from  PL(1 1,1));  the  ply 
temperature  (T&,  from  PL(50,I));  the  interply  delamination  limit  (A<pdely,  from  PL(60,I));  the  ply 
thermoelastic  properties  stored  in  PL(24  to  26,1)  and  PL(31  to  42,1);  the  ply  extensional  and  coupling 
rigidities,  Acx=ACX  and  Ccx=CPC;  the  local  curvatures  wcf?x=WXX;  the  adjustment  constants 
K%\2tt=  BET(1 ,  7),  Kg\2CT=  ®ET  (2,  7),  BET(1 ,  8),  and  Kp\2CC  =  BET(2,  8);  and  the  load 

conditions  AfCAr=NBS(m). 

The  subroutine  outputs  are  the  modified  distortion  energy  PL(62,I),  Hoffman’s  criterion  PL(71,I), 
the  interply  delamination  criterion  PL(63,I),  and  the  adjacent  ply  relative  rotation  (A<pj,  from 
PL(70,I)). 

Subroutine  EDGSTR. — This  subroutine  computes  the  interlaminar  stresses  azz,  ozy ,  and  o^near  a 
straight  free  edge  region  of  a  finite  width,  infinitely  long  plate  under  uniform  extension.  The 
equations  used  are  based  on  an  approximate  formulation  analogous  to  that  in  reference  18.  The 
calculations  are  performed  in  two  parts.  The  first  part  consists  of  computations  of  decay  lengths  for 


the  interlaminar  stresses.  The  decay  length  is  a  measure  of  a  free  edge  region  in  which  the 
interlaminar  stresses  may  be  significant.  This  is  achieved  in  the  main  program.  The  second  part  uses 
this  information  to  compute  the  interlaminar  stresses  in  the  subroutine  EDGSTR.  The  pertinent 
equations  are  discussed  in  the  following  paragraphs.  Note  that  in  the  case  of  hybrid  composite  plies, 
the  calculations  are  repeated  not  only  for  the  primary  composite  but  also  for  the  secondary  composite 
by  using  the  appropriate  ply  constituent  properties.  The  primary  and  the  secondary  composites  are 
distinguished  by  using  the  letters  P  and  S,  respectively,  in  the  Fortran  variables.  In  the  case  of  biaxial 
loading,  this  subroutine  is  bypassed  as  there  are  no  free  edges. 

Part  ]' — Decay  length  or  boundary  layer  width  computations .  The  interlaminar  stresses  near  the 
free  edge  are  assumed  to  decay  exponentially.  The  decay  length  is  calculated  with  the  aid  of  the 
following  equations: 

M-T('i) 


where 

a  =  f/?1  (0.001) 

and 


M  = 


—  \\l  -  —  - 1~] 

E(yy  L  '  4(1  -  k^kf  J 


) 


1/2 


The  calculations  are  repeated  for  each  layer.  Quantities  tb  and  X,  are  stored  in  arrays  YPL  and 
PLMDAY.  These  quantities  pertain  to  the  free  edge  parallel  to  the  load  axis  X.  The  corresponding 
quantities  for  the  load  axis  parallel  to  Y  are  stored  in  arrays  XPL  and  PLMDAX.  These  are 
computed  by  replacing  E^y  with  E in  the  preceding  equations.  For  the  intraply  hybrid  composite, 
the  respective  arrays  for  the  secondary  composite  are  denoted  by  YSL,  SLMDAY,  XSL,  and 
SLMDAX.  Note  that  the  letter  P  is  replaced  by  S.  This  notation  is  followed  consistently  throughout 
the  text.  The  labeled  common  block  ILAB6  is  used  to  store  and  pass  these  data  to  subroutine 
EDGSTR. 

Part  2.— Interlaminar  stress  computations.  In  the  EDGSTR  subroutine,  the  ply  stresses  PL(67,I) 
to  PL(69,I)  are  transformed  to  the  structural  coordinate  system  x,  y,  and  z.  These  stresses  are  stored 
in  the  matrix  SIGMA  (3,1)  for  each  layer.  The  interlaminar  stresses  (afe)  are  computed  with  the  aid  of 
the  following  relations: 


for  i  =  N(_  jto  N/2  + 1 


2  N( 

<J(yy 


The  interlaminar  shear  stresses  {a^y}  and  are  calculated  by 


7+1 

oiyy^i 

i  oc  j=N( 

atzy  ?a  _  j) 


fori  =  Ntto  ^  + 1 
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and 


for  i=Ngto  +  1 


In  these  equations,  the  computations  are  started  from  the  top  layer  =  After  the  midplane  is 
approached  (i-N/ 2+1),  the  calculations  are  repeated  starting  from  the  bottom  layer  (/  =  1)  and 
continued  until  i  becomes  (TV/-  1). 

The  interlaminar  stresses  are  stored  in  the  arrays  YSZZP,  SZYP,  and  SZXP  for  the  primary 
composite  and  in  the  arrays  YSZZS,  SZYS,  and  SZXS  for  the  secondary  composite.  They  are, 
however,  made  dimensionless  by  dividing  by  the  applied  normal  stress  o^x. 

Subroutine  STRCNF. — This  subroutine  calculates  the  stress  concentration  factors  around  a 
circular  hole  due  to  membrane  loading.  The  equations  used  are  taken  from  reference  19  and  are 
strictly  applicable  for  infinite  plates.  Three  factors  are  computed  in  the  subroutine  and  are  defined  by 
the  following  equations: 
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Klxx~ 


°99 

°XXOO 


K ,  =  aee 
Al  yy  a — 

uyyoo 


K\xy~ 


Gee 

Gxyoo 


Quantities  <JXXOO,  oyyOQy  and  oxyO0  are  the  applied  stresses,  and  oqq  is  the  hoop  stress  at  any  angle  0  from 
the  load  axis.  The  stress  concentration  factors  are  stored  in  the  local  arrays  XK1,  XK3,  and  TEMP. 
The  expressions  for  KXxx ,  KXyy,  and  KXxy  are  the  following: 


In  the  preceding  expressions,  Ectt  and  Ecn  are  the  composite  moduli  in  the  tangential  and  radial 
directions  at  angle  0.  Angle  6  is  measured  from  the  x-axis  for  KXxx  and  KXxy  and  from  thejy-axis  for 
KXyy.  The  program  rearranges  the  computed  KXyy  values  so  that  they  correspond  to  the  same  location 
as  those  of  KXxx  and  KXxy. 

Subroutine  NUDIFS. — In  this  subroutine,  the  Poisson’s  ratio  differences  between  the  adjacent 
plies  and  the  composite  are  computed  around  a  circular  hole  at  5°  intervals.  The  products  of  the 
differences  and  the  corresponding  stress  concentration  factors  are  computed  next.  These  products  are 
expected  to  provide  insight  into  the  probable  delamination  locations.  It  is  assumed  that  onset  of 
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delamination  is  likely  to  occur  at  the  locations  for  which  the  product  of  Poisson’s  ratio  mismatch 
with  the  corresponding  stress  concentration  factor  is  a  maximum.  Accordingly,  these  products  are 
computed  at  5°  intervals  and  the  maxima  are  calculated.  Two  sets  of  tables  are  the  output  from  this 
subroutine.  The  first  table  comes  out  optionally  if  the  boolean  NONUDF  is  set  to  FALSE.  It  contains 
all  the  details  of  the  computations.  The  second  table  consists  of  the  summary  of  results,  with  notes  on 
the  maxima  and  the  locations.  The  following  are  the  programmed  equations. 

At  any  angle  d  the  Poisson’s  ratio  is  computed  by 


vcrt  ~  ^crr 


vcxy  _  f  I 
-  ECXX  k 


+  2vt 


cxy 


+  - 


*cyy 


&cxy  ' 


I  cos2  0  sin2  0 


The  ply  Poisson’s  ratio  is  given  by 


(%/] =E\ 


trr 


M  -  ( 1  +  -1 - L  )cos2  0sin2 

_E(U  '  eQ2  Gftl' 


The  difference  in  Poisson’s  ratio  between  the  ith  and  (i  +  l)th  plies  is  given  by  (v$ 1  -  v'rt),  and  the 
difference  with  respect  to  the  composite  is  given  by  -  vcrt).  These  are  stored  in  the  arrays  A2  and 
A3,  respectively.  The  products  of  Kixx>  K\yy,  and  K\xy  with  A3  are  computed  next  and  are  stored  in 
the  arrays  A5,  A6,  and  A7,  respectively.  The  maxima  and  their  location  in  each  of  the  four  quadrants 
(0-90,  90-180,’  180-270,  and  270-0)  are  computed  by  calling  the  subroutine  AMAXF  for  the  three 
arrays  A5,  A6,  and  A7.  The  values  of  stress  concentration  factors  are  passed  through  the  labeled 
common  block  ILAB8  from  the  subroutine  STRCNF. 

Subroutine  MSCBFL  (A INF).— A  complete  laminate  failure  stress  analysis,  based  on  first-ply 
failure  and  the  maximum  strength  criteria,  is  performed  in  this  subroutine.  The  inputs  to  this  routine 
are  the  ply  allowables  Sfnc,  ShiT’  $(120  S(22T’  an<^  ^fl2S  an<^  the  Pjy  influence  coefficient  matrix 
AINF.  The  ply  stress  allowables  are  generated  by  the  INHYD  routines  and  are  stored  in  the  ply 
properties  array  PL.  These  are  accessed  through  the  labeled  common  block  ILAB2.  The  ply  stress 
influence  coefficients  are  generated  by  COMSA  and  the  main  program  and  are  passed  to  the  present 
routine  by  the  subroutine  argument. 

The  failure  stress  for  a  particular  ply  due  to  a  specific  loading  is  given  by  the  ratio  of  the  allowable 
strength  to  the  ply  stress  influence  coefficient.  For  example,  the  failure  stress  due  to  a  tensile  load  is 
given  by 


«=  $11 T 
c  Factl' 

where  Factl'  is  the  stress  influence  coefficient  for  rth  ply  due  to  unit  tensile  loading,  S\UT  [s  the 
strength  allowable  for  ith  ply  in  longitudinal  tension,  and  S', is  the  failure  stress  for  the  ith  ply  due  to  a 
tensile  loading.  The  failure  stresses  are  stored  in  the  matrix  FAILD.  In  the  case  of 
temperature/moisture  presence,  the  allowable  strengths  are  updated  to  take  into  account  temperature 
or  moisture  stresses;  the  failure  stresses  are  computed  with  and  without  the  effects  of  temperature- 
and  moisture-induced  stresses  for  comparison.  The  program  considers  primarily  five  different 
loadings,  longitudinal  compression  and  tension,  transverse  compression,  and  tension  and  inplane 
shear. 

After  the  failure  load  computations  for  each  ply  are  determined,  the  active  failure  mode  and  the 
corresponding  failure  strength  for  each  type  of  loading  are  determined  by  calling  the  subroutine 
AMINF.  This  subroutine  returns  the  value  of  the  minimum  failure  load,  the  ply  number,  and  the 
failure  mode  as  output.  The  output  from  this  subroutine  is  printed  under  the  heading  LAMINATE 
FAILURE  STRESS  ANALYSIS. 

Subroutine  MCRSTR— This  subroutine  generates  the  microstresses  in  the  ply  constituents  due  to 
the  inplane  loading.  These  are  stored  in  the  ply  microproperty  arrays  PLMP  and  PLMS  for  the 


primary  and  the  secondary  composites.  The  ply  constituent  properties  and  the  applied  loads  are 
inputs  to  this  subroutine.  They  are  accessed  with  the  aid  of  the  common  blocks  PBANK,  MFBANK, 
ILAB2,  ILAB5,  and  ILAB9.  The  PLMP  and  PLMS  each  contain  41  entries  which  are  explained  in 
the  following  list: 


Code 

name 

Algebraic  notation 

Fortran  variable 

PLM(1,I) 

'  °m\\L 

SMIL 

PLM(2,I) 

amllT 

SM1T 

PLM(3,I) 

ofin 

SF1L 

PLM(4,I) 

OfllT 

SF1T 

PLM(5,I) 

(A) 

am22L 

SM2AL 

PLM(6,I) 

(A) 

am22T 

SM2AT 

PLM(7,I) 

(B) 

am22L 

SM2BL 

PLM(8,I) 

(B) 

am22T 

SM2BT 

PLM(9,I) 

JB) 

V/22L 

SF2BL 

PLM(10,I) 

JB) 

af22T 

SF2BT 

PLM(11,I) 

L 

SM3AL 

PLM(12,I) 

am33T 

SM3AT 

PLM(13,I) 

ffm33L 

SM3BL 

PLM(14,I) 

am33T 

SM3BT 

PLM(15,I) 

(B) 

a/33L 

SF3BL 

PLM(16,I) 

(B) 

°/33T 

SF3BT 

PLM(17,I) 

(A) 

aml2 

SM12A 

PLM(18,I) 

(B) 

aml2 

SM12B 

PLM(19,I) 

(B) 

12 

SF12B 

PLM(20,I) 

(A) 

cml3 

SM13A 

PLM(21,I) 

(B) 

aml3 

SM13B 

PLM(22,I) 

(B) 
af  13 

SF13B 

PLM(23,I) 

(A) 

<Jm23 

SM23A 

PLM(24,I) 

b 

SM23B 

PLM(25,I) 

JB) 

a/23 

SF23B 

PLM(26,I) 

°mll  ^ 

SMI  IDT 

PLM(27,I) 

°/ll 

SF11DT 

PLM(28,I) 

(A) 

am22 

SM2ADT 

PLM(29,I) 

JB) 

am22 

JB)  / 
af22  l 

Microstresses  due 
•  to  temperature 

SM2BDT 

PLM(30,I) 

gradient  ATI 

SF2BDT 

PLM(31,I) 

(A)  \ 

°m33  1 

SM3ADT 

PLM(32,I) 

JB) 

ffm33 

SM3BDT 

PLM(33,I) 

JB) 

ff/33  j 

SF3BDT 

PLM(34,I) 

) 

°m\\ 

SM11DM 

PLM(35,I) 

a/n 

SF11DM 

PLM(36,I) 

(A)  j 
am22  J 

SM2ADM 

PLM(37,I) 

JB)  f 
°m22  y 

Microstresses  due 

SM2BDM 

PLM(38,I) 

JB)  / 
°J22 

to  moisture  Mj> 

SF2BDM 

PLM(39,I) 

JA) 

SM3ADM 

PLM(40,I) 

JB) 

am33 

SM3BDM 

PLM(41,I) 

JB) 

ff/33  j 

SF3BDM 

In  this  list,  entries  26  to  41  are  suppressed  automatically  if  the  temperature  gradients  and  moisture 
contents  are  not  present.  The  superscripts  A  and  B  refer  to  two  regions  as  described  in  figure  4. 

The  microstresses  are  calculated  with  the  aid  of  the  following  equations:  (For  notation  and  sign 
conventions,  see  figs.  4  and  6.) 

Ply  microstresses  due  to  a  longitudinal  stress  an  ]  are  given  by 

am  1 1  =  (Pm/Ef  1 1  )(T(\  i 


°/ii  =(Efil/En[)on] 


(A)  _ 
am22~ 


( vm 


*'m)(Em/Ei  j)<jn  i 


(B)  JB) 
am22  aJ22 


1  -  VkT 


(A)  _  (A) 
am33  °m22 
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Ply  microstresses  due  to  a  transverse  stress  oq 2  are  given  by 


am 


_  ( vn2^m\ 
am\\~  \^m  Em  / 

af\  1  =  {^f\2~vl\2^^ja!22 


ffm22  -  (Em/ ^2)fft22 


CTm22  ~  (Em/ E 2)<J  (22 


a/22  =  (Em/E2)°a2 


where  E2  is  given  by 


^}Em 

E2  =  (\-Jkf)Em  + 

1-V^ 

°m33  =  (vm/ vdi)(Em/E (22)0(22 


(B)  - 
am33~ 


1  -VE 


(5)  (5) 

°/33  “aw33 


Ply  microstresses  due  to  inplane  shear  stress  (a^)  are  given  by 

ffml2  =  (Gw/G12)fffl2 


aml2  =  (Gfl  2/G\2)°tt2 


°/12  =  (Gfl2/G12)(Tfl2 


where  G12  is  given  by 


Ply  microstresses  due  to  through-the-thickness  shear  stress  om  are  given  by 
(A) 

<23  =  (C  m  /G  (22)0122 

am23  =  (G23/ Gm)aa3 
where  G2 3  is  given  by 


Ply  microstresses  due  to  temperature  gradient  AT)  are  given  by 
o fa  u  (njfl  1  ^ /w)^ 
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(B)  -iB) 
°/33  ~°fl2 


AT(-  T(—  Tcu 

Ply  microstresses  due  to  moisture  M(  are  given  by 


am\\  ~  (Pa l  “ &m)M(Em 


Jfu-PtnMeEfn 


°m22  ~  (.&0.2 


JB)  _  l-V^.(A) 

(Sm22  r—  °m22 

'Fk f 


(A)  (A) 

am33  ~  °m22 


Jfi)  „(B) 

am33  am22 


(A)  (A) 

°/33  =  «^/33 


Subroutine  MINCOF. — This  subroutine  generates  the  microstress  influence  coefficients  for  each 
different  material  system  used  in  the  layup.  The  equations  used  are  similar  to  those  programmed  for 
MCRSTR.  However,  the  influence  coefficients  are  based  on  the  application  of  unit  load  in  a  specific 
direction  or  unit  temperature  difference  or  unit  moisture  content.  The  influence  coefficients  are 
stored  in  the  matrix  PINF.  This  matrix  has  17  entries.  They  are  described  in  the  following  list: 


Entry 

Algebraic  notation 

Fortran  variable 

PINF(1,K,NLD) 

amll 

SMI  1 

PINF(2,K,NLD) 

(A) 

am22 

SM22A 

PINF(3,K,NLD) 

m 

am22 

SM22B 

PINF(4,K,NLD) 

JA) 

cm\2 

SM12A 

PINF(5,K,NLD) 

(B) 

°m\2 

SM12B 

PINF(6,K,NLD) 

JA) 

CTml3 

SM13A 

PINF(7,K,NLD) 

(B) 

am\2 

SM13B 

PINF(8,K,NLD) 

JA) 

°m2'i 

SM23A 

PINF(9,K,NLD) 

JB) 

ffm23 

SM23B 
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PINF(10,K,NLD) 

(A) 

am33 

SM33A 

PINF(11,K,NLD) 

<B) 

°m  33 

SM33B 

PINF(12,K,NLD) 

afU 

SF1 1 

PINF(13,K,NLD) 

JB) 

af22 

SF22B 

PINF(14,K,NLD) 

JB) 

°/33 

SF33B 

PINF(15,K,NLD) 

°f\2 

SF12 

PINF(16,K,NLD) 

°/13 

SF13 

PINF(17,K,NLD) 

JB) 

<7/23 

SF23B 

The  dimension  K  varies  from  1  to  NMS,  where  NMS  is  the  number  of  material  systems.  NLD 
varies  from  1  to  7.  The  expression  NLD  =  1  to  5  refers  to  unit  applied  stresses  in  11,  22,  12,  13,  and 
23,  respectively.  The  expression  NLD  =  6  corresponds  to  unit  temperature  loading,  and  the 
expression  NLD  =  7  corresponds  to  unit  moisture  loading. 

The  microstress  influence  coefficients  are  computed  for  secondary  composites  and  optionally 
computed  for  intraply  hybrid  composites.  These  are  stored  in  the  matrix  SINF. 

Subroutines  AMAXF,  AMINF,  LOGO,  and  LOGOI.—  These  subroutines  perform  several 
auxiliary  duties.  AMAXF  finds  the  maximum  value  of  a  one-dimensional  array  and  its  location. 
AMINF  finds  the  minimum  value  of  a  one-dimensional  array  and  its  location.  These  two  subroutines 
are  utilized  by  MSCBFL  and  NUDIFS  in  conjunction  with  searching  for  failure  loads  and  the 
probable  locations  of  delamination.  LOGO  is  a  subroutine  to  generate  the  ICAN  emblem  for  the 
output.  The  description  of  the  material  and  the  structural  coordinate  system  by  appropriate  figures  is 
generated  by  the  subroutine  L0G02. 

Subroutine  INHYD  — This  subroutine  generates  the  composite  ply  properties,  necessary  for  the 
laminate  response  analysis.  The  inputs  to  this  routine  are  the  constituent  properties  which  are 
supplied  in  the  appropriate  format  by  the  subroutines  IDGED  and  BANKRD.  INHYD  calls  the 
subroutines  FIBMT,  COMPP,  and  HTM  to  perform  the  micromechanics  analysis,  including  the 
analysis  of  hygrothermal  effects.  The  ply  properties  are  stored  in  the  array  PROPS,  which  is  accessed 
by  the  main  program  through  the  labeled  common  block  PBANK.  INHYD  is  called  once  for  each 
different  material  system  by  the  main  program.  The  outputs  of  INHYD  show  the  properties  of  the 

fiber,  matrix,  and  composite.  . 

The  fiber  and  matrix  properties  for  the  primary  composite  are  read  in  from  the  input  provided  by 
IDGER.  These  are  stored  in  arrays  PF  and  PM.  Similarly,  the  arrays  SF  and  SM  are  used  to  store  the 
properties  of  secondary  composite  constituents  if  the  composite  is  of  the  hybrid  type.  The  program 
then  checks  for  temperature  and  moisture.  The  properties  of  the  matrix  are  updated  for  the  presence 
of  temperature  and  moisture.  The  following  are  the  equations  programmed  to  account  for  the 
hygrothermal  property  degradation: 

The  wet  glass  transition  temperature  is  computed  from 


Tgwr= (0.005M}-  AMe+  \)Tgdr 

where  Tgwr  is  the  wet  glass  transition  temperature,  Tgdr  is  the  dry  glass  transition  temperature  for  the 
resin  matrix,  and  M( is  the  percentage  of  moisture  by  weight. 

The  reduction  factors  Xmp  and  Xtp  are  computed  from 

Xmp  =  ^(Tgwr-Tu)/(Tgdr-T0) 

X tp  =  i  !Xmp 

where  T0  is  the  reference  temperature  (70  °F),  and  Tu  is  the  use  temperature. 


The  moduli  and  strengths  of  the  matrix  are  multiplied  by  Xmp  to  obtain  the  new  properties  for  the 
matrix.  The  density  is  given  by 

Pmw~  Pm  "b  3pmkmM(/ 1 00 

The  thermal  properties,  such  as  heat  capacity,  thermal  expansion  coefficient,  and  thermal 
conductivity  are  multiplied  by  the  second  factor  X,p  to  account  for  the  hygrothermal  conditioning. 

After  the  property  arrays  PF,  PM,  SF,  and  SM  are  properly  filled,  the  program  chooses  either 
FIBMT  or  HTM  subroutines  to  perform  micromechanics.  The  subroutine  HTM  is  chosen  if 
temperature/moisture  effects  are  to  be  taken  into  consideration.  Otherwise,  FIBMT  is  chosen  for  dry 
room  temperature  property  computations.  The  outputs  from  these  routines  are  primary  and 
secondary  composite  ply  properties.  They  are  stored  in  the  arrays  P  and  S,  respectively.  These 
properties  are  made  common  to  subroutine  COMPP  through  the  common  blocks  ILAB1  and 
ILAB3.  The  subroutine  COMPP  is  called  by  INHYD  for  hybrid  composites  to  compute  the  hybrid 
composite  ply  properties.  These  properties  are  stored  in  the  array  H.  One  of  the  arrays  P,  S,  or  H  are 
passed  to  ICAN  via  common  block  PBANK  and  the  array  PROPS.  For  example,  if  the  ply  is  made  of 
100  percent  primary  composite  only,  the  array  PROPS  is  assigned  to  have  the  same  entries  as  P,  etc. 

The  subroutine  INHYD  also  calls  FLEXX,  which  performs  a  flexural  strength  analysis.  However, 
these  are  only  for  additional  information  and  are  not  used  by  ICAN  at  the  present  time. 

Subroutine  FIBMT  (C,  F,  M,  VF,  VM,  VP,  KV,  IFLAG). — This  subroutine  generates  properties 
of  a  ply  by  using  the  constituent  properties  which  are  supplied  from  the  subroutine  INHYD.  The 
constituent  properties  are  stored  in  the  arrays  F  and  M;  F  contains  the  fiber  properties,  and  M 
contains  the  matrix  properties.  The  composite  properties  are  stored  in  the  array  C,  which  is  returned 
to  INHYD.  The  theory  behind  the  programmed  equations  is  discussed  in  reference  13.  The  following 
is  a  description  of  each  entry  in  the  arrays  C,  F,  and  M,  with  the  corresponding  algebraic  notation: 


Composite  Properties  Array  C(I) 


Entry 

Description 

Notation 

C(l) 

elastic  moduli 

An 

C(2) 

elastic  moduli 

C(3) 

elastic  moduli 

Em 

C(4) 

shear  moduli 

Gin 

C(5) 

shear  moduli 

G123 

C(6) 

shear  moduli 

Al  3 

C(7) 

Poisson’s  ratio 

vl\2 

C(8) 

Poisson’s  Ratio 

vm 

C(9) 

Poisson’s  Ratio 

vl\2 

C(10) 

thermal  expansion  coefficient 

«m 

C(11) 

thermal  expansion  coefficient 

af22 

C(12) 

thermal  expansion  coefficient 

aB3 

C(13) 

density 

Pe 

C(14) 

heat  capacity 

Q 

C(1 5) 

heat  conductivity 

K(n 

C(16) 

heat  conductivity 

Kn  2 

C(17) 

heat  conductivity 

Ki 33 

C(18) 

strength 

Sflir 

C(19) 

strength 

All  c 

C(20) 

strength 

S(22T 

C(21) 

strength 

sa.2C 

C(22) 

strength 

Se\2 

C(23) 

moisture  diffusivity 

All 

C(24) 

moisture  diffusivity 

% 

C(25) 

moisture  diffusivity 

A 33 

C(26) 

moisture  expansion  coefficient 

All 

C(27) 

moisture  expansion  coefficient 

13(22 

C(28) 

moisture  expansion  coefficient 

C(29) 

flexural  moduli 

Em 

C(30) 

flexural  moduli 

E(22 

C(31) 

strengths  (flexural) 

0(23 

C(32) 

strengths  (flexural) 

0(1  IF 

C(33) 

strengths  (flexural) 

0(22F 

C(34) 

strengths  (flexural) 

0(12 

C(35) 

ply  thickness 

‘( 

C(36) 

interply  thickness 

6S 

C(37) 

interfiber  spacing 

Fiber  Properties  Array 

Entry 

Description 

Notation 

F(l) 

elastic  moduli 

E/n 

F(2) 

elastic  moduli 

EJ22 

F(3) 

shear  moduli 

Gfl2 

F(4) 

shear  moduli 

Gf22 

F(5) 

Poisson’s  ratio 

vfl2 

F(6) 

Poisson’s  ratio 

vf23 

F(7) 

thermal  expansion  coefficient 

«/ii 

F(8) 

thermal  expansion  coefficient 

a/22 

F(9) 

density 

Pf 

F(10) 

number  of  fibers  per  end 

Nf 

F(ll) 

fiber  diameter 

df 

F(12) 

heat  capacity 

Cf 

F(13) 

heat  conductivity 

KfU 

F(14) 

heat  conductivity 

Kfn 

F(15) 

heat  conductivity 

Kfi  3 

F(16) 

strength 

SfT 

F(17) 

strength 

Matrix  Properties  Array 

SfC 

Entry 

Description 

Notation 

M(l) 

elastic  modulus 

F 

M(2) 

shear  modulus 

Gm 

M(3) 

Poisson’s  ratio 

vm 

M(4) 

thermal  expansion  coefficient 

Oim 

M(5) 

density 

Pm 

M(6) 

heat  capacity 

c 

M(7) 

heat  conductivity 

Km 

M(8) 

strength 

S mT 

M(9) 

strength 

S mC 

M(10) 

strength 

SmS 

M(1 1) 

moisture  coefficient 

@m 

M(12) 

diffusivity 

Dm 

The  following  are  the  programmed  equations  for  the  entries  in  array  C 


Normal  moduli: 


Efi  i  —  kfEfl  i  +  kmEm 

£  _  _ _ 

G2  \-JFf(\-Em/Em) 

Em=E0.2 

Shear  moduli: 


Poisson’s  Ratio: 

Vn2  =  '>m  +  kf{vjX2-vm) 


vn3~  vn2 


Coefficients  of  thermal  expansion: 

afi  l  k nXioifnE m/ Ej\  i)  —  otf\  j] 
O!01  = - 

’♦‘-(ffr1) 


<*122  =  “mi1  -'Tty 

“33  =  “{22 

Density: 

P?~  Pfkf~^~  Pmkm 

Heat  capacity: 

q  _  (kfCfPf+  kmCmPrn) 

1  pt 

Heat  conductivities: 

E$\  i  ~  kfKf\  1  T 

_  Vkf 

Ka2  =  {\-4kf)Km  + - Km 

Km  =  Km 

In  the  preceding  equations,  Km  should  be  replaced  by 
Ktn'J~k'v 

Km-(l-^TjKm+ - 

if  there  are  voids.  The  quantity  Kv  is  the  void  conductivity. 

Strengths: 

^flir  Sfr(kf+  kmEm/  Ef\\) 


1  +kjvmEfn 


En  i  [  E, 


'Jfn~tKm{Em~EfUSj 


+amkf 


The  longitudinal  compressive  strength  is  computed  based  on  three  different  criteria,  rule  of  mixtures, 
fiber  microbuckling,  and  delamination.  The  minimum  of  the  three  estimates  is  returned  as  Sfuc.  The 
equations  for  the  three  cases  are 


Snit'(rule  of  mixtures)  =  Sfc(kf+  kmEm/ £/i i) 
Sue  (delamination)  =  (13Sfl2  +  5mc) 


S(1  ic  (fiber  microbuckling)  = 


M'-k) 


The  transverse  strengths  are  calculated  from 


SQ2T=  Sot7-(FACT/DENOM) 


Sq2C  —  Snlc/Dl:  N  O  M 


„  _  i  1  ~b  G m/ G f\2)F2G (\2$ ms\  p a rr 

7  ~  ^  t?  '  J 


where  Fx  and  F2  are  defined  by  the  equations: 


The  variable  DENOM  is  a  Fortran  variable  given  by 


DENOM  = 


-[>-VF/( 


F 

1 


\+<p(<p-\)  +  -  (ip-1)2 


where  <p  is  given  by 


The  Fortran  variable  FACT  takes  the  value  km  if  IFLAG  is  unity.  Otherwise  FACT  takes  the  value 
unity.  This  variable  is  introduced  to  correlate  the  strengths  of  HMS  and  Kevlar  fiber  composites  with 
the  experimentally  observed  values.  The  main  program  INHYD  checks  for  these  fibers  and  assigns 
the  appropriate  values  for  IFLAG.  IFLAG  is  set  at  zero  for  other  fibers. 
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Moisture  diffusivities: 


1 

Dm  =  V-'fk})Dm 

DB2~Da2 

Moisture  expansion  coefficients: 

/5fl  i =  fimkmEm/En  j 

=  ~  V^/)(l  +kfmEfn/. En 2) 

&B3  =  0(22 

Flexural  moduli  ( EeuF ,  Enopj: 
EeuF~E(n 
EaiF~Ea2 

Flexural  Strengths: 

(fri  +  ^G^ 

SG3jp= - 

Gn/l 


Sfl2F=  !-55fl2 


Ply  thickness:  A  default  value  of  0.005  is  set  for  tg.  This  is  overridden  by  the  user  specified  value  in 
the  1CAN  main  program. 

Interply  thickness  and  interfiber  spacing: 


8f=8g 
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In  the  preceding  equations,  FACT  is  a  Fortran  variable  which  is  given  by 
FACT  =  5/6/ 

for  Kevlar  and  HMS  fibers.  For  all  other  fibers  FACT  =  1 . 

Subroutine  FLEXX  (C).— The  entries  C(32)  and  C(33)  of  the  ply  property  array  C  are  generated  in 
this  subroutine.  They  are,  respectively,  the  longitudinal  flexural  strength  and  the  transverse  flexural 
strength.  The  longitudinal  flexural  strength  is  given  by 
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2.5Sflir 


The  transverse  flexural  strength  is  given  by 
2-5^G2r 

Simp  = - 


Subroutine  COMPP  (IPFLAG,  ISFLAG). — This  subroutine  is  called  by  INHYD  to  generate  the 
properties  of  a  hybrid  ply.  The  equations  are  based  on  the  theory  proposed  in  reference  13.  The 
properties  are  stored  in  the  array  H.  The  entries  are,  however,  the  same  as  those  of  array  C  given  in 
the  description  for  subroutine  FIBMT.  The  inputs  to  this  routine  are  the  primary  composite 
properties  array  P,  the  secondary  composites  property  array  S,  and  the  percentage  of  the  secondary 
composite  ksc.  The  equations  are  the  following: 


Elastic  normal  moduli: 

En  m=Em  (P)  +  [En ,  (S)  -  Ea  j  (P)]  ksc 


Em(P)  _ 

£f22(H)  =  l+ksclEe22(P)/Ea2(S)-l] 


EB3(H)=Ee33(P)  +  [EB3(S)-EB3(P)]ksc 
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Poisson’s  ratios: 


Heat  capacity: 

Q(H)  =  {(1  -  ksc)[C/P)kj{P)p/P)  +  Cm(P)km(P)pJP)} 
+  ksc\Cj{S)kj(S)p/S)  +  Cm(S)km(S)pm(S)} 

+  [kj(P )k ^(P )  +  ksckp(S)]MpmS(Cmsi^  / Pf(H) 


where  pmst  and  Cmst  are  the  moisture  density  and  heat  capacity,  respectively. 


Heat  conductivities: 


ATfll(H)  =  (l-  ksc){kj(P)Ki\  j(P)  +  km(¥)Km(V)}  +  ksc[k/S)Kni(S)  +  kJS)Km(S)} 


Ka2(P)  = 


Kqi(  S)  = 


^P22(H)  = 


(i  -  ^<p)Km(P) + vr^pj^(P) 

\-^k/P)[\-KJP)/Kl22(P)} 

(1  -  yfkj$)Km(S)  +  yfkj$)Km(S) 
l-yfk/S)[l-Km(S)/Kl22(S)] 

Kn2(V) 

1  ~ksc[l  -Ka2(S)/Ka2( P)] 


^03(H)=^G2(H) 


The  void  conductivity  Kv  with  moisture  content  M  is  given  by  Kv  -  MATmst.  If  there  are  voids  in  the 
primary  composite,  Afm(P)  in  the  preceding  equations  for  heat  conductivities  is  replaced  by 


Km(P)  =  [l-^kJP)]Km(P)  + 


'fkjV)Km(Y) 

[1  -  VOP)][l  -Km(P)/K„ 


Similarly,  for  the  secondary  composite,  A^S)  is  replaced  by 


Km(S)  =  [1  -  \fkJS)]Km(S)  - 


VMS)^m(S) 

[1  -  VT7S)][1  -Km(S)/Kv] 


Strengths. — The  longitudinal  strengths  are  based  on  the  rule  of  mixtures: 


Sfl  i  r(H)  -  i  t<P)(1  -  ksc )  +  S#  i  j{S)kSt 


Sn  j  C(H)  =  Sn  j  C(P)(  1  -  ksc)  +  Sn  j  c(S)kSi 


The  following  are  a  few  intermediate  variables  defined  for  convenience  in  the  evaluation  of 
transverse  strengths: 


SmC=rnin[SmC(P)  and  SmC(S)] 

SOT7-=min[Sw7<P)  and  Sm7<S)] 

SmS=rnin[SmS(P)  and  SmS(S)] 

FACT  1  =  km{ P)  (for  HMS  and  Kevlar  fibers) 
FACT  2  =  km(S)  (for  HMS  and  Kevlar  fibers) 
FACT  1  =  FACT  2  =  1  (for  all  other  fibers) 
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The  transverse  and  the  shear  strengths  of  hybrid  composites  are  given  by 


Se 22t<H)  =££22(H) 


I  (1  -Ar^FACTl  kscFACT2 

UG2(P)DENOMP  £m(S)DENOMS  J  " 


d-^c) 


Saic{K)  -  Ea2^)  [_£G2(P)DEN0MP  +  £G2(S)DENOMS 


Sfl2(H)  = 


(1  -Ar5C)gg 
2  Gni\  Gfl2(P) 


^ scQs  P  ACT2  1 

Gn2(  S)  L 


FACT1[,-V¥(-^)] 

-V¥(*-SS)]| 

- \smS 

/  ZO\  \  1 


L_  ..x 


Moisture  diffusivity: 


l-V^P)  +  Ar„(P) 

D'l  P)=  - r~ - : - 

[l->fkjP)]D£P) 

,o  _  i-Vms)+ms) 

T>fi  i(H)  =  (1  -  ksc)km(V)DH?)  +  ksckm(S)D'i(S) 

Dai{ H)  =  (1  - ksc) [l  -  2-Jk/F)]D'^P)  +  -  VtyS)]z>KS) 


£>G3(H)=£>G2(H) 


Moisture  expansion  coefficients: 


^m(P)  = 


0m(S)  = 


km(P)(3JP)Em(P) 

'  Ay(JP)E/11(P)  +  A:/„(P)^m(P) 

Arw(S)/3w(S)Em(S) 
kf(S)Efn(S)  +  km(S)Em(S) 


Pn  2(P)=i8m(P) 


[i-v^p)],; 


A:/<P)A:/w(P)^ni(P) 

1  +  %i(P)  +  Arw(P)[£m(P)-£/11(P)] . 
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/3m(S)  =  MS)[l-^<S)] 


fcAS)km(S)Eni(S) 


Ef\ i  (S)  +  km(S)\Em(S) 


(S) _ ] 

-Em(S)]  ) 


[d  -  ksc)km{V)Pm{V)Em{V)  +  ksckm(S)(3m(S)Em(S)\ 

£m(H) 

/353(H)  =  [-^,3(H)£f33(H)/3fll(H)  +  (1  -  ^c)£'G3(P)[^e2(P)  +  *WP)/3fll(P)]  + 

+  ,fl3(S)/3fll(S)]jAo3(H) 

«'fl2(P)  =  "m(P)  +  ^P)[^2(P)  -  "m(P)] 
j/fl2(S)  =  J'm(S)  +  k/S)[vfi2(S)  ~  ^(S)] 
vB2( P)  =  ^(P)  +  ^P)[^3(P)  -  "m(P)] 

jrf)  =  ^(S)  +  A:/S)[^23(S)  -  i'm(S)] 

/3  (22(H)  =  (1  -  ksc)  [/3G2(P )0  +  */02(P))  +  ^12(P)^fl  l(P)]  +  ^sc[/3f22(s)(1  +  "B2(S)) 

+  f'fH^iSfliCS)]  -  J'fl2(H)/3fli(H)  -  ^53(H)/353(H) 

Flexural  moduli: 

£ai/*H)  =  £m(H) 

^G2f(H)  =i?52(H) 


Flexural  strengths: 


WH)  = 


2Gm(H))  °G3(P) 


d-*»)Q»  r,_-\/^<e71-  52©) 

I  Gm(P)  L  X  x  V  Gej(P)/J 


cG3(S) 


t-V¥(' 


-VtyS)(i 


M s) 

gG3(S)> 


Gm(S) 

G/23(S)  / 


Sfl2S/?(H)  =  1.5Sa2(H) 
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Fiber  volume  ratio: 


*/■  H)  =  k/P)  +  ksc[kj(  S)  -  kj{  P)] 


Subroutines  BANKRD  and  IDGER . — These  two  subroutines  do  preprocessing  to  generate 
compatible  input  data  to  the  subroutine  INHYD.  The  subroutine  BANKRD  is  called  first  by  the 
ICAN  main  program.  The  input  to  this  routine  is  primarily  the  data  supplied  on  the  material  card 
MATCRD  by  the  user.  These  cards  indicate  the  coded  names  for  the  fiber  and  matrix,  the  volume 
ratios  of  primary  and  secondary  composites,  and  their  respective  fiber,  and  the  matrix  and  void 
volume  ratios.  The  subroutine  BANKRD  has  its  own  data  base  containing  the  properties  of  fibers 
and  matrices  of  commonly  used  materials.  This  data  base  is  assigned  to  input  unit  8.  It  is  named 
FBMTDATA.BANK.  The  output  of  BANKRD  are  the  arrays  PFP,  PFS,  PMP,  and  PMS.  The 
entries  in  PFP  and  PFS  are  the  fiber  properties  of  primary  and  secondary  composites.  The  entries  in 
PMP  and  PMS  are  the  matrix  properties  of  primary  and  secondary  composites.  These  arrays  are 
made  common  to  the  main  program  and  the  subroutine  IDGER  through  the  labeled  common  block 
MFBANK.  The  entries  of  PF  and  PM  arrays  are  explained  in  the  following  list: 

Fiber  Property  Arrays  PFP  and  PFS 
Entry  Description  Notation 


1 

not  used 

.... 

2 

fiber  density 

pf 

3 

normal  moduli 

Em 

4 

normal  moduli 

Ej22 

5 

Poisson’s  ratio 

vf\2 

6 

Poisson’s  ratio 

vf2  3 

7 

shear  moduli 

Gf\l 

8 

shear  moduli 

G/23 

9 

thermal  expansion  coefficient 

«/ll 

10 

thermal  expansion  coefficient 

am 

11 

heat  conductivity 

Km 

12 

heat  conductivity 

Kfn 

13 

heat  capacity 

cf 

14 

strengths 

Sfr 

15 

strengths 

Sfc 

16 

not  used 

17 

not  used 

— 

18 

not  used 

— 

19 

not  used 

— 

20 

number  of  fibers  per  end 

Nf 

21 

fiber  diameters 

df 
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Matrix  Property  Arrays  PMP  and  PMS 


Entry 

Description 

Notation 

1 

not  used 

— 

2 

density 

Pm 

3 

normal  modulus 

Em 

4 

Poisson’s  ratio 

vm 

5 

coefficient  of  thermal  expansion 

am 

6 

heat  conductivity 

Km 

7 

heat  capacity 

c 

8 

tensile  strength 

SmT 

9 

compressive  strength 

SmC 

10 

shear  strength 

SmS 

11 

allowable  tensile  strain 

emT 

12 

allowable  compressive  strain 

emC 

13 

allowable  shear  strain 

emS 

14 

allowable  torsional  strain 

emTOR 

15 

void  conductivity 

Kv 

16 

glass  transition  temperature 

Tgdr 

The  coded  names  for  the  fiber  and  matrix  are  stored  in  the  matrix  CODES  by  the  main  program. 
The  entries  in  CODES  are  explained  as  follows: 


CODES(l,l,I) 

CODES(l,2,I) 

CODES(2,l,I) 

CODES(2,2,I) 


coded  name  of  primary  fiber 
coded  name  of  primary  matrix 
coded  name  of  secondary  fiber 
coded  name  of  secondary  matrix 


The  subroutine  IDGER  takes  the  information  generated  by  BANKRD  and  arranges  it  in  a  proper 
format  for  the  subroutine  INHYD.  These  data  are  transferred  to  input  unit  7  prior  to  calling 
INHYD.  These  data  are  purged  at  the  end  of  the  program  execution. 


Data  Base  FBMTDATA.BANK. 

The  constituent  properties  data  base  is  a  unique  feature  of  the  computer  code  ICAN.  Its  primary 
aim  is  to  reduce  the  burden  on  the  user  by  preparing  properly  formatted  data  for  the  program.  The 
user  only  needs  to  specify  the  coded  names  for  the  fiber  and  matrix.  The  format  of  the  data  is 
explained  in  this  section  so  as  to  enable  the  user  to  introduce  new  contents  or  to  modify  existing 
entries  as  appropriate  to  his/her  needs.  Data  for  four  fibers  and  three  matrices  are  provided  in  the 
present  package. 

The  fiber  properties  are  arranged  in  five  physical  cards  of  80  column  length.  The  first  card  contains 
a  four-character  code  name  of  a  fiber  in  format  A4.  The  second  to  the  fifth  cards  start  with  a  two- 
letter  mnemonic  to  indicate  the  type  of  properties  that  follow.  The  format  on  any  of  these  cards  is 
A4,  7E10.3,  except  for  the  second  card.  The  second  card  is  in  the  format  A3,  16,  7E10.3.  The 
mnemonics  FP,  FE,  FT,  and  FS  stand  for  fiber  physical,  elastic,  thermal,  and  strength-related 
properties,  respectively.  The  entries  on  these  cards  are  explained  as  follows: 


card  1  four  character  coded  name  for  fiber 
card  2  FP;  Nf,  djy  pj 

card  3  FE;  Eph  Ep2,  yp2,  Y/23>  Gf\2>  GJ23 
card  4  FT;  oljw,  otp2,  Kf\\*  Kfih  Gf 

card  5  FS;  SjTi  SfC  (The  remaining  entries  are  open  for  future  modifications.) 

The  matrix  properties  are  arranged  next  after  the  line  OVER  END  OF  FIBER  PROPERTIES.  The 
properties  have  essentially  the  same  format  as  those  for  fiber  property  cards.  There  are,  however,  six 
physical  cards  for  each  matrix  material.  The  mnemonics  used  are  MP,  ME,  MT,  MS,  and  MV.  They 
stand  for  matrix  physical,  elastic,  thermal,  strength-related,  and  miscellaneous  properties, 
respectively.  The  format  for  the  first  card  is  A4,  and  the  format  for  the  rest  of  the  cards  is  A3, 
7E10.3.  The  entries  in  each  card  are  as  follows: 

card  1  four  character  coded  name  for  matrix 

card  2  MP;  pm 

card  3  ME;  Em,  vm,  am 

card  4  MT ;  Km,  Cm 

card  5  MS,  Smj,  Smc>  ^mS>  ^mT>  emO  e mS>  e mTOR 

card  6  MV;  Kv>  Tgdr 

The  data  base  presently  contains  properties  for  T-300  (T300),  AS  graphite  (AS--),  S-Glass  (SGLA), 
and  HMS  (HMSF)  fibers.  The  available  matrix  materials  are  high-modulus,  high-strength  (HMHS), 
intermediate-modulus,  high-strength  (IMHS),  and  intermediate-modulus,  low-strength .  (IMLS) 
matrices,  which  are  epoxy-type  resins.  The  complete  list  of  properties  is  shown  in  appendix  C. 


Lewis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  October  11,  1985 
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Appendix  A 
List  of  Code  Identifiers 


Engineering  Fortran  Comment 

symbol  symbol  code 


■^cx 

ACX 

composite  axial  stiffness;  generated  in  subroutine 
GPCFD2 

RAC 

reduced  axial  stiffness;  computed  in  subroutine  GPCFD2 

BIDE 

Boolean 

true  if  interply  effects  are  included;  input 

Ccc 

CPC 

composite  coupling  stiffness;  generated  in  subroutine 
GPCFD2 

Cel 

RESF 

string  with  force  variables  in  BLOCK  DATA 

Ce2 

DISP 

string  with  displacement  variables  in  BLOCK  DATA 

COMSAT 

Boolean 

true  if  COMSA  is  executed;  input 

CSANB 

Boolean 

true  if  membrane  and  bending  symmetry  exists;  input 

Ecx 

FTC 

composite  flexural  rigidities;  generated  in  subroutine 
GPCFD2 

D§x 

RDC 

reduced  bending  rigidities;  computed  in  subroutine 
GPCFD2 

Df 

DIAF 

filament  equivalent  diameter;  input 

Dv 

DISV,  DISVI 

displacement  vectors;  DISVI  is  either  read  in  main 
program,  or  is  generated  in  subroutine  COMSA 

Ef>  Ecf 

ECF 

filament  elastic  constants;  input 

Ef\\,a\,mn 

EF11,EL11,EM11 

filament,  ply,  and  matrix  normal  moduli;  filament  and 
matrix  moduli  input 

Gfl2,n2,mU 

EF12,EL12,EM12 

filament,  ply,  and  matrix  shear  moduli;  filament  and 
matrix  shear  moduli  input 

E,EC( 

ECL 

ply  elastic  constants;  generated  in  subroutine  INHYD 

Em>Ecm 

ECM 

matrix  elastic  constants;  generated  in  subroutine  INHYD 

Hj 

PL(9,I) 

interply  distortion  energy  coefficient;  generated  in  main 
program 

Hkc 

CHK 

array  of  constituent  heat  conductivities;  input 

hc 

HHC 

composite  heat  capacity  stored  in  PC(18)  and  PC(54) 

u 

IJ 

index;  generally  ply  or  interply 

Ecxx,cyy,cxy 

HK1 1,22,33 

composite  two-dimensional  heat  conductivities  in  PC(51) 
to  PC(53) 

Kc\\,c22,c22 

HK1 1,22,33 

composite  three-dimensional  heat  conductivities  along 
the  material  axes  in,  PC(15)  to  PC(17) 

Kf,y 

KF,V 

apparent  fiber  and  void  volume  ratios;  input 

CHK 

see  Hkc 

E\xx,  1  yy,  Ixy 

XK1,XK2,XK3 

stress  concentration  factors  generated  in  STRCNF 

kf,m 

KFB,MB 

actual  fiber  and  matrix  volume  ratios 

kf!,vl 

KFL,VL 

ply  apparent  fiber  and  void  volume  ratios 

Esc 

LSC 

array  of  limiting  conditions;  input 

Mcx 

MBS 

applied  moment;  input 

McTtx 

MSDT 

thermal  moments;  generated  in  GPCFD2 

McM(x 

MSDH 

hygral  moments;  generated  in  GPCFD2 

m 

M 

load  condition  index 

E^cx 

NBS 

applied  membrane  loads;  input 

E^cMfX 

NSDH 

hygral  force;  generated  in  GPCFD2 

NcTtx 

NSDT 

thermal  force;  generated  in  GPCFD2 

Nf 

NFPE 

number  of  filaments  per  end;  input 

N( 

NL 

number  of  plies;  input 

Nt 

NLC 

number  of  load  conditions;  input 

N 

iyms 

NMS 

number  of  material  systems;  input 

Npc 

NPC 

string  PROPC  length;  input 

"pi 

NPL 

string  PROP  length;  input 

NONUDF 

Boolean 

T  (true)  if  Poisson’s  ratio  difference  chart  is  to  be 
suppressed 

Pc 

PC 

composite  properties  array;  generated  in  GACD3  and 
GPCFD2 

Pep 

PROPC 

string  PROPC;  composite  property  identifiers  in 
GDCFD2 

Pi 

PL 

ply  property  array;  portions  generated  in  all  parts  of  the 
program 

P(p 

PROP 

string  PROP;  ply  properties  identifiers  in  main  program 

Qf,i,p,r,s 

QF,I,P,R,S 

indices  to  print  out  string  PROP 

R 

R 

transformation  matrix;  GACD3,  GPCFD2,  and  COMSA 

RINDV 

Boolean 

T  (true)  if  displacements  are  read  in;  input 

^flir>etc- 

PL(51)  to  PL(59,I) 

ply  limit  stresses;  generated  in  GLLSC 

h 

TL 

ply  thickness;  input 

Web 

wxx 

composite  local  curvatures  relative  to  the  structural  axes 

0ic 

CTE 

composite  coefficient  of  thermal  expansion;  three- 
dimensional  in  PC(12)  to  PC(14),  two-dimensional  in 
PC(48)  to  PC(50) 

afXm 

vaf,al,am 

filament,  ply,  and  matrix  thermal  coefficients  of 
expansion;  input 

$t,vt 

VCF 

correlation  factors  for  ply  thermoelastic  properties  and 
strain  magnification  factors;  set  to  unity  in  COMSA 

$h 

BTA 

correlation  factors  for  ply  heat  conductivity;  set  to  unity 
in  COMSA 

& 

BET 

correlation  factors  for  ply  strength;  set  to  unity  in 
COMSA 

5f 

PL(8,I) 

interply  layer  thickness;  generated  in  INHYD 

ecsx 

UX 

reference  plane  membrane  strain;  solved  in  terms  of  Ncx 
or  input 

EPS,PL(74)  to  PL(66,I) 

ply  strains;  generated  in  COMSA 

@cs 

THCS 

angle  between  composite  material  and  structural  axes; 
input 

9a,dk 

THLC 

angle  between  ply  material  and  composite  axes;  input 

vf\2,  fl2,  m\2 

NUF12,L12,M12 

filament,  ply,  and  matrix  Poisson’s  ratios;  input 

7T 

PIE 

constant;  input 

Pf.m.i 

RHOF,M,L 

filament  and  matrix  weight  density;  input  and  generated 
in  FIBMT,  HTM,  and  COMPP 

Of,  om 

SF,  SM 

microstresses  in  fibers  and  matrices  generated  in 

MCRSTR 

k 

XPL,XSL,YPL,YSL 

boundary  zone  decay  length;  generated  in  the  main 
program  and  paired  to  EDGSTR 

SIGL,PL(67)  to  PL(69,I) 

ply  stress;  generated  in  COMSA 

Appendix  B 
Sample  Input/Output 


Item  1 


xxxxxxxx 

xxxxxxxx 

XX 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


XX 

SSSSSSSSSS 

ccccccccc 

AAAAAAAA 

NNN 

NNN 

KX 

XX 

SSSSSSSSS5SS 

ccccccccccc 

A AAAAAA AA A 

NNNN 

NNN 

XX 

XX 

SSS5 

ccc 

AAA 

AAA 

NNNNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNNNNN 

NNN 

XX 

XX 

SS5S 

ccc 

AAA 

AAA 

NNNHNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA AAAAAAAA A 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA AAAAAAAA A 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNNNNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNNNNN 

XX 

XX 

SSSS 

ccc 

AAA 

AAA 

NNN 

NNNNN 

XX 

XX 

SS5SSSSSSSSS 

ccccccccccc 

AAA 

AAA 

NNN 

NNNN 

XX 

XX 

SSSSSSSSSS 

ccccccccc 

AAA 

AAA 

NNN 

NNN 

XX 

XX 

XX 

XX 

XX 

XX 

INTEGRATED 

COMPOSITES  ANALYZER 

BY  P 

.  L  .  N  . 

MURTHY 

XX 

XX 

NASA 

LERC. CLEVELAND 

XX 

XX 

SEPTEMBER 

1983 

XX 

xxxxxxxx 

xxxxxxxx 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


Item  2 


I  C  A  N :  COORDINATE  SYSTEMS 


/  /I 

/  /A  | 

/  /ABj 

/  /ABB| 

/.  /ABBA/ 

/.  .  /ABBA/ 

/.  .  .  THETA  /ABBA/ 

/.  .  .  .  /ABBA/ 

/ .  /ABBA/ 

- - - ABBA/ 

|  000000000000000000000000000000 ! BBA/ 

I XXXXXXXXXXXXXXXXXXXXXXXXXXXXXX i BA/ 

I XXXXXXXXXXXXXXXXXXXXXXXXXXXXXX ! A/ 

1 000000000000000000000000000000  j/ 


/  /I 

/  /B/- 

/  /B/ 

/  /B/ 

Y  /  /B/ 

/  /B/ 

/  /B/ 

/  /B/ 

/  /B/ 

- B/ 

!  000000000000000000000000000000 | / 


LAMINATE  STRUCTURAL  AXES 


PLY  MATERIAL  AXES 


Item  3 


I  C  A  H 

INPUT  DATA  ECHO 


FOUR  PLY  SYMMETRIC  LAMINATE.  ICAN  SAMPLE  INPUT  DATA. 


STDATA 

4 

1 

2 

T 

COMSAT 

F 

CSANB 

F 

BIDE 

F 

RINDV 

T 

NONUDF 

PLY 

1 

1 

70.00 

70.0 

.0  0.0 

.010 

PLY 

2 

2 

70.00 

70.0 

.0  90.0 

.  005 

PLY 

3 

2 

70.00 

70.0 

.  0  90.0 

.005 

PLY 

4 

1 

70.00 

70.0 

.0  0.0 

.010 

MATCRDAS — IMLS 

.55 

.02 

AS— IMLS 

0.0  .57 

.03 

MATCRDSGL AHMHS 

.55 

.01 

AS— IMHS 

0.4  .57 

.01 

PLOAD 

1000  . 

0.0 

0.0 

0  .  0 

NX,  NY/ 

NXY , THCS 

PLOAD 

0.0 

0.0 

0.0 

MX , MY  f  1 

MXY 

PLOAD 

0.0 

0.0 

DMX/QX 

,DMY/QY,PRSS 

Item  4 

J  M  M  A 

R  Y  0  F 

I  N  P 

U  T  D 

A  T  A 

FOUR  PLY  SYMMETRIC  LAMINATE.  ICAN  SAMPLE  INPUT  DATA. 


NUMBER 

OF 

CASE  CONTROL 

LAYERS 

DECK 

NL  = 

4 

NUMBER 

OF 

LOADING 

CONDITIONS  NLC  = 

1 

NUMBER 

OF 

MATERIAL 

SYSTEMS 

NMS  = 

2 

COMSAT 

T 

CSANB 

F 

BIDE 

F 

RINDV 

F 

NONUDF 

T 

_  „  _ 

LAMINATE 

CONFIGURATION 

-  -  - 

PLY 

NO 

MID 

DELTAT 

DELTAM 

THETA 

T-NESS 

PLY 

1 

1 

0.000 

0.0% 

0.0 

0.010 

PLY 

2 

2 

0.000 

0.0% 

90.0 

0.005 

PLY 

3 

2 

0.000 

0  .  0% 

90 . 0 

0.005 

PLY 

4 

1 

0  .000 

0  .  0% 

0.0 

0  .010 

-  -  - 

COMPOSITE  MATERIAL  SYSTEMS 

■  -  - 

MATCRD 

MID 

PRIMARY 

VFP 

VVP 

SECONDARY 

VSC 

VFS 

vvs 

MATCRD 

1 

AS  — IMLS 

0.55 

0.02 

AS  — IMLS 

0.00 

0.57 

0.03 

MATCRD 

2 

SGLAHMHS 

0.55 

0  .01 

AS  — IMHS 

0.40 

0.57 

0.01 

LOADING  CONDITIONS 


PRESCRIBED  LOADS  FOR 

THE  LOAD 

CONDITION  1 

INPLANE  LOADS 

NX 

= 

1000.0000 

LB/IN 

NY 

0.00  00 

LB/IN 

NXY 

= 

0.0000 

LB/IN 

3ENDING  LOADS 

MX 

= 

0.0000 

LB. IN/IN 

MY 

= 

0.0000 

LB. IN/IN 

MXY 

= 

0.0000 

LB. IN/IN 

TRANSVERSE  LOADS 

DMX/QX 

s 

0  .  00  00 

LB/IN 

DMY/QY 

= 

0  .0000 

LB/IN 

TRANSVERSE  PRESSURE 

PU 

= 

0.0  000 

LB/SQ.  IN 

TRANSVERSE  PRESSURE 

PL 

0  .0000 

LB/SQ.  IN 

56 


Item  5(a) 

— >  CONSTITUENT  PROPERTIES:  ECHO  FROM  DATA  BANK.  <-- 


PRIMARY  FIBER  PROPERTIES;  AS 

1  ELASTIC  MODULI 

2 

3  SHEAR  MODULI 

4 

5  POISSON’S  RATIO 

6 

7  THERM.  EXP.  COEF. 

8 

9  DENSITY 

10  HO.  OF  FIBERS/END 

11  FIBER  DIAMETER 

12  HEAT  CAPACITY 

13  HEAT  CONDUCTIVITY 

14 

15 

16  STRENGTHS 

17 


PRIMARY  MATRIX  PROPERTIES; 

1  ELASTIC  MODULUS 

2  SHEAR  MODULUS 

3  POISSON’S  RATIO 

4  THERM.  EXP.  COEF. 

5  DENSITY 

6  HEAT  CAPACITY 

7  HEAT  CONDUCTIVITY 

8  STRENGTHS 

9 

10 

11  MOISTURE  COEF 

12  DIFFUSIVITY 


—  FIBER 

EFP1  0.3100E  08 
EFP2  0.2000E  07 
GFP12  0.2000E  07 
GFP23  0.1000E  07 
NUFP12  0.2000E  00 
NUFP23  0.2500E  00 
CTEFP1  -0.5500E-06 
CTEFP2  0.5600E-05 
RHOFP  0.6300E-01 
NFP  0.1000E  05 
DIFP  0.30  0  0E-03 
CFPC  0.1700E  ao 
KFP1  0.5800E  03 
KFP2  0.5800E  02 
KFP3  0.5800E  02 
SFPT  0.4000E  06 
SFPC  0.4000E  06 


IMLS  MATRIX.  DRY  RT .  PROPERTIES. 

EMP  0.5000E  06 
GMP  0.1773E  06 
HUMP  0.4100E  00 
CTEMP  0 . 57  Q0E-G4 
RHOMP  0.4600E-01 
CMPC  0.2500E  00 
KMP  0.1250E  01 
SMPT  0.7000E  04 
SHPC  0.2100E  05 
SMPS  0.7000E  04 
BTAMP  0.4000E-02 
DIFMP  0.2000E-03 


Item  5(b) 

PRIMARY  COMPOSITE  PROPERTIES;  55/  43  AS — /IMLS 

BASED  ON  MICROMECHANICS  OF  INTRAPLY  HYBRID  COMPOSITES:  ELASTIC  AND  THERMAL  PROPERTIES. 


FIBER 

VOLUME  RATIO  -  0.550 

MATRIX 

VOLUME  RATIO 

VOID 

CONDUCTIVITY  -  0.22499990E  00 

1 

ELASTIC  MODULI 

EPC1 

0.1726E  08 

2 

EPC2 

0.1127E  07 

3 

EPC3 

Q.1127E  07 

4 

SHEAR  MODULI 

GPC12 

0.5470E  06 

5 

GPC23 

0.3238E  06 

6 

GPC13 

0.5470E  06 

7 

POISSON’S  RATIO 

NUPC12 

0.2945E  00 

8 

NUPC23 

0 . 482 IE  00 

9 

NUPC13 

0.2945E  00 

10 

THERM.  EXP.  COEF. 

CTEPC1 

0.1418E-06 

11 

CTEPC2 

0 .2464E-04 

12 

CTEPC3 

0.2464E-04 

13 

DENSITY 

RHOPC 

0 . 544 3E“0 1 

14 

HEAT  CAPACITY 

CPC 

0.1991E  00 

15 

HEAT  CONDUCTIVITY 

KPC1 

0.3195E  03 

16 

KPC2 

0.3702E  01 

17 

KPC3 

0.3702E  01 

18 

STRENGTHS 

SPC1T 

0.2228E  06 

19 

SFC1C 

0.8764E  05 

20 

SPC2T 

0.5006E  04 

21 

SPC2C 

0.1502E  05 

22 

SPC12 

0.5126E  04 

23 

MOIST.  DIFFUSIVITY 

DPC1 

0 .8600E-04 

24 

DPC2 

0 .5168E-04 

25 

DPC3 

0.5168E-04 

26 

MOIST.  EXP.  COEF. 

BTAPC1 

0 . 4  98  IE- 04 

27 

BTAPC2 

0.1452E-02 

23 

BTAFC3 

0.1452E-02 

29 

FLEXURAL  MODULI 

EPC1F 

0.1726E  08 

30 

EPC2F 

0.1127E  07 

31 

STRENGTHS 

SFC23 

0.3983E  04 

32 

SPC1F 

0.1572E  06 

33 

SPC2F 

0.9387E  04 

34 

SPCSB 

0.7689E  04 

35 

PLY  THICKNESS 

TPC 

0.5000E-02 

36 

INTERPLY  THICKNESS 

PLPC 

0.5850E-04 

37 

INTERFIBER  SPACING 

PLPCS 

0 .5850E-04 

VOID  VOLUME  RATIO  -  0.020 
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Item  5(c) 

—  >  CONSTITUENT  PROPERTIES  *  ECHO  FROM  DATA  BANK.  <— 
PRIMARY  FIBER  PROPERTIES;  SGL A  FIBER 


1 

ELASTIC  MODULI 

EFP 1 

0.1240E  08 

2 

EFP2 

0.1240E  08 

3 

SHEAR  MODULI 

GFP12 

0.5170E  07 

4 

GFP23 

0.5170E  07 

5 

POISSON* S  RATIO 

NUFP12 

0.2000E  00 

6 

NUFP23 

0.2000E  00 

7 

THERM.  EXP.  COEF. 

CTEFP1 

0.280  0E-05 

8 

CTEFP2 

0.28C0E-05 

9 

DENSITY 

RHOFP 

0.9000E-01 

10 

NO.  OF  FIBERS/END 

NFP 

0.2040E  03 

11 

FIBER  DIAMETER 

DIFP 

0.3600E-03 

12 

HEAT  CAPACITY 

CFPC 

0.1700E  00 

13 

HEAT  CONDUCTIVITY 

KFP1 

0.7500E  01 

14 

KFP2 

0.7500E  01 

15 

KFP3 

0.7500E  01 

16 

STRENGTHS 

SFPT 

0.3600E  06 

17 

SFPC 

0.3000E  06 

PRIMARY 

MATRIX  PROPERTIES; 

HMHS  MATRIX. 

DRY  RT.  PROPERTIES 

1 

ELASTIC  MODULUS 

EMP 

0.7500E  06 

2 

SHEAR  MODULUS 

GMP 

0.2778E  06 

3 

POISSON  *  S  RATIO 

NUMP 

0.3500E  00 

4 

THERM.  EXP.  COEF. 

CTEMP 

0.4000E-04 

5 

DENSITY 

RHOMP 

0.4500E-01 

6 

HEAT  CAPACITY 

CMPC 

0.2500E  00 

7 

HEAT  CONDUCTIVITY 

KMP 

0 . 1250 E  01 

8 

STRENGTHS 

SMPT 

0.2000E  05 

9 

SMPC 

0.5000E  05 

10 

SMPS 

0.1500E  05 

11 

MOISTURE  COEF 

BTAMP 

0.4000E-02 

12 

DIFFUSI VITY 

DIFMP 

0.2000E-03 

Item  5(d) 

PRIMARY 

COMPOSITE  PROPERTIES; 

55/  44  SGLA/HMHS 

BASED  ON  MICROMECHANICS  OF  INTRAPLY  HYBRID  COMPOSITES:  ELASTIC  AND  THERMAL  PROPERTIES. 

FIBER  VOLUME  RATIO  -  0.550  MATRIX  VOLUME  RATIO  -  0,440  VOID  VOLUME  RATIO  -  0.010 
VOID  CONDUCTIVITY  -  0.22499990E  00 


1 

ELASTIC  MODULI 

EPC1 

0.7150E  07 

2 

EPC2 

0.2473E  07 

3 

EPC3 

0.2473E  07 

4 

SHEAR  MODULI 

GPC12 

0.9314E  06 

5 

GPC23 

0.5792E  06 

6 

GPC13 

0.9314E  06 

7 

POISSON'S  RATIO 

NUPC12 

0.2675E  00 

3 

NUPC23 

0.3778E  00 

9 

NUPC13 

0.2675E  00 

10 

THERM.  EXP.  COEF. 

CTEPC1 

0.4488E-05 

11 

CTEPC2 

0.1580E-04 

12 

CTEPC3 

0.1580E-04 

13 

DENSITY 

RHOPC 

0 . 6  930E-0 1 

14 

HEAT  CAPACITY 

CPC 

0.1929E  00 

15 

HEAT  CONDUCTIVITY 

KPC1 

0.4675E  01 

16 

KPC2 

0.2750E  01 

17 

KPC3 

0.2750E  01 

18 

STRENGTHS 

SPC1T 

0.2076E  06 

19 

SPC1C 

0.1730E  06 

20 

SPC2T 

0.1256E  05 

21 

SPC2C 

0.3140E  05 

22 

SPC12 

0.1047E  05 

23 

MOIST.  DIFFUSIVITY 

DPC1 

0.8800E-04 

24 

DPC2 

0.5168E-04 

25 

DPC3 

0.5168E-04 

26 

MOIST.  EXP.  COEF. 

BTAPC1 

0.1846E-03 

27 

BTAPC2 

0.1379E-02 

28 

BTAPC3 

0.1379E-02 

29 

FLEXURAL  MODULI 

EPC1F 

0.7150E  07 

30 

EPC2F 

0.2473E  07 

31 

STRENGTHS 

SPC23 

0.6510E  04 

32 

SPC1F 

0.2359E  06 

33 

SPC2F 

0.2243E  05 

34 

SPCSB 

0.1570E  05 

35 

PLY  THICKNESS 

TPC 

0. 5000 E” 02 

36 

INTERPLY  THICKNESS 

PLPC 

0. 7020 E- 04 

37 

INTERFIBER  SPACING 

PLPCS 

0.7020E-04 
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Item  5(e) 

— >  CONSTITUENT  PROPERTIES:  ECHO  FROM  DATA  BANK. 
SECONDARY  FIBER  PROPERTIES;  AS—  FIBER 


1 

ELASTIC  MODULI 

EFS1 

0.310CE  05 

2 

EFS2 

0.2000E  07 

3 

SHEAR  MODULI 

GFS12 

0.2000E  07 

4 

GFS23 

0.1000E  07 

5 

POISSON’S  RATIO 

NUFS12 

0.2000E  00 

6 

NUFS23 

0 . 250  0E  0  0 

7 

THERM.  EXP.  COEF. 

CTEFS1 

-Q.5500E-06 

8 

CTEFS2 

0 . 5600E-05 

9 

DENSITY 

RHOFS 

0.6300E-01 

10 

NO.  OF  FIBERS/END 

NFS 

0.1000E  05 

11 

FIBER  DIAMETER 

DIFS 

0.300  0E" 0  3 

12 

HEAT  CAPACITY 

CFSC 

0.1700E  00 

13 

HEAT  CONDUCTIVITY 

KFS1 

0.5800E  03 

1  4 

KFS2 

0.5800E  02 

15 

KFS3 

0.5800E  02 

16 

STRENGTHS 

SFST 

0.4000E  06 

17 

SFSC 

0.4000E  06 

SECONDARY 

MATRIX  PROPERTIES; 

IMHS  MATRIX 

:.  DRY  RT.  3 

1 

ELASTIC  MODULUS 

EMS 

0.5000E  06 

2 

SHEAR  MODULUS 

GMS 

0.1852E  06 

3 

POISSON’S  RATIO 

HUMS 

0 . 350 0E  00 

4 

THERM.  EXP.  COEF. 

CTEMS 

0 . 36Q0E-Q4 

5 

DENSITY 

RHOMS 

0 .4400E-01 

6 

HEAT  CAPACITY 

CM  SC 

0.2500E  00 

7 

HEAT  CONDUCTIVITY 

KMS 

0.1250E  01 

8 

STRENGTHS 

SMST 

0.1500E  05 

9 

SMSC 

0.3500E  05 

10 

SMSS 

0.1300E  05 

11 

MOISTURE  COEF 

BTAMS 

0 .4000E-02 

12 

DIFFUSIVITY 

DIFMS 

0 .200  0E-03 

Item  5(f) 

SECONDARY  COMPOSITE  PROPERTIES;  57/  42  AS — /IMHS 

BASED  ON  MICROMECHANICS  OF  INTRAFLY  HYBRID  COMPOSITES:  ELASTIC  AND  THERMAL  PROPERTIES. 


FIBER  VOLUME  RATIO  -  0.570  MATRIX  VOLUME  RATIO  -  0.420  VOID  VOLUME  RATIO  -  0.010 

VOID  CONDUCTIVITY  -  0.22499990E  00 


1 

ELASTIC  MODULI 

ESC1 

0.1788E  08 

2 

ESC2 

0.1153E  07 

3 

ESC3 

0.1153E  07 

4 

SHEAR  MODULI 

GSC12 

0 . 588 0E  06 

5 

GSC23 

0.3458E  06 

6 

GSC13 

0.5880E  06 

7 

POISSON’S  RATIO 

NUSC12 

0.2645E  00 

8 

MUSC23 

0.4294E  00 

9 

NUSC13 

0.2645E  00 

10 

THERM.  EXP.  COEF. 

CTESC1 

-0 . 1280E-06 

11 

CTESC2 

0 . 1605E-04 

12 

CTESC3 

0 . 1605E-04 

13 

DENSITY 

RHOSC 

0 .5439E-01 

14 

HEAT  CAPACITY 

CSC 

0.1972E  00 

15 

HEAT  CONDUCTIVITY 

KSC1 

0 . 331  IE  03 

16 

KSC2 

0.3918E  01 

17 

KSC3 

0.3918E  01 

18 

STRENGTHS 

SSC1T 

0.2307E  06 

19 

SSC1C 

0.1568E  06 

20 

SSC2T 

0.1026E  05 

21 

SSC2C 

0.2394E  05 

22 

SSC12 

0.9369E  04 

23 

MOIST.  DIFFUSIVITY 

DSC1 

0 .8400E-04 

24 

DSC2 

0.4900E-04 

25 

DSC3 

0 .4900E-04 

26 

MOIST.  EXP.  COEF. 

BTASC1 

0 .4698E-04 

27 

BTASC2 

0 .1319E-02 

28 

3TASC3 

0 . 1319E-02 

29 

FLEXURAL  MODULI 

ESC  IF 

0.1788E  08 

30 

ESC2F 

0.1153E  07 

31 

STRENGTHS 

SSC23 

0.7424E  04 

32 

SSC1F 

0.2334E  06 

33 

SSC2F 

0.1796E  05 

34 

SSCSB 

0.1405E  05 

35 

PLY  THICKNESS 

TSC 

0 .5000E-02 

36 

INTERPLY  THICKNESS 

PLSC 

0.5215E-04 

37 

INTERFIBER  SPACING 

PLSCS 

0.5215E-04 
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Item  5(g) 

HYBRID  COMPOSITE  PROPERTIES;  60/40  SGLA/HMHS/AS — /IMHS 

BASED  ON  MICROMECHANICS  OF  INTRAPLY  HYBRID  COMPOSITES:  ELASTIC  AND  THERMAL  PROPERTIES. 


PRIMARY  COMPOSITE  VOLUME  RATIO  -  0.600  SECONDARY  COMPOSITE  VOLUME  RATIO  -  0.400 


1 

ELASTIC  MODULI 

EHC1 

0.1144E  08 

2 

EHC2 

0  .  16  96E  07 

3 

EHC3 

0.1945E  07 

4 

SHEAR  MODULI 

GHC12 

0.7551E  06 

5 

GHC23 

0 .  456  IE  06 

6 

GHC13 

0 . 7  94  IE  06 

7 

POISSON'S  RATIO 

NUHC12 

0 . 26  6  3E  00 

8 

NUHC23 

0.3985E  00 

9 

NUHC13 

0.2689E  00 

10 

THERM.  EXP.  COEF. 

CTEHC1 

0.1603E-05 

11 

CTEHC2 

0 .1601E-04 

12 

CTEHC3 

0 . 1634E-04 

13 

DENSITY 

RHOHC 

0 .6334E-01 

14 

HEAT  CAPACITY 

CHC 

0.1943E  00 

15 

HEAT  CONDUCTIVITY 

KHC1 

0.1352E  03 

16 

KHC2 

0.2305E  01 

17 

KHC3 

0.2305E  01 

18 

STRENGTHS 

SHC1T 

0.2168E  06 

19 

SHC1C 

0.1665E  06 

20 

SHC2T 

0.9915E  04 

21 

SHC2C 

0.2314E  05 

22 

SHC12 

0.1195E  05 

23 

MOIST.  DIFFUSIVITY 

DHC1 

0 .8736E-04 

24 

DHC2 

0 .5117E-04 

25 

DPC3 

0 . 51 17E-04 

26 

MOIST.  EXP.  COEF. 

BTAHC1 

0 . 9858E-0  4 

27 

BTAHC2 

0 .8565E-03 

28 

BTAHC3 

0 . 1455E-02 

29 

FLEXURAL  MODULI 

EHC1F 

0.1144E  08 

30 

EHC2F 

0.1696E  07 

31 

STRENGTHS 

SHC23 

0.1019E  05 

32 

SHC1F 

0.2355E  06 

33 

SHC2F 

0.1735E  05 

34 

SHCSB 

0.1793E  05 

35 

PLY  THICKNESS 

THC 

0 .5000E-02 

36 

INTERPLY  THICKNESS 

PLHC 

0.5215E-04 

37 

INTEP.FIBER  SPACING 

PLHCS 

0.5215E-04 

38 

FI3ER  VOL.  RATIO 

VFH 

0 . 558 0E  00 

39 

MOISTURE  CONTENT 

M 

0  .  0  0  00 

40 

MATRIX  VOL.  RATIO 

VMH 

0.4320E  00 
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Item  6 


1 

2 

3 

4 

5 

6 


3-D  COMPOSITE 

STRAIN  STRESS 

TEMPERATURE 

MOISTURE  RELATIONS  -  STRUCTURAL 

AXES 

-1- 

-2- 

-3- 

-4- 

-5- 

-6- 

-DT- 

0 . 6976E-07 

-0 .5952E-08 

-0 . 2727E-07 

0.0000 

0.0000 

0.3255E-13 

0 .1102E-05 

0. 

-0.5952E-08 

0  . 1 962E-06 

-0.8485E-07 

0.0000 

0.0000 

-0.1464E-11 

0 . 5805 E- 05 

0, 

-0 . 2727E-07 

-0.8485E-07 

0. 5614 E- 06 

0.0000 

0.0000 

0 .6682E-12 

0 .2839E-04 

0 

0.0000 

0  .0000 

0.0000 

0 .2139E-05 

0. 5229 E- 12 

0.0000 

0.0000 

0  .  0  000 

0.0000 

0 .0000 

0 .5229E-12 

0 . 1935E-05 

0.0000 

0.0000 

0 . 3255E-1 3 

-0 . 1464E-11 

0 .6682E-12 

0 .0000 

0 .0000 

0.1622E-05 

-0.4330E-10 

-0 

3-D  COMPOSITE  STRE5S  STRAIN  RELATIONS  -  STRUCTURAL  AXES 


-1- 

-2- 

-3- 

-4- 

-5- 

-6- 

1 

0.1473E  08 

0.8093E  06 

0.8377E  06 

0  .000  0 

0.0000 

0 .8960E-01 

2 

0 .8093E  06 

0.5499E  07 

0.8703E  06 

0.0000 

0 .0000 

0.45S7E  01 

3 

0.8377E  06 

0.8703E  06 

0.1953E  07 

0.0000 

0.0000 

-0.3603E-01 

4 

0  .  0000 

0 .0000 

0  .  0  0  0  0 

0.4676E  06 

-0.1263E  00 

0  .  0  0  0  0 

5 

0.0000 

0.0000 

0.0000 

-0.1263E  00 

0.5167E  06 

0.0000 

6 

0 .8960E-01 

0.4587E  01 

-0.3603E-01 

0.0000 

0  .0  000 

0.6164E  06 

MAT 9  CARD  FOR  MSC/NASTRAN  SOLID  ELEMENTS 


ril  GI2  G13  G1 4  » G15  » G16  »  G22»  G23 »  G24  ,  G25 » G26  » G33  *  G34  » G35  »036  » G44»  G45  i.  G46  »  G55 

n  n*  n  8 0  927  925 E  06  0  83770038E  06  0 .39597344E-01  0.00000000 

S'^S65379E  01  0  00000000  0.00000000  0.19533170E  °7-0 . 36027569E-01 

"o. 00000000  0.00000000  0 . 46757519E  06-0 . 12633294E  00  0.51670638E  06 


G56,G66 

0  .  00000000 
0.00000000 


0 . 549876 90E  07  0.87032406E 
0.00000000  0 . 61636813E 


Item  7 

COMPOSITE  PROPERTIES 


rn^pn^TTF~ PROPERTIES” -"VALID*" ONLY  FOR  CONSTANT  TEMPERATURE  THROUGH  THICKNESS 
COMPOSHE  PROPERTIES  pVAUiu  ^pERTIES  ABquT  MATERIAL  AXES 

LINES  33  TO  62  2-D  COMPOSITE  PROPERTIES  AB0U^STRUCTURAL^AXES _ 


1 

RHOC 

0.5740E-01 

2 

TC 

0 . 3000E-01 

3 

ecu 

0.1473E  08 

4 

CC12 

0.8093E  06 

5 

CC1 3 

0.8377E  06 

6 

CC22 

0.5499E  07 

7 

CC23 

0.8703E  06 

8 

CC33 

0.1953E  07 

9 

CC44 

0.4676E  06 

10 

CC55 

0.5167E  06 

11 

CC66 

0.6164E  06 

12 

CTE11 

0 .11Q2E-05 

13 

CTE22 

0 .53Q5E-05 

14 

CTE33 

0.2839E-04 

15 

HK1 1 

0.2138E  03 

16 

HK22 

0.4755E  02 

17 

HK33 

0.3080E  01 

18 

HHC 

0.1975E  00 

19 

ECU 

0.1433E  08 

20 

EC22 

0.5098E  07 

21 

EC33 

Q.1731E  07 

22 

EC23 

0  .  46 7 6 E  06 

23 

EC3 1 

0 .5167E  06 

24 

EC  12 

0.6164E  06 

25 

MUC12 

0.8531E-01 

26 

NUC21 

0  .  3034E-01 

27 

N'JCl  3 

0.3908E  00 

28 

NUC31 

0 .4856E-01 

29 

NUC23 

0.4326E  00 

30 

NUC32 

0.1511E  00 

31 

ZCGC 

0 . 1500E-01 

32 

B2DEC 

0.0000 

33 

CC11 

0.1215E  08 

34 

CC12 

0.3746E  06 

35 

CC 1 3 

0 . 1067E  00 

36 

CC22 

0.4610E  07 

37 

CC23 

0.4055E  01 

38 

CC33 

0.6164E  06 

39 

ECU 

0.1212E  08 

40 

EC22 

0.4598E  07 

41 

EC12 

0.6164E  06 

42 

NUC12 

Q.8126E-01 

43 

NUC21 

0 . 3084E-01 

44 

CSN13 

-0.3615E-06 

45 

CSN31 

-0. 1839E-07 

46 

CSN23 

0 .6574E-05 

47 

CSN32 

0 .8811E-06 

48 

CTE11 

0 .1102E-05 

49 

CTE22 

0 . 5805E-05 

50 

CTE12 

-0 .4330E-10 

51 

HK11 

0.2138E  03 

52 

HK22 

Q . 4755 E  02 

53 

HK12 

-0 .5785E-04 

54 

HHC 

0.1975E  00 

55 

DPCU 

0 .3573E-02 

56 

DPC22 

0 . 3578E-02 

57 

DPC33 

0.2000E  05 

58 

DPCU 

0.0000 

59 

BTACU 

0  . 1009E-03 

60 

BTAC22 

0 .3370E-03 

61 

BTAC33 

0.  1859E-02 

62 

BTACU 

-0 .2425E-08 

-DM- 

1009E-03 
.3370E-03 
.  1859E-02 
0.0000 
0  .  0  000 
.  2425E-08 


06 

06 
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Item  8 

T-FORCES 


0  .0  00  0 
0.0000 
0.0000 


0.0000 

0.0000 

0.0000 


Item  10 


SOME  USEFUL  DATA  FOR  F.E.  ANALYSIS 
COMPOSITE  THICKNESS  FOR  F.E.  ANALYSIS  =  G.30000E-01 

PROPERTIES  FOR  F.E.  ANALYSIS  Ell , E12, E13, E22, E23, E33  PROPERTIES  SCALED  BY  10**-6 
0.32533E-01  -0.67069E-02  -0.29839E-07  0.21747E  00  0.14296E-05  0.16224E  01 

BENDING  EQUIVALENT  PROPERTIES  NUCXY,  NCYX,  ECXX,  ECYY ,  GCXY 
0.22261E  00  0.20153E-01  0.16708E  08  0 . 15126  E  07  0.55473E  06 

NASTRAN  MEMBRANE  EQUIVALENT  ELASTIC  COEFFICIENTS  011,012,013,322,023,033 
0.12147E  08  0.37462E  06  0.10669E  00  0.46099E  07  0.40551E  01  0.61637E  06 

NASTRAN  BENDING  EQUIVALENT  ELASTIC  COEFFICIENTS  G1 1 , 012 , G1 3 , 022 , G23 , 033 
0.16734E  08  0.33824E  06  0.11854E-01  0.15194E  07  0.45056E  00  0.55473E  06 


FORCES 


FORCE  DISPLACEMENT  RELATIONS 


DISPL 


NX 

0.3644E  .06 

0.1124E  05 

0.3201E-02 

NY 

0.1124E  05 

0.1383E  06 

0.1217E  00 

NXY 

0 . 3201E-02 

0.1217E  00 

0.1849E  05 

-0 

MX 

0.0000 

O.CCCO 

-0.7276E-1I 

0 

MY 

0.0000 

-0 .1373E-03 

-0.26 1 9E-0  9 

0 

MXY 

-Q.7276E-U 

-0.2619E-09 

0.0000 

0 

REDUCED  STIFFNE5S  MATRIX 

0  . 36441E  06  0.11238E  05  0.32006E-02 
0.11238c  05  0 . 13330  E  06  0.12165E  00 
0. 32006E-02  0.12165E  00  0.18491E  05 


0  .0000 
0  .  0000 


00 


0.0000 

-0.1373E-03 

-0.2619E-09 


0.7610E  00 
0.3419E  01 
0.1014E-05 


-0.7276E-11  UX 
-0.2619E-09  V Y 

0.0000  VXPUY 


0.2667E-07  UXX 
0.1014E-05  WYY 
0.1248E  01  WXY 


Item  9 

REDUCED  BENDING  REGIDITIES 

0.37763E  02  0.76104E  00  0.26672E-07 

0.76104E  00  0 . 34186 E  01  Q.10138E-05 

0.26672E-07  0.10138E-05  0.12481E  01 


Item  11 

DISPLACEMENT  FORCE  RELATIONS  COMBINED  FORCES 


1 

Q.2751E-02 

-1- 

0 .2751E-05 

-2“ 

-0 . 2236  E“06 

-3“ 

0.9946E-I2 

-4- 

0.1818E-12 

-5- 

-0.9021E-11 

-6- 

-0.2356E-16 

0.1000E  04 

2 

“0 .2236E-03 

-0 .2236E-06 

0.7249E-05 

-0.4765E-10 

-0.5895E-11 

0.2925E-09 

0 . 128  3E“ 1 4 

0.0000 

3 

0 .9946E-09 

0.9946E-12 

-0 .4765E-10 

0 .5408E-04 

“0 . 3466E-16 

0.2237E-14 

-0.1181E-19 

0  .  0  000 

4 

0.1818E-09 

0.1818E-12 

-0 .5895E-11 

-0.3466E-16 

0.2660E-01 

-0 . 5922E“02 

0.4241E-08 

0.0000 

5 

-0  .  9021E-08 

-0.9021E-11 

0.2925E-09 

0.2237E-14 

-0.5922E-02 

0.2938E  00 

-0 .2385E“06 

0.0000 

6 

-0.2356E-13 

-0 .2356E-16 

0 . 1283E-14 

“0 . 1181E-19 

0. 4241 E- 08 

“0 .2385E“06 

0.8012E  00 

0.0000 

NOTE:  THE  DISPLACEMENTS  ARE  REFERENCE  PLANE  MEMBRANE  STRAINS  (UX  ,  VY  ,  VXPUY)  AND  CURVATURES  (UXX  ,  WYY  ,  WXY) 


1-FORCES 


0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 
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Item  12 


PLY  HYGROTHERMOMECHANICAL  PROPERTIES/RESPONSE 


FOR  LOAD  CONDITIONS 


MEMBRANE 

LOADS  NBS(X,Y,XY-M) 

ARE 

1000  . 

0  . 

0 

BENDING 

LOADS  MBS( X, Y , XY-M) 

ARE 

0  . 

0. 

0 

QXZ , QYZ 

AND  APPLIED  PRESSURES 

ARE 

0. 

0. 

0 

NOTE  :  NO  MOISTURE  OR  TEMPERATURE 

LAYER  PROPERTIES,  RONS-PROPERTY ,  COLUMNS-LAYER 


PLY  NUMBER 
MATERIAL  SYSTEM 

ORIENTATION 

1 

AS— /IMLS 
AS— /IMLS 

0  .  0 

2 

SGLA/HMHS 
AS— /IMHS 

90  .  0 

3 

SGLA/HMHS 
AS— /IMHS 

90  .  0 

4 

AS— /IMLS 

AS— /IMLS 

0.0 

1 

KV 

0 .2000E-01 

0  .  10  0  0E-01 

0 . 1000E-01 

Q.2000E-01 

2 

KF 

Q.5500E  00 

0.5580E  00 

0.5580E  00 

0 .5500E  00 

3 

KFB 

0.5390E  00 

0.5524E  00 

0.5524E  00 

0.5390E  00 

4 

KM 

0 . 450  0E  00 

0.4420E  00 

0.4420E  00 

0.4500E  00 

5 

KMB 

0.4410E  00 

0.4376E  00 

0.4376E  00 

0.4410E  00 

6 

RHOL 

0 .5443E-01 

0  .  6  334E-01 

0 .6334E-01 

0 .5443E-Q1 

7 

TL 

0 . 1000E-01 

0 .5000E-02 

Q.5QG0E-Q2 

Q.1000E-01 

8 

DELTA 

0 .5850E-G4 

0 .5215E-04 

0 .5215E-04 

0 .5850E-04 

9 

ILDC 

0  .0  0  0  0 

0.0000 

0.0000 

0 .0000 

10 

ZB 

0 . 5000E-02 

0.1250E-01 

0 .1750E-01 

0 .2500E-01 

11 

ZGC 

-0 . 1000E-01 

-0 .2500E-02 

0.2500E-02 

0  . 1000E-01 

12 

THCS 

0  .0  0  00 

0.0  00  0 

0 .0000 

0  .0000 

13 

THLC 

0  .0  0  0  0 

0.1571E  01 

0 . 157  IE  01 

0  .  0000 

14 

THLS 

0  .0  00  0 

0.1571E  01 

0 . 157  IE  C  l 

0.0  000 

15 

sen 

0 . 210  IE  08 

0.1323E  08 

0.1323E  08 

0.21Q1E  08 

16 

SCI  2 

0.7797E  06 

0.8684E  06 

0.8684E  06 

0.7797E  06 

17 

SCI  3 

0.7797E  06 

0.8684E  06 

0.8684E  06 

0 .7797E  06 

18 

SC22 

0.1631E  07 

0.2166E  07 

0.2166E  07 

0.1631E  07 

19 

SC23 

0.8713E  06 

0.9537E  06 

0.9537E  06 

0 .8713E  06 

20 

SC33 

0.1776E  07 

0.2307E  07 

0.2307E  07 

0.1776E  07 

21 

SC44 

0.3238E  06 

0 . 456  IE  06 

0 . 456 IE  06 

0.3238E  06 

22 

SC55 

0.5470E  06 

0 . 7 55 IE  06 

0.7551E  06 

0.5470E  06 

23 

SC66 

0.5470E  06 

0.7551E  06 

0 . 7  55  IE  06 

0.5470E  06 

24 

CTE11 

0.1418E-06 

0.1603E-05 

0 . 1603E-05 

0 . 1418E-06 

25 

CTE22 

0 .2464E-04 

0. 1601 E- 04 

0 . 1601E-04 

0 .2464E-04 

26 

CTE33 

0 .2464E-04 

0 . 1601E-04 

0 . 1601E-04 

0 .2464E-04 

27 

HK1 1 

0.3195E  03 

0.1352E  03 

0.1352E  03 

0.3195E  03 

28 

HK22 

0.3702E  01 

0.2305E  01 

0 .2305E  01 

0.3702E  01 

“29 

HK33 

0.3702E  01 

0.2305E  01 

0.2305E  01 

0.3702E  01 

30 

HCL 

Q.1991E  00 

0.1943E  00 

0 . 1943E  00 

0 . 1991E  00 

31 

ELI  1 

0.1726E  08 

Q.1144E  08 

Q.1144E  08 

0.1726E  08 

32 

EL22 

0.1127E  07 

0.1696E  07 

0 . 1696E  07 

0.1127E  07 

33 

EL33 

0.1127E  07 

0.1696E  07 

0.1696E  07 

0.1127E  07 

34 

GL23 

0.3238E  06 

0 . 456 IE  06 

0 . 456  IE  06 

0.3238E  06 

35 

GL1 3 

0.5470E  06 

0 . 7 55 IE  06 

0.7551E  06 

0.5470E  06 

36 

GL1 2 

0.5470E  06 

0.7551E  06 

0.7551E  06 

Q.5470E  06 

137 

NUL12 

0.2945E  00 

0.2663E  00 

0.2663E  00 

0.2945E  00 

38 

NUL21 

0  .  1922E-01 

0 . 3947E-01 

0 . 3947E-01 

0 . 1 922E-01 

39 

NUL13 

0.2945E  00 

0.2663E  00 

0.2663E  00 

0.2945E  00 

40 

NUL31 

0  . 1922E-Q1 

0.3947E-01 

0.3947E-01 

0 . 1922E-01 

41 

NUL23 

0 . 482 IE  00 

0.3985E  00 

0.3985E  00 

0 . 482 IE  00 

42 

NUL32 

0.4821E  00 

0.3985E  00 

0.3985E  00 

0 . 48 2 IE  00 

43 

DPL1 

0 .8600E-04 

0.8736E-04 

0 .8736E-04 

0.8600E-04 

44 

DPL2 

0 .5168E-Q4 

0.5117E-04 

0 .5117E-04 

0 .5168E-04 

45 

DPL3 

0 .5168E-04 

0.5117E-04 

0 .5117E-04 

0.5168E-04 

46 

BTAL1 

0  .4981E-04 

0 . 9858E-04 

0 . 9853E-04 

0 .4931E-04 

47 

BTAL2 

0  . 1452E-02 

0 .8565E-03 

0 .8565E-03 

0 .1452E-02 

48 

BTAL3 

0.1452E-02 

0.1455E-02 

0 . 1455E-02 

0 .1452E-02 

49 

ILMFC 

0  .0  00  0 

0.8405E  02 

0.8916E  02 

0.8405E  02 

50 

TEMPD 

0  .  0  00  0 

0.0000 

0  .0000 

0  .0000 

51 

LSC11T 

0.2228E  06 

0.2168E  06 

0.2168E  06 

0.2228E  06 

52 

LSC11C 

0.8764E  05 

0.1665E  06 

0.1665E  06 

0.8764E  05 

53 

LSC11D 

0.8764E  05 

0.1665E  06 

0.1665E  06 

0.8764E  05 

54 

LSC22T 

0.5006E  04 

0.9915E  04 

0.9915E  04 

0.5006E  04 

55 

LSC22C 

0.1502E  05 

0.2314E  05 

0.2314E  05 

0.15O2E  05 

56 

LSC12 

0.5126E  04 

0.1195E  05 

0.1195E  05 

0.5126E  04 

57 

LSC23 

0.3983E  04 

0.1019E  05 

Q.1019E  05 

0.3983E  04 

58 

LSCC23 

0  .00  0  0 

0.6196E  05 

Q.1351E  06 

0.4417E  05 

59 

LSCC13 

0  .0  0  00 

0.7303E  05 

0.8147E  05 

0.5238E  05 

60 

LSCDF 

0  .0  0  00 

0 .4164E-03 

0.3925E-03 

0 .4164E-03 

61 

KL12AB 

0.9858E  00 

0.9075E  00 

0.9075E  00 

0.985SE  00 

62 

MDEIE 

0.9646E  00 

0 , 78 0 IE  00 

0 . 78  0  IE  00 

0.9646E  00 

63 

RELP.OT 

0.0  00  0 

0.1000E  01 

0.1000E  01 

0.1000E  01 

64 

EPS11 

0 .2751E-Q2 

-0 .2236E-03 

-0.2236E-03 

0 .2751E-02 

65 

EPS22 

-0 .2236E-03 

0 .2751E-02 

0 . 275 IE- 0  2 

-0 . 2236E-0  3 

66 

EPS12 

0 . 9946E-09 

-0 .8536E-08 

-0.8536E-08 

0 . 9946E-09 

67 

SIG11 

0.4769E  05 

-0.1329E  04 

-0.1329E  04 

0.4769E  05 

68 

SIG22 

0.6647E  03 

0.4614E  04 

0.4614E  04 

0.6647E  03 

69 

SIG12 

0 .5441E-03 

-0 .6445E-02 

-0 .6445E-02 

0 .5441E-03 

70 

DELFI 

0 .0000 

-0.4765E-08 

0.0000 

0.4765E-08 

71 

HFC 

0.1121E  01 

0.6393E  00 

0.6393E  00 

0.1121E  01 

72 

MPCTGE 

0  .0  000 

0.0000 

0.0000 

0  .  0000 

73 

SIG13 

0.0  0  00 

0.0000 

0.0000 

0.0000 

74 

SIG23 

0.0  0  00 

0  .0000 

0  .  0000 

0  .0000 

75 

SIG33 

0  .0  000 

0  .  0000 

0  .  0000 

0  .  0000 
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Item  13 


D  E.T  A 

I  L  S  OF 

POISSON  RATIO 

M  I  S  M  A 

T  C  H 

FOISSON”S 

RATIOS  OF  THE  COMPOSITE 

ANUCXY 

0.0813 

ANUCYX 

0 .0308 

ANUCSX 

-0.0  000 

AMUCSY  = 

0  .00  00 

NO. 

THETA 

ANULXY 

ANULSX  ANULSY 

POIDFN 

POIDFS 

1 

0  .  0 

0.2945 

0  .  000  0  0.0  00  0 

0.2132 

-0  .  0  0  0  0 

2 

90  .  0 

0 .0  395 

-0.0000  -0.0000 

-0 .0418 

-0  .0  000 

3 

90.0 

0 .0395 

-0 . 0  000  -0  .  00  00 

-0  .  0418 

-0  .  0  000 

4 

0  .  0 

0.2945 

0.0000  0.0000 

0.2132 

-0.0000 

Item  14 


PLY 


1 

1 

2 

2 

3 

3 

4 
4 


FREE  EDGE  STRESSES 


THETA 

SIGXX 

SIGYY 

SIGXY 

YDCAY  LENGTH 

SIGZY 

SIGZZ 

SIGZX 

0  .  0 

0  .  143E 

01 

0. 199E-01 

0.163E-07 

0.119E  00 

0 . 116E-01 

0 .337E-02 

0 .412E-08 

0  .  0 

0  .  143E 

01 

0  . 199E-01 

0 . 163E-07 

0  .000 

0  .  116E-01 

0.337E-02 

0 . 412E-08 

90.0 

0  .  1  38E 

00 

-0 . 399E-01 

-0 .326E-07 

0.298E  00 

-0 . 979E-08 

0 .803E-03 

0 .248E-03 

90.0 

0  . 138E 

00 

-0 . 399E-01 

-0 .326E-07 

0.345E  00 

-0 .845E-08 

0 .200E-03 

0 .248E-08 

90.0 

0  .  138E 

00 

-0 . 399E-01 

-0 .326E-07 

0.298E  00 

-0 .574E-08 

0 . 803E-03 

0 . 248E-08 

90.0 

0  .  138E 

00 

-0.399E-01 

-0 .326E-07 

0.345E  00 

-0 .495E-08 

0 .200E-03 

0 .248E-Q8 

0  .  0 

0 . 143E 

01 

0  . 199E-01 

0.163E-07 

0.119E  00 

0 .116E-01 

0 .337E-02 

0 . 412E-08 

0  .  0 

0.143E 

01 

0.199E-01 

0  . 163E-07 

0.000 

0.116E-01 

0 .337E-02 

0 . 412E-08 

NOTE: 

NOTE: 

NOTE: 


THE  INTERLAMINAR  STRESSES  ARE  BETWEEN  PLIES  (1-1)  AND  (I). 

IF  THE  PLY  HO  IS  REPEATED  THEN  THE  SECOND  ONE  INDICATES  STRESSES  IN  THE  SECONDARY 
FOR  ANGLE  PLY  LAMINATES  SIGYY  IS  0.  CONSEQUENTLY  SIGZY  AN  D  SIGZZ  ARE  COMPUTED  AS 
TO  OBTAIN  NONTRIVIAL  SIGZY  AND  SIGZZ,  ONE  MUST  SPECIFY  A  THIN  INTERPLY  LAYER. 

THE  INTERPLY  LAYER  THICKNESS  MAY  BE  OBTAINED  FROM  THE  PLY  PROPERTY  TABLE. 


COMPOSITE. 

ZERO. 


Item  15(a) 

microstresses 


FOR  LOAD  CONDITIONS 
MEMBRANE  LOADS  NESCX,Y,XY-M)  ARE 

BENDING  LOADS  MES(X, Y,XY-M)  ARE 

QXZ,QYZ  AND  APPLIED  PRESSURES  ARE 
NOTE  •  NO  MOISTURE  OR  TEMPERATURE 


(NOTE:  R0W5-PR0PERTY »  COLUMNS-LAYER) 


1000  . 
0. 
0  . 


0. 
0  . 
0  . 


0. 
0  . 
0. 


0. 


MAlERmESYSTEM  AS-/ks  SG  LA/HMHS 

AS— /IMLS  AS--/IMHS 

ORIENTATION  0.0  90.0 


1 

SMIL 

1 

SMIL 

2 

SM1T 

2 

SM1T 

3 

SF1L 

3 

SF1L 

4 

SF1T 

4 

SF1T 

5 

SM2AL 

5 

SM2AL 

6 

SM2AT 

6 

SM2AT 

7 

SM2BL 

7 

SM2BL 

8 

SM2BT 

8 

SM23T 

9 

SF2BL 

9 

SF2BL 

10 

SF2BT 

10 

SF2BT 

11 

SM3AL 

11 

SM3AL 

12 

SM3AT 

12 

SM3AT 

13 

SM3BL 

13 

SM3BL 

14 

SM3BT 

14 

SM3BT 

15 

SF3BL 

15 

SF3BL 

16 

SF3BT 

16 

SF3BT 

17 

SM12A 

17 

SM12A 

18 

SM12B 

18 

SM12B 

19 

SF12B 

19 

SF12B 

20 

SM13A 

20 

SM13A 

21 

SM13B 

21 

SM13B 

22 

SF13B 

22 

SFI33 

23 

SM23A 

23 

SM23A 

24 

SM23B 

24 

SM23B 

25 

SF23B 

25 

SF23B 

0.1381E  04 
0.0000 
0.2669E  03 
0.0000 
0.8564E  05 
0.0000 
-0.2185E  03 
0.0000 
0.1595E  03 
0 .0000 
0.3445E  03 
0.0000 
-0.1662E  05 
0.0000 
0.7763E  03 
0.0000 
-0.I662E  05 
0  .0000 
0.7763E  03 
0.0000 
0.1595E  03 
0.0000 
-0.2126E  02 
0.0000 
-0.1662E  05 
0.0000 
0.23X6E  03 
0.0000 
-0.1662E  05 
0  .  0000 
0.23I6E  03 
0.0000 
0.2137E-03 
0.0000 
0.6592E-03 
0.0000 
Q.6592E-03 
0.0000 
0.2137E-03 
0.0000 
0.6592E-03 
0.0000 
0.6592E-03 
0  .  000  0 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0  .  000  0 


-0.8714E  02 
-0.5809E  02 
0.1534E  04 
Q.1561E  04 
-0.1441E  04 
-0.3602E  04 
-0 . 4087 E  03 
-0.2406E  04 
-0 . 7294E  01 
-0.4862E  01 
0.1706E  04 
0 . 2324 E  04 
0.4632E  03 
0.4314E  03 
0.3858E  04 
0.7882E  04 
0.4632E  03 
0.4314E  03 
0.3858E  04 
0.7882E  04 
-0.7294E  01 
-0.4862E  01 
-0.9887E  02 
-0.6591E  02 
0.4632E  03 
0.4314E  03 
0.1607E  04 
0.1497E  04 
0.4632E  03 
0.4314E  03 
0  . 1607E  04 
0 . 14 97 E  04 
-Q.2348E-02 
-0.2439E-02 
-0 . 6  382E-02 
-0. 9945 E- 02 
-0.6382E-02 
-0.9945E-02 
-0 . 2348E-02 
-0.2439E-02 
-0.6382E-02 
-Q.9945E-02 
-0 .6382E-02 
-0 . 9945E-Q2 
0  .  0000 
0 .0000 
0 .0000 
0 .0000 
0.0000 
0.0000 


3 

SGLA/HMHS 
AS--/IMHS 
90  .  0 


-0.8714E  02 
-0.5809E  02 
0.1534E  04 
0.1561E  04 
-0.1441E  04 
-0.3602E  04 
-Q.4087E  03 
-0.24Q6E  04 
-0.7294E  01 
-0.4862E  01 
0.1706E  04 
0.2324E  04 
Q.4632E  03 
0.4314E  03 
0.3858E  04 
0.7882E  04 
0.4632E  03 
0.4314E  03 
0.3358E  04 
0.7882E  04 
-0.7294E  01 
-0.4862E  01 
-0 . 9887 E  02 
-0 . 6591 E  02 
0.4632E  03 
0.4314E  03 
0.1607E  04 
0.1497E  04 
0.4632E  03 
0.4314E  03 
0.1607E  04 
0 . 14 97  E  04 
-0 .2348E-02 
-0. 2439 E- 02 
-0 . 6  382E- 02 
-0 . 9945E-Q2 
-0.6382E-02 
-0.9945E-02 
-0.2348E-02 
-0.2439E-02 
-0.6382E-02 
-0 . 9945E-02 
-0.6382E-02 
-0  .  9945E-02 
0  .  00  0  0 
0  .0  0  0  0 
0  .  0  0  0  0 
0  .0  00  0 
0  .  0  0  0  0 
0.0000 


4 

AS— /IMIS 
AS— /IMLS 
0  .  0 

0.1381E  04 
0.0000 
0.2669E  03 
0.0000 
0.8564E  05 
0.0000 
-0 . 2185 E  03 
0.0000 
0.1595E  03 
0.0000 
0.3445E  03 
0.0000 
-0.1662E  05 
0.0000 
0.7763E  03 
0.0000 
-0.1662E  05 
0  .0  0  0  0 
0.7763E  03 
0  .  0  0  0  0 
0.1595E  03 
0  .0  0  0  0 
-0.2126E  02 
0  .  0  000 
-Q.1662E  05 
0.0000 
0.2316E  03 
0  .  0000 
-0.1662E  05 
0  .  0  0  0  0 
0.2316E  03 
0 .0000 
0 .2137E-03 
0.0000 
0  . 6592E-03 
0.0000 
0.6592E-03 
0.0000 
0.2137E-03 
0.0000 
0.6592E-03 
0 .0000 
0 .6592E-03 
0  .0  0  0  0 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 


NOTATION:  S  —  STRESS  (SIGMA) 

M  -  MATRIX  AND  F  -  FIBER  ,VCQ 

1,2,3  —  DIRECTIONS  FOR  STRESSES  -  PLY  MATERIAL  AXES 

L,T  -  DIRECTIONS  OF  PLY  STRESSES 

A  — -  REGION  CONTAINING  NO  FIBERS 

B  -  REGION  CONTAINING  FIBERS  AND  MATRIX 

EXAMPLE*.  SM2AL  STANDS  FOR  TRANSVERSE  NORMAL  STRESS 

IN  REGION  A  DUE  TO  A  LOAD  IN  THE  LOGITUDINAL 
DIRECTION 


Item  15(b) 

MICROSTRESS  INFLUENCE  COEFFICIENTS 


THE  FOLLOWING  ARE  THE  MICROSTRESS  INFLUENCE  COEFFICIENTS  FOR  THE  PRIMARY  COMPOSITE  AS— /IMLS  SYSTEM 


INF.  COEF. 

SIGMA  1 1 
LBS/SQ. IN 

SIGMA22 
LtfS/SQ .  IN 

5IGMA12 

LBS/SQ.  IN 

SIGMA13 
LBS/SQ. IN 

SIGMA23 
LBS/SQ. IN 

DELTA  T 

1  DEG  F 

DELTA  M 
IX 

1  SMI  1 

0  .  0290 

0 .4015 

0.0000 

0.0000 

0  .000  0 

-28.4291 

-1975.0911 

2  SM22A 

0.0033 

0.5183 

0.0000 

0.0000 

0.0000 

-16.1814 

-1274.0061 

3  SM22B 

-0.3484 

1.1678 

o.ocoo 

0  .00  00 

0  .0000 

5.6376 

443.8633 

4  5M12A 

0.0000 

0.0000 

0.3927 

0.0000 

0  .0000 

0.0000 

0  .0000 

5  SM12B 

0 .0000 

0.0000 

1.2116 

0.0000 

0.0000 

0  .0  000 

0  .00  00 

6  SMI3A 

0.0000 

0.0000 

0.0000 

0.3927 

0  .0000 

0.0000 

0.0000 

7  SM13B 

0.0000 

0.0000 

0.0000 

1.2116 

0.0000 

0  .  0  00  0 

0.0000 

3  SM23A 

0  .0000 

0  .  ocoo 

0 .0000 

0.0000 

0.5475 

0.0000 

0.0000 

9  SM23B 

0.0000 

0  .  ocoo 

0 .0000 

0.0000 

1.4043 

0  .0000 

0  .0  0  0  0 

10  SM33A 

0 .0033 

-0.0320 

0 .0000 

0.0000 

0.0000 

-16.1814 

-1274.0061 

11  SM33B 

-0 . 3484 

0.3484 

0 .0000 

0.0000 

0.0000 

5.6376 

443.8633 

12  SF1 1 

1.7955 

-0.3288 

0.0000 

0.0000 

0.0000 

21.4465 

1544.1631 

13  SF22B 

-0.3484 

1.1678 

0.0000 

0.0000 

0.0000 

5.6376 

443.8633 

14  SF33B 

-0.3484 

0.3484 

0 .0000 

0.0000 

0.0000 

5.6376 

443.8633 

15  SF12 

0  .0000 

0.0000 

1.2116 

0.0000 

0.0000 

0  .0  000 

0.0000 

16  SF1 3 

0 .0000 

0.0000 

0.0000 

1.2116 

0.0000 

0.0000 

0.0000 

17  SF23B 

0.0000 

0.0000 

0.0000 

0.0000 

1.4043 

0.0000 

0.0000 

NOTE:  TO  OBTAIN 

STRESSE5  OR  THE 

THE  ABSOLUTE  VALUE  OF  THE  MICROSTRESSES  THE 
TEMPERATURE  GRADIENT  OR  THE  MOISTURE  CONTENT. 

INF.  COEF. 

.  SHOULD  BE  MULTIPLIED  BY  THE 

APPROPRIATE 

EXPLANATION:  SM22B 

, FOR  EXAMPLE, 

STANDS  FOR  TRANSVERSE  NORMAL 

STRESS  INFLUENCE 

COEFFICIENT 

IN  REGION  B 

THE  FOLLOWING 

MICROSTRESS  INFLUENCE  COE 

ARE  THE  MICROSTRESS  INFLUENCE  COEFFICIENTS  FOR  THE  PRIMARY 

FFICIENTS 

COMPOSITE  SGLA/HMHS  SYSTEM 

INF.  COEF. 

SIGMA11 
LBS/SQ. IN 

SIGMA 2 2 
LBS/SQ. IN 

SIGMA12 

LBS/SQ. IN 

SIGMA13 
LBS/SQ. IN 

SIGMA23 
LBS/SQ. IN 

DELTA  T 
1  DEG  F 

DELTA  M 
IX 

1  SMI  1 

0.0655 

0.3325 

0.0000 

0.0000 

0.0000 

-28.7981 

-2926.0596 

2  SM22A 

0.0055 

0.569S 

0  .0  0  0  0 

0.0000 

0.0000 

-17.9927 

-2357.6494 

3  SM22B 

-0 . 3484 

0.8363 

0.0030 

0.0000 

0  .  0000 

6.2687 

821.4043 

4  SMI 2 A 

0.0000 

0.0000 

0.3643 

0.0000 

0.0000 

0.0000 

0.00  0  0 

5  SMI 23 

0 .0000 

0  .0  0  0  0 

0.9902 

0.0000 

0.0000 

0 .0000 

0.00  0  0 

6  SM13A 

0.0000 

0.0000 

0 .0000 

0.3643 

0.0000 

0 .0000 

0  .0000 

7  SMI 33 

0. 0000 

0.0000 

0.0000 

0.9902 

0.0000 

0 .0000 

U  .  U  U  U  U 

8  SM23A 

0.0000 

0  .  0  0  0  0 

0.0000 

0.0000 

0.6091 

0.0000 

0  .  0  0  00 

9  SM23B 

0.0  0  0  0 

0.0000 

0.0000 

0.0000 

2.0423 

0  .0  000 

0.0000 

10  5  M  3  3  A 

0 .0055 

-0  .  0214 

0.0000 

0.0000 

0.0000 

-17  .  9927 

-2357.6494 

11  SM333 

-0 . 3484 

0 . 3484 

0.0000 

0.0000 

0.0000 

6.2687 

821.4043 

12  SF11 

1.0837 

-0.0386 

0.0000 

0.0000 

0.0000 

-14.8484 

1222.4417 

13  3F22B 

-0.3484 

0.8363 

0 .0000 

0.0000 

0.0000 

6.2687 

821.4043 

14  3F33B 

-0.3484 

0.3484 

0.0000 

0.0000 

0.0000 

6.2687 

821.4043 

15  S FI 2 

0.0000 

0.0000 

0.9902 

0.0000 

0.0000 

0.0000 

0.0000 

16  SF1 5 

0  .0  00  0 

0.0000 

0.0000 

0.9902 

0.0000 

0.0000 

0.0000 

17  SF23B 

0  .0  000 

0.0000 

0.0000 

0.0000 

2.0423 

0.0000 

0.0000 

NOTE:  TO  OBTAIN  THE  ABSOLUTE  VALUE  OF  THE  MICROSTRESSES  THE 
STRESSES  OR  THE  TEMPERATURE  GRADIENT  OR  THE  MOISTURE  CONTENT. 

INF.  COEF.  SHOULD  BE  MULTIPLIED  BY  THE 

APPROPRIATE 

EXPLANATION:  SM22B- 

, FOR  EXAMPLE, 

STANDS  FOR  TRANSVERSE  NORMAL  STRESS  INFLUENCE 

COEFFICIENT 

IN  REGION  B 

MICROSTRES 

S  I  N  F  L 

UENCE  COEF 

F  I  C  I  E 

NTS 

THE 

FOLLOWING 

ARE  THE  MICROSTRESS  INFLUENCE 

COEFFICIENTS 

FOR  THE  SECONDARY 

COMPOSITE 

AS—  /IMHS  SYSTEM 

INF. 

COEF. 

SIGMA11  SIGMA22 

SIGMA 12 

SIGMA13 

SIGMA23 

DELTA  T 

DELTA  M 

LBS/SQ. IN  LBS/SQ. IN 

LBS/SQ. IN 

LBS/SQ. IN 

LBS/SQ. IN 

1  DEG  F 

1% 

1 

SMI  1 

0 .0437 

0.3384 

0  .  0  0  0  0 

0.0000 

0.0000 

-17.1987 

-1950.7063 

2 

SM22A 

0.0037 

0.5036 

0  .0  00  0 

0.0000 

0.0000 

-9.9952 

-1571.7661 

3 

SM22B 

-0.3245 

1.7084 

0.0000 

0.0000 

0.0000 

3.2438 

510  .  0889 

4 

SMI  2  A 

0  .  0  0  00 

0 .0000 

0.3784 

0.0000 

0  .  0000 

0.0000 

0.0000 

5 

SM12B 

0  .  0  0  0  0 

0.0000 

1.5430 

0.0000 

0.0000 

0.0000 

0.0000 

6 

SM13A 

0  .  0  0  0  0 

0.0  00  0 

0.3784 

0.0000 

0.0000 

0.0000 

0.0000 

7 

SMI  33 

0.0  0  0  0 

0  .  000  0 

1.5430 

0.0000 

0.0000 

0  .0  0  0  0 

0.0000 

3 

SM23A 

0.0  0  0  0 

0  .  000  0 

0.0000 

0  .0  0  0  0 

0.4060 

0  .  0  0  0  0 

0.0000 

9 

SM23B 

0  .0  000 

0 .0000 

0  .0  0  0  0 

0.0000 

1.0551 

0.0000 

0.0000 

10 

SMS  3  A 

0 .0037 

-0.0143 

0  .  00  0  0 

0.0000 

0.0000 

-9.9952 

-1571.7661 

11 

SM333 

-0.3245 

0.3245 

0.0000 

0.0000 

0  .  0000 

3.2438 

510.0389 

12 

S  FI  1 

2.7093 

-0.5215 

0  .0  0  0  0 

0  .  0000 

0.0000 

66 . 7289 

3056.1045 

13 

SF22B 

-0.3245 

1.7084 

0  .0  0  0  0 

0  .00  0  0 

0.0000 

3.2438 

510.0889 

14 

SF33B 

-0.3245 

0.3245 

0.0000 

0.0000 

0  .0000 

3.2438 

510.0889 

15 

S  FI  2 

0  .0000 

0.0000 

1.5430 

0.0000 

0.0000 

0.0000 

0.0000 

16 

SF1 3 

0.0000 

0  .000  0 

1.5430 

0.0000 

0.0000 

0  .0  000 

0  .000  0 

17 

SF23B 

0  .0  00  0 

0.0000 

0.0000 

0.0000 

1.0551 

0.0000 

0.0000 

NOTE:  TO  OBTAIN  THE  ABSOLUTE  VALUE  OF  THE  MICROSTRESSES  THE  INF.  COEF.  SHOULD  BE  MULTIPLIED  BY  THE  APPROPRIATE 

STRESSES  OR  THE  TEMPERATURE  GRADIENT  OR  THE  MOISTURE  CONTENT. 


EXPLANATION:  SM22B  ,  FOR  EXAMPLE,  STANDS  FOR  TRANSVERSE  NORMAL  STRESS  INFLUENCE  COEFFICIENT  IN  REGION  B 
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Item  16 

STRESS  CONCENTRATION  FACTORS 
(AROUND  A  CIRCULAR  HOLE) 


NOTE:  K1XX  — >  STRESS  CONCENTRATION  FACTOR  DUE  TO  SIGMA  XX 

KIYY  — >  STRESS  CONCENTRATION  FACTOR  DUE  TO  SIGMA  YY 

K1XY  — >  STRESS  CONCENTRATION  FACTOR  DUE  TO  SIGMA  XY 


LAYUP  — > 

0  90  90 

0 

THETA 

K1XX 

KIYY 

K1XY 

THETA 

K1XX 

KIYY 

K1XY 

0.0 

-0.6160 

3.8562 

0.0000 

180.0 

-0.6160 

3.8562 

0.0002 

5.0 

-0.5709 

3.6729 

-1.1975 

185.0 

-0.5709 

3.6729 

-1.1973 

10.0 

-0.4572 

3.2156 

-2.1209 

190.0 

-0.4572 

3.2155 

-2.1208 

15.0 

-0.3168 

2.6650 

-2.6886 

195.0 

-0.3168 

2.6649 

-2.6885 

20.0 

-0.1799 

2.1516 

-2.9777 

200.0 

-0.1799 

2.1515 

-2.9777 

25.0 

-0.0569 

1.7253 

-3.1020 

205.0 

-0.0570 

1.7252 

-3.1020 

30.0 

0.0532 

1.3875 

-3.1493 

210.0 

0.0532 

1.3875 

-3.1493 

35.0 

0.1566 

1.1225 

-3.1741 

215.0 

0.1566 

1.1225 

-3.1741 

40.0 

0.2613 

0.9117 

-3.2080 

220.0 

0.2613 

0.9116 

-3.2080 

45.0 

0.3764 

0.7383 

-3.2701 

225.0 

0.3764 

0.7382 

-3.2701 

50.0 

0.5138 

0.5879 

-3.3730 

230.0 

0.5137 

0.5879 

-3.3730 

55.0 

0.6898 

0.4473 

-3.5263 

235.0 

0.6898 

0.4472 

-3.5262 

60.0 

0.9302 

0.3019 

-3.7367 

240.0 

0.9301 

0.3018 

-3.7367 

65.0 

1.2764 

0.1331 

-4.0029 

245.0 

1.2763 

0.1331 

-4.0028 

70 . 0 

1.7980 

-0.0864 

-4.2958 

250.0 

1.7979 

-0.0865 

-4.2957 

75.0 

2.6025 

-0.3968 

-4.5006 

255.0 

2.6023 

-0.3969 

-4.5006 

30 . 0 

3.7983 

-0.8378 

-4.2743 

260.0 

3.7980 

-0.8379 

-4.2744 

85.0 

5.2286 

-1.3548 

-2.9009 

265.0 

5.2283 

-1.3549 

-2.9013 

90.0 

5.9765 

-1.6233 

-0.0001 

270.0 

5.9765 

-1.6233 

-0.0008 

95.0 

5.2287 

-1.3549 

2.9007 

275.0 

5.2290 

-1.3550 

2.9003 

100.0 

3.7984 

-0.8379 

4.2743 

280.0 

3.7987 

-0.8380 

4.2741 

105.0 

2.6026 

-0.3969 

4.5007 

285.0 

2.6028 

-0.3970 

4.5007 

110.0 

1.7980 

-0.0865 

4.2958 

290.0 

1.7982 

-0.0865 

4.2958 

115.0 

1.2764 

0.1330 

4.0029 

295.0 

1.2765 

0.1330 

4.0030 

120.0 

0.9302 

0.3018 

3.7367 

300.0 

0.9303 

0.3018 

3.7368 

125.0 

0.6899 

0.4472 

3.5263 

305.0 

0.6899 

0.4472 

3.5263 

130.0 

0.5138 

0.5879 

3.3730 

310.0 

0.5138 

0.5878 

3.3731 

135 . 0 

0.3764 

0.7382 

3.2702 

315.0 

0.3764 

0.7382 

3.2702 

140.0 

0.2613 

0.9116 

3.2081 

320.0 

0.2613 

0.9116 

3.2081 

145.0 

0.1567 

1.1225 

3.1741 

325.0 

0.1567 

1.1224 

3.1741 

150.0 

0.0532 

1.3875 

3.1493 

330.0 

0.0532 

1.3874 

3.1493 

155.0 

-0.0569 

1.7252 

3.1021 

335.0 

-0.0569 

1.7251 

3.1021 

160.0 

-0.1799 

2.1515 

2.9777 

340.0 

-0.1798 

2.1514 

2.9778 

165.0 

-0.3168 

2.6648 

2.6886 

345.0 

-0.3168 

2.6647 

2.6887 

170.0 

-0.4571 

3.2155 

2.1210 

350.0 

-0.4571 

3.2153 

2.1211 

175.0 

-0.5709 

3.6728 

1.1976 

355.0 

-0.5709 

3.6728 

1.1974 
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Item  17 


LOCATIONS  OF  PROBABLE  DELAMINATION 


RESULTS  FOR  PLY  NO.  1  ORIENTATION  0 

.0 

MATERIAL  AS- 

-IMLS  AS  — IMLS 

CRITERION 

RANGE 

VALUE 

LOCATION  # 

MAX 

OF 

K1XX# (NUCRT -NULRT) 

0.0  — 

90.0 

0.588 

75.0 

MAX 

OF 

KIYYK(NUCRT-NULRT) 

0.0  — 

90.0 

0.822 

0.0 

MAX 

OF 

K1XY#(NUCRT-NULRT) 

0.0-- 

90.0 

1.411 

60.0 

MAX 

OF 

K1XX#(NUCRT-NULRT) 

90  .  0  — 

180.0 

0.588 

105.0 

MAX 

OF 

K1YY#(NUCRT-NULRT) 

90.0  — 

180.0 

0.822 

180.0 

MAX 

OF 

K1XY#(NUCRT- NULRT) 

90. 0-- 

180.0 

1.411 

120 . 0 

MAX 

OF 

K1XX#( NUCRT-NULRT) 

180.0  — 

270.0 

0.588 

255.0 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

180.0  — 

270.0 

0.822 

180.0 

MAX 

OF 

K1XY#CNUCRT-NULRT) 

180.0  — 

270.0 

1.411 

240.0 

MAX 

OF 

K1XX# (NUCRT-NULRT) 

270.0  — 

0.0 

0.588 

285.0 

MAX 

OF 

K1YYKCNUCRT-NULRT) 

270.0  — 

0.0 

0.822 

0.0 

MAX 

OF 

K1XY# (NUCRT-NULRT) 

270.0  — 

0.0 

1.411 

300.0 

RESULTS  FOR  PLY  NO.  2  ORIENTATION  90. 

.0 

MATERIAL  SGLAHMHS  AS  — IMHS 

CRITERION 

RANGE 

VALUE 

LOCATION  # 

MAX 

OF 

KIXXX(NUCRT-NULRT) 

0.0  — 

90.0 

1.407 

90.0 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

0.0  — 

90.0 

1.087 

15.0 

MAX 

OF 

K1XY#( NUCRT-NULRT ) 

0.0  — 

90.0 

1.476 

25.0 

MAX 

OF 

K1XX# (NUCRT-NULRT ) 

90.0  — 

180.0 

1.407 

90.0 

MAX 

OF 

K1YY#( NUCRT-NULRT) 

90.0  — 

180.0 

1.087 

165.0 

MAX 

OF 

K1XY#(NUCRT-NULRT) 

90.0  — 

180.0 

1.476 

155.0 

MAX 

OF 

KIXXX(NUCRT-NULRT) 

180.0  — 

270.0 

1.4  07 

270.0 

MAX 

OF 

K1YY#(NUCRT-NULRT) 

180.0  — 

270.0 

1.087 

195.0 

MAX 

OF 

K1XY#( NUCRT-NULRT) 

180.0  — 

270.0 

1.476 

205.0 

MAX 

OF 

K1XX#( NUCRT-NULRT) 

270.0  — 

»  0.0 

1.407 

270.0 

MAX 

OF 

KlYY^C NUCRT-NULRT) 

270.0  — 

0.0 

1.087 

345.0 

MAX 

OF 

K1XY#(NUCRT-NULRT) 

270.0  — 

0.0 

1.476 

335.0 

RESULTS  FOR  PLY  NO.  3  ORIENTATION  90. 

0 

MATERIAL  SGLAHMHS  AS  — IMHS 

CRITERION 

RANGE 

VALUE 

LOCATION  # 

MAX 

OF 

KIXXK(NUCRT-NULRT) 

0  .  0  — 

90.0 

1.407 

90.0 

MAX 

OF 

K1YY#(NUCRT-NULRT) 

o.o— 

90.0 

1.087 

15.0 

MAX 

OF 

K1XY#(NUCRT-NULRT) 

0  .  o  — 

90.0 

1.476 

25.0 

MAX 

OF 

K1XX#( NUCRT-NULRT ) 

90.0  — 

180.0 

1.407 

90.0 

MAX 

OF 

KIYY# (NUCRT-NULRT) 

90  .  0  — 

180.0 

1.087 

165.0 

MAX 

OF 

K1XY#( NUCRT-NULRT) 

90 . 0 -- 

180.0 

1.476 

155.0 

MAX 

OF 

K 1XX# ( NUCRT-NULRT ) 

180.0  — 

270.0 

1.407 

270.0 

MAX 

OF 

KIYY# (NUCRT-NULRT) 

180.0  — 

270.0 

1.087 

195.0 

MAX 

OF 

K1XY#( NUCRT-NULRT) 

180.0  — 

270.0 

1.476 

205.0 

MAX 

OF 

K1 XX# (NUCRT-NULRT ) 

270.0  — 

0.0 

1.407 

270.0 

MAX 

OF 

KIYY#(NUCRT~NULRT) 

270.0  — 

0.0 

1.087 

345.0 

MAX 

OF 

K1XY#( NUCRT-NULRT ) 

270.0  — 

0.0 

1.476 

335.0 

RESULTS 

FOR  PLY  NO.  4  ORIENTATION  0. 

0 

MATERIAL  AS— 

■IMLS  AS— IMLS 

CRITERION 

RANGE 

VALUE 

LOCATION  # 

MAX 

OF 

K1XX# (NUCRT-NULRT) 

0.0  — 

90.0 

0.588 

75.0 

MAX 

OF 

KIYY# (NUCRT-NULRT) 

0.0  — 

90.0 

0.822 

0.0 

MAX 

OF 

K1XY#( NUCRT-NULRT) 

0.0  — 

90.0 

1.411 

60.0 

MAX 

OF 

K1 XX# (NUCRT-NULRT) 

90.0  — 

180.0 

0.588 

105.0 

MAX 

OF 

KIYY# (NUCRT-NULRT) 

90.0  — 

180.0 

0.822 

180.0 

MAX 

OF 

K1XY#( NUCRT-NULRT) 

90  .  0  — 

180.0 

1.411 

120.0 

MAX 

OF 

K1XX#( NUCRT-NULRT) 

180 . 0-- 

270.0 

0.583 

255.0 

MAX 

OF 

K1YY#(NUCRT-NULRT) 

180.0  — 

270.0 

0.822 

180.0 

MAX 

OF 

K1XY#( NUCRT-NULRT) 

180 . 0  — 

270.0 

1.411 

240 . 0 

MAX 

OF 

K1XX#( NUCRT-NULRT) 

270.0  — 

0.0 

0.538 

285.0 

MAX 

OF 

KIYY# (NUCRT-NULRT) 

27  0 . 0  — 

0.0 

0.822 

0.0 

MAX 

OF 

K1XY#( NUCRT-NULRT) 

270.0  — 

0.0 

1.411 

300.0 

NOTES:  K1XX  — >  STRESS  CONCENTRATION  FACTOR  DUE  TO  SIGMA  XX 

KIYY  — >  STRESS  CONCENTRATION  FACTOR  DUE  TO  SIGMA  YY 

K1XY  — >  STRESS  CONCENTRATION  FACTOR  DUE  TO  SIGMA  XY 

NULRT  — >  PLY  POISSON  RATIO  IN  R  AND  T  AXE5 

NUCRT  -->  COMPOSITE  POISSON  RATIO  IN  R  AND  T  AXES 
(  R  AND  T  ARE  THE  RADIAL  AND  THE  TANGENTIAL  DIRECTIONS) 
ONLY  5  DEG.  INTERVALS  ARE  CONSIDERED.  THE  ACTUAL  VALUE 
IS  EXPECTED  TO  BE  WITHIN  5  DEG.  OF  THE  PRINTED  RESULT. 


Item  18 


p 

L  Y  S 

TRESS 

AND  ST 

RAIN  I 

N  F  L  U  E 

N  C  E  CO 

E  F  F  I  C  I 

ENTS  A 

R  R  A  Y  S 

PLY  NO.  MATERIAL 

THETA 

RESPONSE 

NX 

NY 

NXY 

MX 

MY 

MXY 

DELTAT 

DELTAM 

SYSTEM 

(UNIT 

LOAD  ...  LB 

./INCH) 

(UNIT  MOMENT. . .LB. 

IN/INCH) 

(1  DEG  F)  (1  X ) 

1  AS— /IMLS 

0.0 

EPS11 

2.7511 

-0.2236 

0.0000 

266.0017 

-59.2159 

0.0000 

1.1801 

105.0192 

EPS22 

-0.2236 

7.2490 

-0 .0000 

-59.2160 

2938.3235 

-0.0024 

5.8164 

342.4932 

EPS12 

0.0000 

-0.0000 

54.0302 

0  .0000 

-0.0024 

8011.9414 

-0  .00  0  0 

-0.0023 

SIGll 

47.6932 

-1.4628 

0.0000 

4598.8906 

-47.6819 

-0.0001 

11.7478 

588.3518 

SIG22 

0.6647 

8.1391 

-0.0001 

21.6664 

3309.6074 

-0.0027 

-20 . 9790 

-1238.7244 

SIG12 

0.0000 

-0.0000 

29.5831 

0.0000 

-0.0013 

4382.7227 

-0 .0  0  00 

-0.0012 

2  SGLA/HMHS 

90.0 

AS— /IMHS 

EPS11 

-0.2236 

7.2490 

0.0000 

-14.8040 

734.5815 

0.0019 

5.8164 

342.4929 

EPS22 

2.7511 

-0.2236 

-0.0  001 

66.5005 

-14.8040 

-0.0025 

1.1801 

105.0192 

EPS12 

-0.0000 

0 .0001 

-54.0802 

-0.0002 

0.0025 

-2002.9873 

0.0001 

0.0029 

SIGH 

-1.3294 

83.7219 

0 .0  002 

-140.8320 

8487.6133 

0.0213 

41 . 9581 

2477.4478 

SIG22 

4.6136 

2.9257 

-0.0001 

107.2302 

309.9329 

-0.0034 

-23.4957 

-1176.7065 

SIG12 

-0.0000 

0.0001 

-40.8337 

-0.0002 

0.0019 

-1512.3735 

0  .  00  00 

0.0022 

3  SGLA/HMHS 

90  .  0 

AS— /IMHS 

EPS  1 1 

-0.2236 

7.2490 

0.0000 

14.8040 

-734.5796 

-0.0019 

5.8164 

342.4927 

EPS22 

2.7511 

-0.2236 

-0 .0001 

-66.5004 

14.8040 

0 .0025 

1.1801 

105.0192 

EPS12 

-0.0000 

0.0001 

-54.0802 

0.0002 

-0.0025 

2002.9846 

0.0001 

0 .0029 

SIGH 

-1.3294 

83.7219 

0.0002 

140.8316 

-8487.5898 

-0.0213 

41.9580 

2477.4451 

SIG22 

4.6136 

2.9257 

-0.0001 

-107.2300 

-309.9319 

0.0034 

-23.4957 

-1176.7065 

SIG12 

-0  .0000 

0.0001 

-40.8337 

0.0002 

-0.0019 

1512.3713 

0.0000 

0.0022 

4  AS--/ IMLS 

0.0 

EPS1I 

2.7511 

-0.2236 

0.0000 

-266.0017 

59.2159 

-0.0000 

1 . 1801 

105.0192 

EPS22 

-0 . 2236 

7.2490 

-0 . 0000 

59.2159 

-2938.3228 

0.0024 

5.8164 

342.4924 

EPS12 

0.0000 

-0  .00  0  0 

54.0802 

-0  .00  00 

0.0  024 

-8011.9414 

-0  .00  0  0 

-0.0023 

SIGH 

47.6932 

-1.4628 

0.0000 

-4598.8906 

47.6819 

0  .0001 

11.7478 

538.3518 

5IG22 

0.6647 

8.1391 

-0.0001 

-21.6664 

-3309.6067 

0.0027 

-20.9790 

-1238.7251 

SIGI2 

0 .0000 

-0.0000 

29.5831 

-0.0000 

0.0013 

-4382.7227 

-0.0000 

-0.0012 

NOTE J  STRAINS  ARE  IN  MICRO  INCH/INCH. 

STRE5SES  ARE  IN  POUNDS/INCH  SQ. 

EXPLANATION  OF  THE  INFLUENCE  COEFFICIENTS 

NX, NY  AND  NXY  ARE  UNIT  LOADS  IN  LB/INCH.  MX, MY  AND  MXY 
ARE  UNIT  MOMENTS  IN  LB. IN/INCH.  DELTAT  IS  A  UNIT  TEMP. 
DIFF.  AND  DELTAM  IS  A  UNIT  PERCENTAGE  OF  MOISTURE  CONTENT. 
TO  OBTAIN  RESPONSE  R  FOR  A  GENERAL  APPLIED  LOAD 
VECTOR  F  USE  THE  FOLLOWING  EQUATION: 

(R)  =  (  AINF  )  X  (F) 

NOTE:  R  IS  A  6X1  COLUMN  VECTOR  DEFINED  BY 

T 

(R)  =  (EPS11  EPS22  EPS12  SIGH  SIGI2  SIG12) 

F  IS  A  8X1  COLUMN  VECTOR  DEFINED  BY 
,  T 

(F)  =  (NX  NY  NXY  MX  MY  MXY  DELTAT  DELTAM) 

AINF  IS  A  (6X8)  MATRIX  CONTAINING  THE  INFLUENCE 
COEFFICIENTS  ARRAYS. 


PLY  S 

TRESS 

I  N  F  L  U  E 

N  C  E  CO 

E  F  F  I  C 

I  E  N  T  S 

ARRAYS 

PLY  NO.  MATERIAL 
SYSTEM 

THETA 

RESPONSE 

NX 

(UNIT 

SCALE 

NY 

LOAD  ...  LB 
FACTOR  = 
Cl/TC) 

NXY 

./INCH) 

33.333 

MX  MY  MXY 

(UNIT  MOMENT. . .LB. IN/INCH) 
SCALE  FACTOR  =  6666.668 

(6/TCKK2) 

DELTAT 
(1  DEG  F) 
13.354 

DELTAM 
(1  X) 
1223.138 

1  AS— /IMLS 

0  .  0 

SIGll 

1.4308 

-0.0439 

0 .0000 

0.6898 

-0.0072 

-0.0000 

0.8797 

0.4810 

SIG22 

0  .0199 

0.2442 

-0 .0000 

0 .0032 

0.4964 

-0 ,0000 

-1.5709 

-1.0127 

SIG12 

0  .  00  0  0 

-0  .  0000 

0.8875 

0.0000 

-0.0000 

0.6574 

-0  .  0  0  0  0 

-0 .0000 

2  SGLA/HMHS 

90.0 

AS  — /IMHS 

SIGll 

-0.0399 

2.5117 

0  .  0000 

-0.0211 

1.2731 

0.0000 

3.1419 

2.0255 

SIG22 

0.1384 

0.0878 

-0 .0000 

0.0161 

0.0465 

-0  .  0  0  0  0 

-1.7594 

-0  .  9620 

SIG12 

-0  .000  0 

0.0000 

-1.2250 

-0.0000 

0.0  0  0  0 

-0.2269 

0.0000 

0  .  0000 

3  SGLA/HMHS 

90.0 

AS— /IMHS 

SIGll 

-0.0399 

2.5117 

0.0000 

0.0211 

-1.2731 

-0.0000 

3.1419 

2.0255 

SIG22 

0  . 1384 

0  .  0878 

-0.0000 

-0.0161 

-0.0465 

0  .  0  0  00 

-1.7594 

-0 .9620 

SIG12 

-0  .  0000 

0  .  0  000 

-1.2250 

0.0000 

-0.0000 

0.2269 

0.0000 

0 .0000 

4  AS  — /IMLS 

0  .  0 

SIGll 

1.4308 

-0.0439 

0.0000 

-0.6898 

0.0072 

0  .  000  0 

0.8797 

0.4810 

SIG22 

0 .0199 

0.2442 

-0.0000 

-0.0032 

-0.4964 

0  .0000 

-1.5709 

-1.0127 

SIG12 

0  .  000  0 

-0.0000 

0.8875 

-0.0000 

0.0000 

-0.6574 

-0.0000 

-0.0000 

NOTE:  THE  MEMBRANE  STRESSES  ARE  NORMALIZED  W.R.T  THE 

AVERAGE  STRESS  DUE  TO  UNIT  LOAD  IN  AN  EQUIVALENT 
HOMOGENEOUS  SECTION.  THE  BENDING  STRESSES  ARE 
NORMALIZED  W.R.T  THE  MAXIMUM  STRESS  DUE  TO  UNIT 
MOMENT.  THE  TEMPERATURE  AND  MOISTURE  STRESSES  ARE 
NORMALIZED  W.R.T  THE  AVERAGE  STRESSES  DUE  TO  UNIT 
TEMPERATURE  DIFFERENCE  AND  UNIT  PERCENTAGE  OF  MOIS¬ 
TURE.  TO  OBTAIN  THE  ABSOLUTE  VALUES  OF  THE  STR¬ 
ESSES  THE  INFLUENCE  COEFFICIENTS  SHOULD  BE  MULTI¬ 
PLIED  BY  THE  INDICATED  SCALE  FACTORS.  THESE  5H0ULD 
BE  MULTIPLIED  BY  THE  CORRESPONDING  LOADS  TO  OBTAIN 
STRESSES  IN  THE  PLIES. 
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Item  19 


LAMINATE  FAILURE  STRESS  ANALYSIS 


LAYUP  — >  0  90  90  0 

LAMINATE  FAILURE  STRESSES  BASED  UPON  FIRST  PLY  FAILURE  CRITERIA  (NO  TEMPERATURE  OR  MOISTURE  STRESSES) 


PLY  NO. 

=  1 

THETA  =  0 . 

.  0  0  MATERIAL 

SYSTEM  =  AS— IMLS 

AS— IMLS 

LOADS 

SLUT 

222.7741 

KSI 

SLUG 

87.6392 

KSI 

SL22T 

5.0065 

KSI 

SL22C 

15.0194 

KSI 

SL12S 

5.1261 

KSI 

FAIL.  LOAD 

KSI 

MODE 

SCXXT  MIN  ( 

155.699 

-61.252 

251.063 

-753.188 

0.000) 

155.699 

SLUT 

SCXXC  MIN  ( 

-155.699 

61.252 

-251.063 

753.188 

0.000) 

61.252 

SL11C 

SCYYT  MIN  ( 

-5076.305 

1997.015 

20.504 

-61.511 

0.000  ) 

20.504 

SL22T 

SCYYC  MIN  ( 

5 076.305 

-1997.015 

-20.504 

61.511 

0.000  ) 

61.511 

SL22C 

SCXYS  MIN  C 

0.000 

0.000 

********** 

********** 

5.776  ) 

5.776 

SL12S 

LAMINATE  FAILURE 

STRESSES  BASED  UPON 

FIRST  PLY  FAILURE  CRITERIA 

(NO  TEMPERATURE  OR 

MOISTURE  STRESSES) 

PLY  NO. 

=  2 

THETA  =  90.00 

MATERIAL 

SYSTEM  =  SGLAHMHS  AS— IMHS 

LOADS 

SLUT 

216.8321 

KSI 

SL11C 

166.5112 

KSI 

SL22T 

9.9151 

KSI 

SL22C 

23.1353 

KSI 

SL12S 

11.9513 

KSI 

FAIL.  LOAD 

KSI 

MODE 

SCXXT 

MIN 

(  -5436.801 

4175.063 

71.638 

-167.155 

0.000) 

71.633 

SL22T 

SCXXC 

MIN 

(  5436.801 

-4175.063 

-71.638 

167.155 

0.000) 

167 . 155 

SL22C 

SCYYT 

MIN 

C  86.330 

-66.295 

112.967 

-263.559 

**********) 

86.330 

SLUT 

SCYYC 

MIN 

(  -86.330 

66.295 

-112.967 

263.589 

**********) 

66.295 

SL11C 

SCXYS 

MIN 

(********** 

********** 

********** 

********** 

9.756) 

9.756 

SL12S 

LAMINATE  FAILURE 

STRESSES  BASED  UPON 

FIRST  PLY  FAILURE  CRITERIA 

(NO  TEMPERATURE  OR 

MOISTURE  STRESSES) 

PLY  NO. 

=  3 

THETA  =  90.00 

MATERIAL 

SYSTEM  =  SGLAHMHS  AS— IMHS 

LOADS 

SLUT 

SL11C 

SL22T 

SL22C 

SL12S 

FAIL.  LOAD 

MODE 

216.8321 

166.5112 

9.9151 

23.1353 

11.9513 

KSI 

KSI 

KSI 

KSI 

KSI 

KSI 

SCXXT 

MIN 

(  -5436.801 

4175.063 

71.638 

-167.155 

0.000) 

71.638 

SL22T 

SCXXC 

MIN 

(  5436.801 

-4175.063 

-71.638 

167.155 

0.000) 

167.155 

SL22C 

SCYYT 

MIN 

(  86.330 

-66.295 

112.967 

-263.589 

**********) 

86.330 

SLUT 

SCYYC 

MIN 

(  -86.330 

66.295 

-112.967 

263.589 

**********) 

66.295 

SL11C 

SCXYS 

MIN 

(********** 

********** 

********** 

********** 

9.756) 

9.756 

SL12S 

NOTE: 

„**********„  inPLIES  -"NOT  APPLICABLE"- 

LAYUP  — >  0  90  90  0 

LAMINATE  FAILURE  STRESSES  BASED  UPON  FIRST  PLY  FAILURE  CRITERIA 

(NO  TEMPERATURE  OR 

MOISTURE  STRESSES) 

PLY  NO. 

=  4 

THETA 

0.00 

MATERIAL 

SYSTEM  =  AS- 

-IMLS  AS— IMLS 

LOADS 

SL11T 

SLUG 

SL22T 

SL22C 

SL12S 

FAIL.  LOAD 

MODE 

222.7741 

87.6392 

5.0065 

15.0194 

5.1261 

KSI 

KSI 

KSI 

KSI 

KSI 

KSI 

SCXXT 

MIN  ( 

155.699 

-61.252 

251.062 

-753.187 

0.000) 

155.699 

SL11T 

SCXXC 

MIN  ( 

-155.699 

61.252 

-251.062 

753.187 

0.000) 

61.252 

SL11C 

SCYYT 

MIN  ( 

-5076 .309 

1997 .017 

20.504 

-61.511 

0.000) 

20.504 

SL22T 

SCYYC 

MIN  ( 

5076.309 

-1997 .017 

-20.504 

61.511 

0.000) 

61.511 

SL22C 

SCXYS 

MIN  ( 

0.000 

0  .000 

********** 

********** 

5.776) 

5.776 

SL12S 
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Item  20 

SUMMARY 


LAMINATE  FAILURE  STRESS  ANALYSIS  -  (NO  TEMPERATURE  OR  MOISTURE  STRESSES) 
(BASED  UPON  FIRST  PLY  FAILURE) 


LOAD  TYPE 

STRESS  KSI 

FAILURE  MODE  PLY  NO. 

THETA 

MATERIAL 

SYSTEM 

SCXXT 

71.638 

SL22T 

3 

90  .  0 

SGLAHMHS 

AS— IMHS 

SCXXC 

61.252 

SL11C 

4 

0.0 

AS— IMLS 

AS— IMLS 

SCYYT 

20.504 

SL22T 

1 

0.0 

AS — IMLS 

AS— IMLS 

SCYYC 

61.511 

SL22C 

1 

0.0 

AS— IMLS 

AS— IMLS 

SCXYS 

5.776 

SL12S 

4 

0.0 

AS— IMLS 

AS — IMLS 

LAMINATE 

FAILURE  STRESS  ANALYSIS 

-  (NO  TEMPERATURE  OR 

MOISTURE  STRESSES) 

(BASED 

UPON  FIBER 

FAILURE) 

LOAD  TYPE 

STRESS  KSI 

FAILURE  MODE 

PLY  NO. 

THETA 

MATERIAL 

SYSTEM 

SCXXT 

155.699 

SLUT 

4 

0.0 

AS— IMLS 

AS— IMLS 

SCXXC 

61.252 

SL11C 

4 

0  .  0 

AS— IMLS 

AS- -IMLS 

SCYYT 

86 . 330 

SLUT 

2 

90  .  0 

SGLAHMHS 

AS— II 1IIS 

SCYYC 

SCXYS 

66.295 

********* 

SL11C 

N/A 

2 

90  .  0 

SGLAHMHS 

AS- -IMHS 

NOTE:  IF  THERE  IS  NO  ANGLE  PLY  "SCXYS"  BASED  UPON  FIBRE  FAILURE  IS  NOT  PREDICTED. 
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Appendix  C 

Resident  Data  Bank  (FBMTDATA.BANK) 


T300 

FP  3000  0.300E-03  0.640E-01 

FE  0.320E  08  0.200E  07  0.200E  00  0.250E  00  0.130E  07  0.700E  06 

FT  -0.550E-06  0.560E-05  0.580E  03  0.580E  02  0.170E  00 

FS  0.350E  06  0.300E  06  0.000  0.000  0.000  0.000 

AS — 

FP  10000  0.300E-03  0.630E-01 

FE  0.310E  08  0.200E  07  0.200E  00  0.250E  00  0.200E  07  0.100E  07 

FT  -0.550E-06  0.560E-05  0.580E  03  0.580E  02  0.170E  00 

FS  0.400E  06  0.400E  06  0.000  0.000  0.000  0.000 

SGLA 

FP  204  0.360E-03  0.900E-01 

FE  0.124E  08  0.124E  08  0.200E  00  0.200E  00  0.517E  07  0.517E  07 

FT  0.280E-05  0.280E-05  0.750E  01  0.750E  01  0.170E  00 

FS  0.360E  06  0.300E  06  0.360E  06  0.300E  06  0.180E  06  0.180E  06 

HMSF  HIGH  MODULUS  SURFACE  TREATED  FIBER. 

FP  10000  0.300E-03  0.703E-01 

FE  0.550E  08  0.900E  06  0.200E  00  0.250E  00  0.110E  07  0.700E  06 

FT  -0.550E-06  0.560E-05  0.580E  03  0.580E  02  0.170E  00 

FS  0.280E  06  0.200E  06  0.000  0.000  0.000  0.000 

OVER  END  OF  FIBER  PROPERTIES. 

I MLS  INTERMEDIATE  MODULUS  LOW  STRENGTH  MATRIX. 

MP  0.460E-01 

ME  0.500E  06  0.410E  00  0.570E-04 

MT  0.125E  01  0.250E  00 

MS  0.700E  04  0.210E  05  0.700E  04  0.140E-01  0.420E-01  0.320E-01  0.320E-01 

MV  0.225E  00  0.420E  03 

IMHS  INTERMEDIATE  MODULUS  HIGH  STRENGTH  MATRIX. 

MP  0.440E-01 

ME  0.500E  06  0.350E  00  0.360E-04 

MT  0.125E  01  0.250E  00 

MS  0.150E  05  0.350E  05  0.130E  05  0.200E-01  0.500E-01  0.350E-01  0.350E-01 

MV  0.225E  00  0.420E  03 

HMHS  HIGH  MODULUS  HIGH  STRENGTH  MATRIX. 

MP  0.450E-01 

ME  0.750E  06  0.350E  00  0.400E-04 

MT  0.125E  01  0.250E  00 

MS  0.200E  05  0.500E  05  0.150E  05  0.200E-01  0.500E-01  0.400E-01  0.400E-01 

MV  0.225E  00  0.420E  03 

OVER  END  OF  MATRIX  PROPERTIES. 
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